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ABSTRACT 


The  use  of  a  digital  comb  filter  for  the  separation  of  two 
speakers  was  formulated  in  previous  efforts.  An  adaptive  pitch 
synchronous  filter  has  been  developed  as  an  alternative.  This 
development  stems  from  the  characteristics  and  structure  of  the 
speech  waveform.   In  the  comb  filtering  development,  there  were 
tradeoffs  between  desired  speaker  distortion  and  undesired  speaker 
separation.   In  the  development  of  the  adaptive  filtering  tech- 
niques, the  results  of  these  tradeoffs  will  be  examined  for 
possible  improvements  over  the  comb  filtering  techniques. 

A  series  of  tests  performed  on  test  signal  inputs  compares 
the  performance  of  both  the  comb  and  adaptive  f i  Itering  systems. 
Included  in  an  appendix  are  several  of  the  main  computer  programs 
used  in  the  computer  implementation.  A  series  of  listening  tests 
will  be  performed  with  these  systems  in  future  work  making  the 
appendix  a  necessary  part  of  the  thesis. 
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CHAPTER  I 

INTRODUCTION 

1.1   Introduction  and  Goals 

The  problem  of  separating  a  signal  waveform  from  a  noise  waveform 
has  been  the  topic  of  unlimited  research  for  many  years.   In  this 
discussion  the  problem  of  interest  is  one  of  separating  a  speech  wave- 
form from  a  "noise  waveform"  which  may  be  the  speech  waveform  from  a 
competing  speaker.  Two  methods  for  speech  enhancement  are  developed 
in  this  thesis:  comb  filtering  and  adaptive  filtering. 

The  remainder  of  Chapter  I  describes  the  speech  waveform  and  its 
characteristics  in  a  manner  that  will  be  helpful  for  terminology  and 
modeling.   In  Chapter  II  a  review  of  the  previous  systems  dealing  with 
this  specific  problem  will  be  covered,  and  the  problem  formulation 
introduction  for  this  thesis  will  be  outlined.  Chapter  III  deals  with 
the  algorithm  that  was  used  in  determining  the  fundamental  frequency 
of  a  speech  waveform  including  some  of  the  background  research  that  was 
developed  on  this  method.   In  Chapter  IV  the  adaptive  filtering  methods 
are  formulated  along  with  a  discussion  of  the  rules  used  and  problems 
that  should  be  encountered  with  the  system.  Chapter  V  concentrates  on 
the  test,  signal  formulation,  processing,  and  results  for  the  different 
filtering  systems  implemented.  This  chapter  also  deals  with  the 
comparison  methods  that  will  be  used  between  the  outputs  of  the 
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different  systems.  The  computer  implementation  of  the  systems  will 
be  discussed  in  Chapter  VI ,  and  Chapter  VII  describes  the  results  of 
the  actual  speech  waveforms  that  were  processed.   Included  in  this 
chapter  are  spectrograms  of  the  inputs  and  outputs  for  the  systems  used, 
The  overall  results  and  conclusions  are  covered  in  Chapters  VI  I  I  and  IX 
respectively.  The  Appendices  contain  computer  programs,  flowcharts, 
and  documentation  for  the  purpose  of  making  the  continuing  work  in  this 
problem  easier  for  those  persons  concerned. 

The  goals  of  this  thesis  can  be  summarized  in  the  following 

manner:  First,  the  work  that  had  previously  been  done  was  to  be 

duplicated.  This  involved  Implementing  the  computer  system  described 

1 
by  Vaden  Shields.   The  purpose  of  this  work  was  to  arrive  at  a 

starting  location  before  any  other  work  was  begun.  The  system  was  also 
implemented  for  the  purpose  of  comparison  with  any  systems  that  were 
generated  by  future  work.  Second,  a  pitch  synchronous  adaptive  filter- 
ing scheme  was  developed  and  implemented  on  the  computer.  Now  with 
these  two  completely  separate  systems,  their  results  could  be  compared 
to  give  the  third  goal.   Finally,  the  ideas  of  comb  filtering  and 
adaptive  filtering  were  examined  in  order  to  find  a  limit  to  their 
effectiveness  on  the  speech  enhancement  problem. 

1.2  The  Speech  Waveform  and  its  Characteristics 

Before  the  speech  waveform  can  be  processed  by  any  type  of  sys- 
tem, the  structure,  characteristics,  and  properties  of  the  waveform 
must  be  understood.  A  thorough  understanding  of  these  items  along 
with  the  definition  of  several  terms  is  necessary  before  the  various 
processing  systems  are  discussed. 
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Speech  production  takes  place  in  an  area  of  the  body  referred  to 
as  the  vocal  tract.  The  vocal  tract  can  be  best  described  as  an 
acoustical  tube  of  nonuniform  cross-sectional  area  that  originates 
at  the  vocal  cords  and  terminates  at  the  lips.  The  vocal  tract  of 
the  average  male  is  17  centimeters  long,  and  is  shorter  in  the  average 
female  adult.  The  nasal  tract  which  lies  between  the  velum  and  nostrils 
is  another  nonuniform  cross-sectional  area  tube  that  may  be  acoustically 
coupled  to  the  vocal  tract  by  the  opening  at  the  velum. 

Within  the  vocal  and  nasal  tracts  lie  the  major  components  for 

speech  production.  The  articulators,  composed  of  the  lips,  jaws,  tongue, 

2 
and  velum  vary  with  the  size  and  shape  of  the  tract.    The  velum  is 

also  referred  to  as  the  soft  palate  and  controls  the  acoustical  coupling 

between  the  vocal  and  nasal  tracts.  The  vocal  cords  I ie  at  the  lower 

end  of  the  vocal  tract.  The  vocal  cords  consist  of  a  pair  of  lips 

made  of  ligament  and  muscle,  and  the  opening  between  the  vocal  cords  is 

referred  to  as  the  glottis.  The  words  subglottal  and  supraglottal 

that  are  used  frequently  in  the  literature  refer  respectively  to  the 

areas  immmediately  below  and  above  the  glottal  opening. 

In  general  there  are  three  methods  of  producing  speech  waveforms. 
First  and  most  important,  voiced  sounds  are  produced  by  exciting  the 
vocal  tract  with  quasi-periodic  pulses  of  air  that  originate  in  the 
lungs  as  a  steady  flow  and  are  chopped  into  pulses  by  the  vocal  cord 
vibrations.  Second,  the  vocal  tract  is  constricted  in  one  area,  and 
then  by  means  of  forcing  the  air  from  the  lungs  through  the  constric- 
tion, a  turbulent  air  flow  is  created.  Fricative  sounds  are  produced 

in  this  manner.  Plosive  sounds  are  created  by  closing  the  vocal  tract 

3 
momentarily,  building  up  a  pressure,  and  then,  releasing  it. 
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Usually  sounds  that  use  the  vocal  cord  vibrations  for  excitation 
are  referred  to  as  voiced,  and  those  that  do  not  use  the  vocal  cords 
are  referred  to  as  unvoiced.  The  sounds  that  are  produced  with  the 
use  of  the  nasal  tract  are  named  nasal  sounds.  For  these  types  of 
sounds  the  nasal  tract  is  coupled  to  the  vocal  tract  at  the  velum. 

The  terms  pitch,  pitch  period,  and  fundamental  frequency  can  all 
be  defined  at  the  same  time.  These  terms  have  been  defined  in  many 
ways  by  different  authors.  G.  Fant  defines  them  in  them  in  the  follow- 
ing manner: 


The  basic  property  of  a  vocal  cord  sound  is  its  peri- 
odicity expressed  by  the  duration  T  of  a  complete 
voice  period  or  by  the  inverse  value  of  the  voice 
fundamental  frequency. 

Frt  =  1/Trt 
o     o 

Fundamental  pitch  and  fundamental  frequency  are  not 
synonymous,  but  these  terms  can  be  used  interchange- 
ably due  to  the  close  one-to-one  correspondence.   In 
more  strict  terminology  pitch  is  a  tonal  sensation 
and  frequency  a  property  of  the  sound  stimulus.  The 
duration  of  a  pitch  cycle  always  varies  somewhat 
from  one  period  to  the  next.  Such  variations  are 
systematic  determining  the  intonation  in  part,  acci- 
dental rather  than  unintentional,  but  nevertheless  of 
importance  for  the  naturalness  of  human  speech.  Only 
speaking  machines  are  capable  of  producing  a  perfect- 
ly monotonic  pitch. 4 


The  fundamental  frequency  of  the  vocal  cords  is  directly  proportional 
to  their  tension,  mass,  and  the  subglottal  air  pressure.  Since  most 
adult  males  have  vocal  cords  that  are  longer  and  thicker  than  adult 
females  the  fundamental  frequency  in  the  male  voice  is  lower  than  in 
the  female  voice. 


For  a  normally  speaking  voice  the  fundamental  frequency  can  range 

5 
between  60  hertz  and  400  hertz.  H.  L.  Shaffer  ,  through  a  series  of 
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experiments,  determined  that  the  average  fundamental  frequency  of  a 
male  speaker  is  125  hertz,  while  the  average  for  a  female  is  192  hertz. 

Another  important  term,  formant,  is  defined  as  the  natural  fre- 
quency of  the  vocal  tract  that  corresponds  to  a  resonance  or  peak  of 
energy  at  various  frequencies  referred  to  as  resonant  frequencies.  The 
resonant  or  formant  frequencies  depend  on  the  shape  of  the  vocal  tract 
and  the  positions  of  the  articulators. 


There  have  been  several  models  proposed  for  the  speech  production 

6 
apparatus.  The  model  proposed  by  Oppenheim  and  Schafer   gives  a  good 

insight  into  the  manner  in  which  the  speech  waveform  is  produced. 
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The  di sere te- time  model  in  figure  1-1  can  be  examined  in  several 
sections.  The  time  varying  digital  filter  corresponds  to  the  vocal 
tract,  and  the  coefficients  of  the  filter  may  be  changed  to  correspond 
to  the  frequency  response  of  the  time  varying  characteristics  of  the 
vocal  tract.  The  vocal  tract  is  usually  slowly  varying  in  normal  speech 
and  can  be  considered  fixed  for  short  periods  of  time  (on  the  order  of 
10  milliseconds).  The  filter  can  be  excited  by  a  train  of  impulses, 
p(n),  or  by  samples  from  a  random  number  generator,  r(n).   In  voiced 
speech  the  impulse  train,  p(n),  is  used  to  excite  the  filter,  and  the 
spacing  between  the  samples  corresponds  to  the  pitch  period.  For  un- 
voiced speech  the  sequence  r(n)  is  used  as  a  noise-type  source  for 
f i Iter  excitation. 

8 
Oppenheim  and  Schafer   describe  a  factor  g(n)  that  accounts  for 

the  fact  that  the  actual  glottal  pulses  are  not  impulses.  This  addi- 
tional factor  can  be  used  with  the  other  two  components  to  describe  the 
speech  waveform,  s  (n),  for  voiced  speech. 

sv(n)  =  p(n)  *  vv  (n)  *  g(n)  (1.1) 

where:  *  denotes  convolution 

The  above  expression  can  also  be  expressed  in  its  frequency  domain 

representation  as: 

S/eJ03)      =     P(e>)   VvCeJw)   G(eJw)  (1.2) 

As  described  above  the  impulse  train,  p(n),  has  a  spacing  between 
pulses  Tp  corresponding  to  the  pitch  period.  P(eJa))  is  also  a  train  of 
impulses  separated  by  intervals  of  fp  =  1/Tp.  This  shows  that  the 
energy  of  the  speech  waveform  will  lie  in  narrow  bands  centered  about 
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the  fundamental  frequency  and  its  harmonics. 

Whenever  the  speech  is  unvoiced,  the  speech  waveform  can  be  des- 
cribed as: 

su(n)  =  r(n)  *  vu(n)  *  g(n)  (1.3) 

Alternatively  as  before: 


Su(eJw)  =  R(eJW)  Vu(eJW)  G(eJW)  (1.4) 


There  is  not  much  that  can  be  said  about  the  properties  of  r(n) 
and  R(eJw)  other  than  it  has  "noise-like"  characteristics.  Therefore, 
for  unvoiced  sounds  the  spectrum  of  the  signal,  su(n),  lacks  any  type 
of  harmonic  structure. 

There  are  certain  other  aspects  of  the  human  speech  system  that 
should  be  considered.  The  ear  which  is  considered  to  be  a  part  of  the 
overall  speech  system  is  insensitive  to  errors  in  phase  of  a  signal, 
while  on  the  other  hand,  the  ear  is  extremely  sensitive  to  errors  in 
the  pitch  epochs  or  pitch  periods.  Many  experiments  have  shown  that 
the  pitch  of  a  speech  waveform  provides  its  naturalness,  and  if  there 
are  any  errors  in  pitch  of  a  processed  signal,  it  becomes  immediately 
apparent. 

This  discussion  has  been  directed  toward  the  explanation  and 
definition  of  certain  aspects  of  the  speech  waveform  that  will  be  used 
and  exploited  in  the  following  chapters.   It  was  intended  to  show  how 
the  speech  waveforms  are  produced  and  how  these  waveforms  can  be 
modeled. 
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CHAPTER  I  I 

REVIEW  OF  THE  LITERATURE 

2.1  Shieldst  System 

Since  the  foundations  of  this  thesis  were  taken  from  another  thesis, 
some  of  the  fundamental  concepts  and  ideas  used  in  a  thesis  by  Vaden 
Shields  will  be  covered  briefly.  Shields  in  his  thesis,  "Separation  of 
Added  Speech  Signals  By  Digital  Comb  Filtering",  suggests  that  the 
harmonic  structure  of  speech  may  be  the  key  to  the  effective  removal  of 
additive  noise.  He  attacks  the  problem  in  the  following  manner.   First, 
he  looks  into  the  problem  of  choosing  an  optimum  unit  sample  response 
for  the  filter.  Second,  he  examines  the  effects  of  unvoiced  speech  and 
filter  interaction.   Finally,  he  looks  into  the  tracking  and  detection 
of  one  signal's  fundamental  frequency  from  a  combined  signal.  The  first 
two  areas  will  be  covered  in  this  summary. 

The  filtering  strategy  used  is  simple.  Since  voiced  speech  has 

its  spectral  energy  concentrated  in  narrow  harmonic  bands,  another 

speech  signal  with  the  same  type  of  spectrum  would  have  its  spectral 

energy  concentrated  in  narrow  harmonic  bands  that  would,  for  the  most 

9 
part,  not  overlap  those  bands  of  the  first  signal.  This  property 

suggested  the  use  of  a  comb  filter  adjusted  to  allow  the  first  signal 
to  pass  and  the  second  to  be  attenuated. 

The  comb  filter  must  have  several  properties  when  it  is  used  in 
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this  application.   First,  the  filter  must  be  time  variant.  This  pro- 
perty evolves  from  the  fact  that  the  pitch  period  of  the  speech  varies 
with  time,  and  the  passbands  of  the  filter  must  be  able  to  move. 
Second,  the  filter  can  only  act  over  a  limited  portion  of  the  speech 
waveform  at  one  time.  This  is  due  to  the  fact  that  the  approximation 
of  periodicity  of  voiced  speech  only  holds  for  a  short  time  segment. 

Also,  the  duration  of  the  impulse  response  of  the  filter  should  be 

10 
finite  and  short  compared  to  changes  in  the  speech  waveform. 

The  property  above  that  requires  the  impulse  response  of  the 

filter  be  finite  implies  that  the  digital  filter  be  nonrecursive  or 

FIR.  There  are  several  ways  to  design  an  FIR  filter  and  the  method 

Shields  uses  is  to  truncate  the  infinitely  long  unit  sample  response 

with  a  finite  duration  data  window.   If  a  filter  were  desired  that 

would  pass  only  certain  frequencies,  the  "teeth"  or  passbands  would 

need  to  be  very  narrow.  One  type  of  filter  for  this  application  would 

be  a  unit  sample  train  in  frequency  with  the  spacing  of  the  samples 

fp  =   1/MT  hertz   for  -  °°  <  f  <  °°     where  M   is  an    integer  and  T   is  the 

11 
sampling  rate.    The  unit  sample  response  of  this  filter  is  a  unit 

sample  train  in  time  with  each  sample  separated  by  (M  -  1)  zeros  for 

12 
_  oo  <  n  <  oo.   This  filter  would  be  set  to  pass  only  harmonics  of  fp. 

As  mentioned  above,  the  infinitely  long  unit  sample  response  of 
the  filter  must  be  truncated  in  order  to  create  the  FIR  filter.  This 
is  done  by  multiplying  the  infinitely  long  unit  sample  response  of  the 
filter  by  the  unit  sample  response  of  the  data  window.  The  multipli- 
cation that  the  above  operation  implies  in  the  time  domain  corresponds 
to  a  convolution  of  the  respective  transforms  in  the  frequency  domain. 
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The  effect  of  this  operation  is  to  create  narrow  harmonically  spaced 
passbands  and  this  is  referred  to  as  a  digital  comb  filter. 

The  shape  of  the  passbands  and  the  stopbands  of  the  comb  filter 

13 
are  a  function  of  the  type  of  window  that  is  used.    Shields  explored 

the  effects  of  four  common  window  functions:  the  rectangular,  the 

Hanning,  the  Hamming,  and  the  Blackman. 

Another  way  of  looking  at  what  is  happening  in  this  process  is  to 
go  back  to  the  step  where  the  filter's  infinite  duration  unit  sample 
response  and  the  data  window  were  being  multiplied  together.   Looking 
at  these  results  gives  the  following  equations: 

Suppose  x(n)  is  a  finite  duration  sequence  with  a  Fourier  Trans- 
form, X(eJw),  where: 

X(eja))  =  j     xCn)  e-ja)n  (2.1) 

n=0 

Now,  suppose  that  x(n)  is  multiplied  by  a  function  f(n),  whose  unit 

sample  response  is  of  the  following  form: 


m-1     j27rrn 
f(n)  =  }_      I       e    m      ;  for  -  °°  <  n  <  °°      (2.2) 
m  r=0 


f(n)  is  a  train  of  unit  samples  that  have  (m-1)  zeros  between  each  unit 
sample.  If  z(n)  if  formed  from  the  product  of  x(n)  and  f(n),  then,  the 
resulting  Fourier  Transform  of  z(n)  is  calculated  as  follows: 

z(n)  =  x(n)  f(n)  (2.3) 

Substituting  into  equation  (2.1) 

N-1 
ZCeJW)  =  I       x(n)  f(n)  e"Jwn  (2.4) 

n=0 


-22- 


Substituting  equation  (2.2)  in  for  f(n)  in  (2.4)  gives: 


N-1 
Z(eJ'w)  =  I       x(n) 
n=0 


m-1     j2iTrn 

1  I       e    m 

m  r=0 


rJwn     (2.5) 


Performing  an  interchange  of  the  summations  in  (2.5)  gives: 


m-1  N-1         -j(o)  -  2ur)n 
Z(eJa))  =  1   I    I   x(n)  e         m 
m  r=0  n=0 


(2.6) 


Equation  (2.6)  may  also  be  written  in  the  following  form: 


2(eJ">  =  i  ^     X{e  j(i°  -  %?   } 
m  r=0 


(2.7) 


Interpreting  (2.7)  may  be  easier  if  considered  in  the  following  form: 

Z(eJw)  =  l{x(eJ°J)  +  X(eJ(w-^>)  +  ... 
m 

+  x(eJ(a)  -  2TT(m-1))}  J 

m      J      (2.8) 

From  equation  (2.8)  it  can  be  seen  that  if  X(eJw)  is  periodic  on  the 

interval  {0,2tt},  then,  Z(eJw)  is  composed  of  a  sequence  of  teeth  that 

are  periodic  on  {0 ,2jr_  },    and  of  the  same  form  as  X(eJw). 

m 

The  difference  equation  that  is  implemented  by  this  method  can  be 
expressed  as: 


K 
y(nT)  =  £    ak  x(nT  +  kmT) 
k=-K 


(2.9) 
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lt  can  be  seen  that  the  output  is  merely  a  weighted  sum  of  the  input 
values  separated  by  mT  seconds.  The  coefficients,  a,,  are  determined 
before  the  processing  starts  and  are  fixed  by  the  input  parameter  K. 
This  dependency  only  on  the  value  of  K  may  be  shown  in  the  following 
manner. 

Using  a  Hamming  Window  with  length,  N  =  LM,  where  L  =  2K  +  1, 
and  L  =  the  number  of  coefficients  used,  the  expression  for  the  window  is 

w(n)  =  0.54  -  0.46  cos  (27rn/N)  for  0<_n<_N  -  1 

=  0  elsewhere.         (2.10) 

Substituting  in  for  the  equivalent  expression  for  N: 

w(n)  =  0.54  -  0.46  cos  (2Tr(iM)/LM)  (2.11) 

where  n  =  iM  are  the  only  points  that  have  non-zero  values.  Now, 

w(iM)  =  0.54  -  0.46  cos  (27Ti/L)  (2.12) 

It  can  be  seen  that  the  coefficients  have  no  dependency  on  the  value 
of  M  or  fundamental  frequency.  This  shows  how  the  coefficients  may  be 
calculated  and  stored  before  processing  and  are  always  the  same  after 
the  parameter  K  has  been  chosen. 

There  are  several  factors  or  properties  of  windows  that  must  be 
considered.   First,  the  longer  the  unit  sample  response  of  the  window, 
the  narrower  the  passbands  of  the  comb  filter  will  be.  This  occurs 
when  the  fundamental  frequency  decreases.  Another  consideration  is 
the  fact  that  a  limit  exists  for  the  duration  of  the  window  due  to 
the  changing  pitch  period.   For  a  fixed  window  length,  the  choice  of 
window  will  change  the  filter  characteristics  such  as;  passband  width 
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14 
and  stopband  attenuation.    In  his  thesis,  Vaden  Shields  studied  the 

choices  for  type  of  window  and  window  length. 


Both  of  these  studies  were  experimental.  The  length  of  the  window 
was  directly  proportional  to  the  parameter  K  and  signals  were  processed 
with  and  without  an  additive  noise  signal  to  determine  the  effects  of 
distortion  and  separation  of  the  original  signal. 

The  experimental  observations  agreed  with  the  theoretical  results. 

As  the  value  of  K  was  increased  (this  corresponds  to  a  longer  window 

length,  more  coefficients,  and  narrower  passbands),  the  separation  of 

two  speakers  improved,  but  the  desired  signal  became  more  distorted. 

Shorter  window  lengths  yielded  less  separation  and  less  distortion. 

Through  a  series  of  listening  experiments  the  value  K  =  3  was  chosen 

as  a  good  compromise  value  from  the  distortion  and  separation  view- 

15 
points. 

The  search  for  the  optimum  window  type  started  with  the  immediate 
elimination  of  the  rectangular  window.  The  rectangular  window  intro- 
duced significantly  more  distortion  than  Hamming,  Hanning,  or  Blackman 
windows.  The  differences  between  these  three  windows  were  slight,  and 

after  listening  tests  the  Blackman  window  was  chosen  as  the  best 

16 
overal I . 

The  second  area  of  Shields'  thesis  dealt  with  the  treatment  of  un- 
voiced speech.  His  first  approach  was  to  turn  off  the  comb  filter  upon 
the  detection  of  a  pitch  period  greater  than  20  milliseconds.  This  ap- 
proach seemed  to  work  effectively  when  used  on  a  signal  without  any 
noise,  but  did  not  effectively  separate  the  desired  signal  from  the  noise, 
During  unvoiced  sections  of  speech,  the  noise  signal  would  "pop  out" 
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of  the  processed  signal.  This  effect  proved  to  be  distracting  to  listen- 

17 
ers  and  caused  an  impairment  to  intelligibility. 


Another  approach  for  the  treatment  of  the  unvoiced  speech  seg- 
ments was  used.  Whenever  an  unvoiced  segment  was  detected,  the  comb 
filter  would  continue  to  use  the  last  valid  pitch  period.  This  ap- 
proach achieved  better  results  for  noise  suppression,  but  it  tended 

18 
to  distort  the  desired  signal.    Shields  suggested  that  other  methods 

or  approaches  for  the  treatment  of  unvoiced  sections  should  be  investi- 
gated due  to  the  fact  that  no  one  method  had  proven  to  be  superior. 

Shields  stated  that  based  on  his  observations  of  the  results 
this  approach  to  speaker  separation  works.  The  performance  of  the 
system  was  measured  by  informal  listening  tests.  The  best  system 
performance  occurred  in  the  recovery  of  a  female  speech  segment  from 
a  male  -  female  combined  signal.  On  the  other  hand,  the  poorest  per- 
formance occurred  in  the  recovery  of  a  male's  speech  segment  from  a 
male  -  female  sum.  This  can  be  explained  in  the  following  manner. 
This  method  of  comb  filtering  allows  for  a  variable  number  of  pass- 
bands  or  "teeth"  which  depend  on  the  pitch  period  M.   If  M  is  the 
value  of  the  pitch  period  at  a  given  time,  then,  the  number  of  passbands 
up  to  the  Nyquist  frequency  (5  ki lohertz)  will  be  M/2.   In  a  female 
voice  the  pitch  period  is  generally  shorter  due  to  many  factors  such 
as;  size  and  shape  of  the  vocal  tract  and  articulators.  This  shorter 
pitch  period  would  cause  fewer  passbands  to  fall  within  the  frequency 
range  of  the  speech  spectrum.  With  fewer  passbands  present  there  is 
less  chance  of  falsely  recovering  an  unwanted  signal.  With  the  roles 

of  the  signal  reversed,  the  opposite  effect  would  occur  explaining  the 

19 
difficulty  in  extracting  the  male's  voice  from  the  combination. 
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During  the  listening  tests  several  factors  were  noticed  that 

caused  problems  in  this  system:   "unitel I igible  unvoiced  segments, 

short  segments  of  complete  distortion  of  the  signal,  and  breakthrough 

20 
of  the  undesired  phrase."    One  cause  of  the  distortion  problem  was 

attributed  to  the  rapid  changes  in  the  pitch  period  over  a  short  in- 
terval of  time.  Another  cause  of  the  distortion  stemmed  from  the 
inaccuracies  in  marking  the  pitch  period  which  was  carried  out  by  hand. 

The  problem  of  incomplete  or  insufficient  attenuation  of  the  un- 
wanted speech  segment  was  most  bothersome  during  periods  of  silence. 
Shields  suggested  that  if  periods  of  silence  could  be  distinguished 

from  unvoiced  segments,  then  unvoiced  segments  and  silent  segments 

21 
could  be  treated  separately. 

This  approach  for  comb  filtering  was  also  tried  with  white  noise. 

Shields  originally  thought  that  if  the  noise  were  wideband  with  most  of 

the  energy  outside  the  passbands,  the  signal-to-noise  ratio  would  be 

improved.   This  method  failed  to  produce  the  desired  result,  and  it 

changed  the  wideband  noise  into  noise  that  was  harmonic.  This  dis- 

22 
torted  the  output  and  produced  a  signal  that  had  a  "reedy"  sound. 

Using  noise  that  has  highly  harmonic,  this  method  showed  some 

promising  results.  The  suppression  and  distortion  were  much  better 

when  the  fundamental  frequency  of  the  noise  was  below  the  average 

pitch  of  the  speaker.  The  "pop  out"  effect  could  still  be  observed 

when  the  passbands  of  the  comb  filter  moved  into  the  range  of  the 

23 
frequencies  where  the  noise  was  present. 

The  remainder  of  Shields'  thesis  deals  with  the  automatic  extrac- 
tion of  pitch  period  from  a  combined  waveform  set  by  jpstral  analysis. 
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Since  another  type  of  pitch  detection  will  be  used  in  this  thesis,  this 
section  of  Shields'  thesis  will  not  be  reviewed. 

The  conclusions  that  Shields  reached  can  be  summarized  in  the 
following  manner.   First,  this  method  will  separate  the  combined 
speech  signals  of  two  speakers.  Second,  this  scheme  will  reduce  the 
effects  of  harmonic  noise.  Third,  possible  areas  of  improvement  in- 
clude the  treatment  of  unvoiced  segments  and  periods  of  silence.  The 
listening  tests  should  be  conducted  in  a  more  sophisticated  manner  in 
order  to  reveal  errors  in  the  system.  Fourth,  the  area  of  pitch 
detection  could  be  improved.  Also,  the  choice  of  coefficients  may  not 
be  optimal,  and  a  better  frequency  response  might  be  obtained  using 
another  set  of  coefficients.  Finally,  Shields  suggests  the  usage  of  a 

second  comb  filter  that  would  use  the  pitch  information  of  the  un- 

24 
wanted  speaker  to  block  or  reject  the  unwanted  speaker. 

2.2  Parsons'  System 

25 
In  a  recent  paper    written  by  Thomas  Parsons,  the  problem  of 

the  automatic  separation  of  the  simultaneous  speech  of  two  talkers  is 

approached  in  a  different  manner. 

Parsons  points  out  that  the  brain  requires  binaural  data  in 
order  to  separate  the  combined  signal,  and  in  some  cases,  a  signal  is 
received  over  a  channel  that  does  not  provide  this  binaural  informa- 
tion. Also,  Parsons  states  that  in  most  schemes,  the  signal  enhance- 
ment exploits  characteristics  of  one  signal  or  statistical  differences 

between  the  signals.  The  properties  or  statistics  of  speech  are  not 

26 
understood  we  I  I  enough. 
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The  method  that  Parsons  uses  to  separate  the  two  talkers  capital- 
izes on  the  harmonic  structure  of  short  segments  of  speech.   In  the 
frequency  domain  the  harmonic  structure  of  the  speech  appears  as  peaks, 
and  the  procedure  involves  removing  the  peaks  of  the  unwanted  speaker, 
and  then,  taking  the  inverse  transform  of  the  remaining  spectrum.  This 
is  a  rather  complicated  procedure  to  implement  automatically  for 
several  reasons.  The  frequency  peaks  from  the  two  speakers  will  over- 
lap in  some  areas.  The  overlap  problem  is  compounded  by  the  pitch 
variations  that  occur  in  natural  speech,  and  this  variation  causes  a 
frequency  modulation  (FM)  distortion  of  the  peaks.  Other  problems 

occur  from  the  event  of  the  pitch  tracks  crossing.  This  may  cause 

27 
errors  that  lead  to  jumping  between  speakers.       - 

Parsons  uses  four  rules  for  detecting  an  overlap  in  the  peaks, 

and  if  any  of  these  criteria  are  met,  then,  the  process  of  separating 

the  peaks  begins.  By  using  the  a  prior?  knowledge  of  the  peak  shapes, 

an  estimate  of  the  peak  may  be  subtracted  from  the  combined  peaks  in 

28 
order  to  resolve  them. 

The  next  item  addressed  is  the  determination  of  the  peak  shapes. 

The  largest  problem  is  the  effect  of  the  frequency  modulation  on  the 

peak  shape.  Parsons  states  that  by  assuming  the  pitch  variation  is 

linear  over  a  short  segment,  the  peaks  at  the  harmonics  can  be 

approximated  by  a  linear  FM  ramp.  He  states  that  for  slow  pitch  rates, 

the  peak  shape  is  approximately  the  same  as  in  the  case  of  a  constant 

pitch  period.  The  frequency  modulation  of  the  pitch  period  is  shown 

as  a  quadrature  component  with  an  amplitude  proportional  to  the  pitch 

rate,  and  the  peak  has  a  shape  that  is  the  second  derivative  of  the 

29 
in-phase  component.    This  approximation  for  the  peak  shape  holds  for 
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frequencies  up  to  3  kilohertz. 

The  separation  of  the  two  speakers  is  done  by  synthesizing  the 
spectrum' from  a  knowledge  of  the  pitch  contour  and  the  amplitude  and 
phase  of  every  one  of  the  harmonics.  The  method  of  synthesis  works 
better  than  the  method  of  subtracting  the  unwanted  speaker's  harmonics 
from  the  combined  spectra.  Parsons  states  that  if  the  subtraction 
methods  are  used,  any  errors  in  parameter  estimation  will  result  in 
incomplete  cancellation  of  the  undesired  speaker's  harmonics. 

Parsons  states  that  this  system  is  still  in  its  beginning  stages, 
and  that  problems  of  unvoiced  speech  still  have  to  be  solved.  Other 
problems  encountered  are  similar  to  those  mentioned  by  Shields.  These 
problems  include  pitch  detection,  tracking  one  speaker's  pitch  contour, 
and  the  areas  of  speech  that  have  pitch  rates  much  faster  than  in 
normal  speech. 

Both  methods  summarized  in  this  chapter  represent  possible  solu- 
tions to  the  speech  enhancement  problems,  and  there  are  probably 
other  existing  methods.  The  purpose  of  this  chapter  was  to  acquaint 
the  reader  with  some  of  the  methods  that  had  been  proposed,  and  not 
an  attempt  to  cover  all  possible  methods  for  speech  enhancement. 

2.3  Problem  Formulation  Introduction 


The  problem  at  this  point  was  to  decide  how  the  methods  for  speech 
enhancement  could  be  improved.  The  system  proposed  by  Shields  was 
used  as  a  guideline  because  his  work  had  been  the  only  source  of 
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i information  on  the  subject  before  Parsons  proposed  his  procedures. 

The  primary  concern  in  this  thesis  was  to  explore  the  methods 
used  for  speech  enhancement  in  voiced  sections.  The  unstructured  un- 
voiced speech  segments  were  studied,  but  because  of  its  lack  of  struc- 
ture, the  procedures  were,  for  the  most  part,  trial  and  error.  Chapter 
IV  formulates  the  method  that  was  proposed  as  an  alternative  to  the 
comb  filter.  The  method  is  developed  around  the  structure  and  time- 
varying  properties  of  the  speech  waveform. 
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CHAPTER  I  I  I 

PITCH  DETECTION  PROCEDURES 

The  need  for  a  simple  and  accurate  algorithm  for  determining 
pitch  period  or  fundamental  frequency  can  not  be  over  emphasized.   If 
the  separation  of  a  speaker  and  noise  is  to  be  done  in  some  small 
amount  of  time  by  the  methods  of  comb  or  adaptive  filtering,  most  of 
the  pitch  detection  algorithms  have  to  be  discarded  because  of  their 
complexity  and  execution  time.  The  accuracy  in  determining  pitch 
period  is  also  very  important.  The  filtering  techniques  that  are 
employed  in  this  thesis  use  the  pitch  period  information  to  set  the 
time-varying  digital  filter.  These  methods  of  filtering  assume  that 
the  pitch  period  information  is  exactly  known  for  a  given  speech  wave- 
form, and  therefore,  errors  in  pitch  detection  downgrade  system  per- 
formance. An  alternative  approach  for  pitch  detection  is  presented  in 
this  chapter. 

3,1   Research  Performed  by  Henke 

30 
In  a  recent  paper,  Henke  ,  discusses  the  features  of  an  accel- 

erometer  signal  that  measures  glottal  movements  may  lead  to  pitch 

period  detection  for  a  speech  waveform.  The  waveforms  shown  in 

figure  3-1  immediately  show  why  an  accelerometer  waveform  is  more 

desirable  from  the  standpoint  of  pitch  epoch  marking  than  the 

corresponding  speech  waveform. 
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SPEECH    WAVEFORM 


ACCELEROMETER  WAVEFORM 
SPEECH  AND  ACCELEROMETER  WAVEFORMS 


Figure  3-1 

Henke  states  that  a  small  accelerometer  can  be  attached  on  the 
midline  in  the  suprasternal  notch  (2  or  3  centimeters  below  the 
glottis)  to  measure  the  varying  pressure  waveform.  He  suggested  that 
the  most  prominent  and  stable  feature,  the  "flyback  stroke",  be  used 
In  pitch  period  detection  algorithms.  The  "flyback  stroke",  shown  in 
figure  3-2,  is  a  feature  of  the  pressure  waveform  that  is  character- 
ized by  a  rapid  change  from  outward  to  inward  acceleration.  This 
feature  occurs  immediately  following  the  maximum  outward  acceleration. 
The  maximum  outward  acceleration  occurs  at  the  instant  of  closure 
of  the  glottis  or  shortly  thereafter.  Therefore,  the  negative  going 
segment  of  the  waveform  that  crosses  the  zero  line  creates  a  stable 
point  for  pitch  epoch  determination.    After  these  pitch  epochs 
have  been  marked,  it  is  a  simple  task  to  determine  the  pitch  period 
or  fundamental  frequency  from  these  markings. 
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FLYBACK  STROKE 


FLYBACK   STROKE 


Figure  3-2 


3.2  Methods  of  Pitch  Detection  Used 


Henke's  method  was  used  in  the  work  carried  out  in  this  thesis. 
A  program  was  written  to  process  a  pressure  waveform  that  was  ob- 
tained from  an  accelerometer  at  the  same  time  as  the  speech  signal 
was  recorded.  This  was  done  in  a  two  channel  interleaving  mode  at  a 
sampling  rate  of  10  ki lohertz  with  the  speech  signal  on  one  channel 
and  the  accelerometer  signal  on  the  other. 

A  program  was  implemented  to  find  all  peaks  and  zero  crossings  in 
a  segment  of  the  accelerometer  waveform.  These  peaks  and  zero  cross- 
ing locations  were  stored  in  two  separate  arrays,  and  a  threshold  was 
used  to  determine  which  peaks  were  the  ones  that  were  associated  with 
the  instance  of  closure  of  the  vocal  cords.  The  first  threshold 
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method  consisted  of  simply  setting  a  constant  threshold,  by  guessing 
before  the  processing  took  place.  This  threshold  was  used  throughout 
the  program  to  determine  which  peaks  were  the  correct  ones.  The 
constant  threshold  approach  yielded  two  types  of  errors  that  can  be 
described  with  figure  3-3.   First,  the  problem  of  false  detection  or 
detecting  a  peak  that  was  not  associated  with  glottal  closure.  Second, 
the  false  rejection  problem  had  the  complementary  effect,  in  that,  a 
true  peak  was  not  large  enough  to  be  detected  by  the  constant  threshold 
The  constant  threshold  method  yielded  about  40  to  50  percent  total 
errors  of  both  types,  and  the  constant  threshold  was  found  to  be 
speaker  and  gain  dependent.  The  errors  incurred  by  this  method 
suggested  a  second  approach  for  the  threshold  determination. 


-THRESHOLD 


FALSE    DETECTION   AND  FALSE   REJECTION 


Figure  3-3 
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The  problems  in  the  first  approach  resulted  from  the  dynamic  range 
of  the  accelerometer  waveform.  The  second  approach  consisted  of  set- 
ting a  variable  threshold.   The  variable  threshold  was  set  by  means  of 
an  energy  measurement  of  the  glottal  signal  over  a  short  time  interval. 
The  energy,  E,  was  measured  in  the  following  manner; 


N-1 


E-  I 


n=0  <- 


x(n) 


(3.1) 


where  N  was  taken  to  be  25.6  milliseconds. 

The  range  of  the  values  for  the  glottal  peaks  was  divided  into 
five  levels,  and  after  a  careful  examination  of  a  few  waveforms  and 
energy  measurements,  these  ranges  were  experimentally  determined.  The 
results  of  the  variable  threshold  approach  were  quite  encouraging.  The 
number  of  errors  made  was  on  the  order  of  20  percent,  and  this  figure 
could  be  improved  if  a  more  careful  statistical  study  were  made  to 
determine  the  optimal  threshold  settings  for  a  given  energy  measure- 
ment value. 

The  idea  of  the  energy  measurement  also  suggested  a  concept  for 
silence  detection.   In  the  thesis  done  by  Shields,  he  suggests  that  if 
silence  could  be  detected,  then  some  other  approach  could  be  taken  to 
improve  separation.  The  silence  detector  was  included  in  the  program 
in  order  that  different  schemes  of  filtering  could  be  tried.  A 
silent  area  was  determined  when  the  value  of  E  was  equal  to  zero. 
[The  silent  areas  were  marked  in  the  table  of  pitch  period  values  by 
a  minus  sign  in  front  of  the  pitch  period  value.J 
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After  the  peaks  were  detected,  zero  crossings  in  the  accelero- 
meter  signal  immediately  following  the  peaks  were  marked.  The  zero 
crossing  was  the  termination  of  the  "flyback  stroke",  and  as  mentioned 
before,  Is  thought  to  be  one  of  the  more  prominent  and  stable  features 
in  the  glottal  waveform.  These  zero  crossings  were  marked  by  using 
the  second  channel  of  the  file,  since  the  signal  only  required  one 
channel . 

Since  the  minimum  pitch  period  is  around  2.0  milliseconds,  the 
algorithm  included  a  rule  that  repositioned  the  pointer  in  the  wave- 
form by  20  samples  after  a  correct  peak  was  detected.  This  helped  to 
speed  up  the  algorithm  somewhat. 

There  was  another  problem  that  is  still  associated  with  this 
method  of  pitch  determination.  The  transitions  just  before  voiced  and 
unvoiced  sections  were  not  of  the  same  form  as  during  voiced  sections 
or  unvoiced  sections.  Even  looking  at  the  section  of  the  waveform 
that  was  considered  a  transition  area  did  not  help  in  the  decision  of 
where  to  put  the  marks  for  the  pitch  period.  Each  signal  had  to  be 
marked  by  hand  in  these  areas,  and  fortunately,  there  were  not  many  of 
them  in  a  waveform. 

The  outputs  of  the  pitch  detection  included  two  separate  files. 
One  file  consisted  of  the  pitch  marks.  This  file  could  be  added  to  a 
speech  signal  for  the  filtering  of  the  signal,  or  it  could  be  added  to 
an  acce lerometer  signal  to  determine  whether  the  pitch  detection 
algorithm  had  made  any  errors.  The  second  output  file  consisted  of 
pitch  period  values  that  were  determined  from  measuring  the  distances 
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between  the  marks  In  the  previous  file.  The  file  is  referred  to  as 
the  pitch  table,  and  it  is  used  in  the  filtering  programs. 

The  methods  used  in  this  section  were  intended  to  make  the  pitch 
detection  as  easy  as  possible.  The  filtering  schemes  used  in  the 
thesis  assumed  that  the  pitch  information  was  known  exactly.  The  pri- 
mary function  of  the  automatic  pitch  marking  program  was  to  eliminate 
as  many  periods  as  possible  that  would  have  to  be  marked  by  hand.   If 
pitch  period  detection  accuracy  had  been  required  in  the  automatic 
marking  program,  more  improvements  would  have  been  made. 


-38- 
CHAPTER  IV 

ADAPTIVE  FILTERING  CONCEPTS 
4. 1   Development 

The  idea  for  the  pitch  synchronous  adaptive  filter  was  developed 
from  the  ideas  of  Siamak  Samsam  as  an  alternative  to  the  conventional 
comb  filtering  techniques. 

The  name  comb  filter  probably  came  about  from  its  frequency 
response  characteristics.  The  frequency  response  of  a  comb  filter, 
shown  in  figure  4-1  consists  of  a  series  of  "teeth"  or  "fingers". 
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The  implementation  used  for  the  comb  filter  in  figure  4-1  was  discussed 
in  section  2.1.  The  basic  idea  behind  this  implementation  was  to  pad 
the  impulse  response  of  a  low-pass  window  function  with  zeros,  using 
the  same  number  of  zeros  between  the  coefficients  as  the  pitch  period 
of  the  signal  at  that  particular  time.  The  spacing  between  the  co- 
efficients of  the  impulse  response  was  uniform  giving  the  frequency 
response  in  figure  4-1. 

In  this  type  of  comb  filtering  implementation  the  length  of  the 
filter  was  on  the  order  of  five  to  nine  times  the  length  of  one  pitch 
period,  and  this  is  equivalent  to  a  range  of  filter  lengths  between 
40  and  140  milliseconds.   In  a  normal  speaking  voice  the  pitch  period 
or  fundamental  frequency  can  change  rapidly  over  several  periods 
depending  on  intonation  or  stress,  and  if  the  filter  length  is  long, 
the  comb  filter  may  not  handle  the  correct  samples  by  the  procedures 
discussed  in  section  2.1.  The  phenomenon  discussed  above  can  pro- 
bably be  more  easily  explained  by  figure  4-2. 


Figure  4-2 
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Suppose  that  the  signal  in  figure  4-2  is  a  speech  waveform  with 
pitch  periods  T  ;  where  p  =  0  to  5.   It  can  be  seen  that  the  values  of 
the  pitch  periods  vary  over  the  length  of  the  filter.   Suppose  that  the 
number  of  zeros  between  coefficients  has  just  been  changed  to  cor- 
respond to  the  pitch  period  T  ,  and  all  spacings  between  the  samples 
are  adjusted  so  that  they  correspond  to  the  period  T0.  This  results 
in  a  misplacement  of  the  coefficients  in  the  waveform  with  respect  to 
the  beginning  of  each  pitch  period.   In  the  case  of  a-?,  its  position 
in  that  particular  period  is  much  different  from  the  position  of  co- 
efficients aQ,  a^,  and  a2» 

Now,  suppose  that  the  pitch  period  T  ,  has  a  smaller  error  due  to 
the  pitch  detection  program.  Since  all  spacings  are  uniform,  the 
error  is  propagated  and  compounded  in  the  subsequent  spacings. 

With  some  insight  into  what  the  comb  filtering  was  trying  to  do, 
the  adaptive  methods  were  formulated.  Suppose  that  the  coefficients 
are  rearranged  so  that  the  spacings  are  no  longer  uniform,  but  in- 
stead, the  spacings  will  be  made  equal  to  the  pitch  period  for  that 
particular  point  in  the  waveform.  Figure  4-3  shows  the  placement  of 
the  coefficients  for  the  adaptive  filtering  methods.  For  this  method 
each  coefficient  is  in  the  same  relative  position  in  its  particular 
period. 

This  method  does  not  really  correspond  to  a  comb  filter  as  shown 
by  the  frequency  response  of  a  sample  adaptive  filter  in  figure  4-4. 
The  most  substantial  thing  that  can  be  said  about  the  structure  of  the 
frequency  response  shown  in  figure  4-4  is  that  it  has  no  general 


-41- 


Figure  4-3 

structure.  The  structure  of  figure  4-4  points  out  the  frequency  res- 
ponse of  the  adaptive  filter  provides  little  insight  into  its 
capabilities  in  the  speech  enhancement  problem. 


Figure  4-4 
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The  output  of  this  system  can  be  expressed  in  the  following 

equation : 

2K 
y(n)  =  I       aj  x(n-|  j ) 
i=0 

where:    aj  =  the  coefficients  or  weightings 

I,  =0 

[2     =  m1 

I  3  =  m  i   +   mo 


!  I   =  m.  +  rru  +  .  .  .  mj  _j 
and,      mj  =  the  pitch  period  of  the  Pn  period  from  the 
originating  point. 
2K  +  1  =  the  number  of  coefficients  used. 


There  is  another  method  of  expressing  the  ideas  given  earlier. 
This  description  is  aided  with  the  use  of  figure  4-5. 

The  method  can  be  thought  of  in  the  following  manner.  First, 
the  waveform  is  broken  into  segments  according  to  the  pitch  epoch  or 
positions  where  the  pitch  marks  occur.  Then,  these  segments  are  lined 
up  as  shown  in  figure  4-5.  Then,  a  weighted  average  is  computed  point 
by  point  as  the  filter  moves  in  the  direction  indicated.  What  this 
operation  intuitively  does  is  to  develop  an  average  impulse  response 
based  on  the  previous  several  periods.  This  technique  works  well  in 
conjunction  with  the  assumption  that  the  impulse  response  of  the  vocal 
tract  is  slowly  varying.  Another  important  point  to  be  said  for  this 
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Figure  4-5 

technique  is  the  fact  that  in  the  separation  problem,  the  weighted 
samples  of  the  desired  speaker  will  add  coherently  while  the  contri- 
butions of  the  undesired  speaker  will  not  add  coherently.   In  the 
comb  filtering  methods  the  same  statement  can  not  be  made  as  strongly 
when  the  pitch  is  changing.   Since  the  averaging  process  in  the  comb 
filter  looks  at  different  relative  positions  in  each  period,  the 
reinforcement  or  coherent  addition  of  the  desired  speaker  is  not  as 
pronounced  as  in  the  adaptive  methods. 


It  should  be  pointed  out  that  both  methods  are  identical,  and  in 
fact,  a  comb  filter  when  the  pitch  period  is  constant.   In  other 
words,  the  adaptive  filter  approaches  a  comb  filter  when  the 
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pitch  period  contour  is  slowly  varying  and  performs  a  comb  filter  opera- 
tion when  the  pitch  period  contour  is  constant.  When  the  pitch  con- 
tour is  fluctuating  rapidly,  the  adaptive  filter  does  not  resemble  a 
comb  filter  in  the  least. 

There  is  one  problem  in  this  method,  and  this  problem  with  a 
possible  solution  is  discussed  in  the  next  section. 

4.2  Overload  Problem 

A  problem  in  the  adaptive  filtering  method  referred  to  as  the 
"overload  problem  or  phenomenon"  can  be  described  with  the  aid  of 
figure  4-6. 


Figure  4-6 
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Suppose  that  there  is  a  rapid  increase  in  the  pitch  period  in  a 
relatively  short  interval  of  time.  This  phenomenon  is  not  rare  in 
normal  speech.   For  some  speakers  a  distinguishable  feature  or 
characteristic  is  the  short  pitch  period  that  occurs  at  the  onset  of 
voiced  speech.   This  short  period  can  be  represented  by  the  segment 
in  figure  4-6  labeled  as  T2.  As  the  filtering  moves  to  the  right  com- 
puting a  point  by  point  average,  the  method  runs  into  problems  at  the 
point  T~.   Up  to  the  point  T2  the  filtering  has  occurred  in  the  cor- 
rect manner.   After  T„  the  coefficient  a~  begins  to  use  samples  in  seg- 
ment two  for  the  averaging.  This  feature  disagrees  with  the  concepts 
of  the  adaptive  filtering  methods  because  the  coefficients  now  are  not 
in  the  same  relative  positions  in  their  respective  segments.  This 
phenomenon  has  been  named  the  overload  problem  and  causes  an  undesirable 
result  in  the  output  waveform.  The  overload  phenomemon  will  be  studied 
more  in  the  next  chapter. 

There  are  several  proposals  for  solutions  to  the  problem  of  over- 
loading.  First,  the  easiest  solution  from  all  viewpoints  is  to  ignore 
the  problem  and  let  the  overloading  occur.   If  the  pitch  period  is 
s lowly  varyi  ng,  the  prob I  em  will  be  minor  and  the  coefficients  will 
not  be  far  away  from  the  same  relative  position  in  their  respective 
periods.  Second,  the  adaptive  filter  could  be  "turned  off"  at  T  and 
pass  the  input  signal  completely  for  the  time  interval  between  T  and 
T0  .  This  idea  originated  from  the  fact  that  the  impulse  response  of 
the  vocal  tract  would  be  of  sufficiently  low  amplitude  in  the  tails 
of  the  impulse  response  so  that  the  filtering  operation  would  not  have 
much  effect  in  this  area.  The  final  solution  proposal,  and  the  one 
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that  was  used  in  this  thesis  was  one  that  padded  the  short  segments  in 
the  filter's  length  with  zeros  in  order  to  make  the  length  of  the 
short  segments  equal  to  the  one  in  the  front  of  the  filter. 

Figure  4-7  explains  this  procedure.  The  shortest  period,  seg- 
ment 3,  will  have  zeros  added  in  the  area  indicated. 


ZEROS  ADDED 


CORRECTING  OVERLOAD    PROBLEM 


Figure  4-7 


When  the  filter  reaches  the  end  of  segment  4,  the  same  problem  occurs 
and  zeros  are  again  added  to  pad  out  this  segment.  In  effect  this  is 
the  same  procedure  as  "turning  off"  the  individual  coefficients  asso- 
ciated with  the  segments  that  have  been  padded  with  zeros.   In  figure 
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4-7  this  would  involve  setting  a2  equal  to  zero  after  T~  is  reached 
and  a-^  equal  to  zero  after  T^  is  reached. 

There  is  one  other  factor  to  be  considered  in  this  method.  The 
gain  associated  with  the  filter  has  been  fixed  to  be  a  unity  gain 
system  by: 


2K 
I         ak  =  1  (4.2) 

k  =  0 


When  a  particular  coefficient  is  set  to  zero  or  "turned  off",  the 
remaining  coefficients  must  be  rescaled  so  that  the  output  is  not 
attenuated. 

The  third  procedure  for  handling  the  overload  problem  makes 
sense  from  an  intuitive  standpoint.   It  says  that  if  a  coefficient 
can  not  be  placed  in  the  same  relative  position  in  the  segment  or 
period  due  to  the  shortness  of  that  period,  then,  the  contribution  of 
that  coefficient  will  only  provide  erroneous  results,  and  therefore, 
it  should  not  be  considered. 

4.3  Rules 

As  the  filtering  algorithms  become  more  complicated  and  complex, 
the  list  of  rules  for  handling  various  situations  also  grows.  Some 
of  the  more  important  rules  used  will  be  discussed  in  this  section. 

The  discussion  thus  far  has  dealt  only  with  the  procedures  for 
handling  voiced  speech.  Normal  speech  waveforms  contain  unvoiced 
sections  also,  and  these  sections  must  also  be  handled.  There  are 
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three  basic  areas  or  situations  involved  in  the  processing  of  unvoiced 
speech  sections.  These  situations  are: 

1.  The  detection  of  an  unvoiced  section,  and  transition 
from  voiced  speech  procedures  to  unvoiced  speech 
procedures. 

2.  The  procedures  that  are  used  when  the  filter  is  entirely 
in  the  unvoiced  section. 

3.  The  detection  of  a  voiced  section,  and  transition  from 
unvoiced  speech  procedures  to  voiced  speech  procedures. 

The  first  situation  mentioned  above  was  handled  by  using  the  last 
voiced  pitch  period  as  shown  in  figure  4-8. 

In  figure  4-8  the  speech  waveform  has  been  omitted  for  clarity, 
but  the  sections  are  marked  as  voiced  and  unvoiced.  The  values  labeled 
m|  are  equal  to  the  pitch  period  of  that  particular  period.   In  figure 
4-8  (a)  the  filter  position  is  shown  just  prior  to  an  unvoiced  area 
detection.   Figure  4-8  (b)  shows  the  spacing  after  the  first  coeffi- 
cient has  entered  the  unvoiced  area.   In  this  figure  two  of  the 
spacings  have  a  spacing  of  m<  while  the  other  two  spacings  are  related 
to  their  particular  periods.   In  figure  4-8  (c)  another  spacing 
change  is  shown  so  that  three  spacings  are  equal  to  m^ .  This  pro- 
cedure is  repeated  until  the  entire  filter  lies  in  the  unvoiced  area. 

Now  the  second  situation  has  been  entered,  and  another  scheme  is 
begun.  Since  the  method  of  attenuating  the  input  in  an  unvoiced 
section  worked  satisfactorily  in  the  system  proposed  by  Shields,  it 
was  used  in  this  system  also.   In  effect,  the  filtering  was  terminated, 
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Figure  4-8 

and  the  input  waveform  was  multiplied  by  a  constant  less  than  one  in 
order  to  produce  the  output  waveform. 

In  the  third  situation  the  procedure  was  basically  the  inverse  of 
the  method  used  in  the  first  situation.  When  a  voiced  section  was 
detected,  the  filter  was  initialized  to  a  constant  spacing  using  the 
first  period  value.  As  the  filter  moved  further  into  the  voiced 
sections,  the  spacings  were  changed  to  conform  to  the  pitch  periods. 


In  summarizing  this  chapter  the  most  prominent  feature  of  the 
adaptive  filtering  approach  is  its  conformity  with  the  speech  waveform. 
The  performance  of  this  system  will  be  evaluated  in  later  chapters. 
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CHAPTER  V 

TEST  SIGNAL  SECTION 

5.1  Development  of  Test  Signals 

It  was  decided  that  before  any  new  methods  for  speech  enhancement 
were  developed,  a  method  of  testing  the  various  schemes  should  first 
be  undertaken.  Even  though  psychological  listening  tests  were  to  be 

conducted  in  the  future,  there  were  just  too  many  variables  in  the  system 
to  be  handled  in  a  complex  listening  test.   In  order  to  assign  some 
types  of  performance  indicies,  and  in  order  to  fix  some  of  the  vari- 
ables, a  test  input  signal  was  generated  and  stored  for  future  use. 

The  purpose  of  the  test  input  signal  was  to  alleviate  the  un- 
certainty of  characteristics  of  the  speech  waveform,  such  as,  the 
pitch  period  and  the  impulse  response.  When  the  pitch  period  and  the 
impulse  response  of  a  voiced  section  of  speech  are  known,  then,  some 
types  of  input  and  output  waveform  comparisons  can  be  made  in  order  to 
determine  optimal  filter  properties. 

A  test  input  signal  was  formulated  in  the  following  manner: 
First,  it  was  decided  that  enough  comparison  data  could  be  obtained 
from  a  relatively  short  segment  of  the  input  signal.  This  feature 
would  provide  results  that  were  accurate,  and  in  addition,  they  could 
be  obtained  quickly.  Second,  it  was  decided  that  the  impulse  response 
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of  the  generation  filter  would  be  time-invariant  over  the  duration  of 
the  signal.  This  assumption  was  chosen  in  order  to  make  the  compari- 
son study  easier,  and  it  conforms  to  the  properties  of  speech,  in  that, 
the  impulse  response  of  the  vocal  cords  over  short  time  segments  can 
be  considered  constant.  Third,  it  was  decided  that  the  pitch  period 
would  vary  sinusoidal ly  about  some  mean  pitch  period. 

The  input  signal  formulation  can  be  described  by  the  following 
diagrams  and  equations: 
Let  xCn)  =  the  test  input  signal. 

v(n)  =  the  impulse  response  of  the  generating  filter. 

w(n)  =  the  nonuniformly  spaced  impulse  train  corresponding  to 
the  pitch  pulses  or  the  excitation. 
Then,  x(n)  =  v(n)  *  wCn)  or  (5.1) 


Figure  5-1 


The  impulse  response  of  the  generating  filter,  v(n),  was  chosen 
to  be  a  damped  sinusoid  of  the  following  form: 
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Figure  5-2 

where:  v(n)  =  e*"n/a   sin  2irwn     0  <   n  <_  L-1  (5.2) 

The  model  seems  to  be  a  very  good  model  for  the  impulse  response  of 
the  vocal  cords.  The  parameters  a  and  w  were  variable  and  were 
chosen  so  that  the  input  signal  that  was  generated  would  be  a  good 
model  of  speech. 

The  length  of  the  impulse  response  was  chosen  to  be  2T  ,  where 
TQ  is  the  value  of  the  mean  pitch  period.   It  was  felt  that  in  order 
to  model  speech  as  closely  as  possible,  the  length  of  the  impulse 
response  of  v(n)  should  be  longer  than  the  pitch  periods  encountered. 
This  would  cause  some  overlap  of  the  impulse  response,  v(n),  in  the 
input  waveform. 


The  nonuniformly  spaced  impulse  train,  w(n),  was  formulated  in 
the  following  manner:  The  spacing  between  the  impulses,  (the  pitch 
period)  was  chosen  to  vary  sinusoidal ly  about  some  mean  pitch  period 
TQ  .   The  figure  5-3  shows  the  procedure  for  this  development.   It 
can  be  seen  that  this  nonuniformly  spaced  impulse  train  is  periodic 
on  the  interval  K),  N  j  . 
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where:  6-    =   6  sin  2  tt 


Figure  5-3 

for  i  =  0  to  N-1 


(5.3) 


Now  with  these  two  components  the  input  signal  x(n)  was  gen- 
erated. This  gave  a  waveform  whose  impulse  response,  and  pitch 
contour  were  exactly  known,  and  the  various  methods  of  speech 
enhancement  could  be  tried  on  this  signal. 

5.2  Processing  of  Test  Signals 


With  the  test  signals  that  have  been  described  in  the  previous 
section,  several  processing  methods  may  now  be  more  closely  examined. 
The  filtering  schemes  used  in  this  section  were  designed  so  that  they 
resembled  the  actual  filtering  systems  in  all  respects  except  for  the 
length  of  the  input  waveforms  used.  Since  the  test  input  waveform  was 
much  shorter  than  the  actual  speech  waveforms,  it  could  be  completely 
stored  in  core  memory,  and  the  amount  of  time  needed  for  the  signal 
processing  was  short. 

There  were  three  filtering  systems  implemented  for  the  test  signal 
input  waveform.  These  systems  included  the  system  proposed  by  Shields, 
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the  adaptive  filtering  systems  proposed  in  Chapter  IV  with  the  over- 
load problem,  and  the  adaptive  filtering  system  without  the  overload 
problem.   For  the  remainder  of  the  thesis  the  adaptive  system  with  the 
overload  problem  will  be  referred  to  as  the  adaptive  overload  system. 
The  adaptive  system  that  compensates  for  this  problem  will  be  referred 
to  simply  as  the  adaptive  system. 

Another  problem  that  received  much  thought  was  the  problem  of  com- 
paring these  system  in  some  manner  to  determine  the  good  and  bad  points 
of  both  systems.  Several  tests  were  decided  upon  that  were  related  to 
the  problems  of  speech  enhancement.   In  the  thesis  done  by  Shields,  the 
tradeoffs  between  desired  speaker  distortion  and  undesired  speaker 
separation  were  described.  These  tradeoffs  were  also  discussed  in  Chap- 
ter II.   It  was  decided  that  two  separate  tests  could  be  conceived  for 
this  area. 

The  first  test  will  be  referred  to  as  the  input  pass  test.   Figure 
5-4  describes  the  manner  in  which  this  was  implemented.   By  using  the 
enhancement  system  in  this  manner,  the  output  waveforms  should  exactly 
resemble  the  input  waveform,  or  in  other  words,  the  system  should  act 
as  an  identity  system.  By  examining  the  inputs  and  the  outputs  of  the 
systems,  the  amount  of  desired  speaker  distortion  could  be  determined 
for  various  numbers  of  coefficients  or  windows  used. 

The  second  test,  described  by  figure  5-5,  will  be  referred  to  as 
the  input  reject  test.   In  this  type  of  filtering  the  pitch  information 
is  provided  from  another  input  waveform  that  is  totally  uncorreiated 
with  the  test  input  signal.  The  idea  behind  this  type  of  processing 
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Figure  5-4 

is  to  examine  the  amount  of  attenuation  introduced  into  the  output 
waveform.   Ideally,  the  output  would  be  zero,  because  the  enhancement 
system  would  be  passing  another  waveform  with  a  different  pitch  con- 
tour.  This  test  was  designed  in  order  to  measure  the  amount  of 
speaker  separation  that  resulted  from  the  various  systems. 
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Figure  5-5 
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Another  test  was  devised  in  order  to  examine  the  capability  of 
each  system  when  white  noise  was  added  to  the  input  signal.  Figure 
5-6  shows  the  implementation  of  this  system. 
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Figure  5-6 

Shields  stated  that  the  system  proposed  in  his  work  seemed  to  turn 

white  noise  into  a  highly  harmonic  noise  in  the  output  waveform  that 

32 
was  very  annoying.     This  test  was  designed  to  determine  if  the 

adaptive  filter  had  a  problem  with  white  noise,  and  if  it  did,  the 

systems  could  be  compared  to  determine  which  one  did  the  best  in  these 

ci  rcumstances. 


Finally,  a  test  could  be  performed  on  the  sum  of  two  waveforms 
that  would  duplicate  the  problems  of  the  overall  speech  enhancement 
system.   It  should  be  pointed  out  that  this  test  was  merely  a 
combination  of  tests  one  and  two,  and  since  the  systems  were  linear, 
superposition  would  have  given  the  same  results. 
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In  all  cases  it  was  assumed  that  the  entire  test  signal  was 
periodic  with  a  period  of  twenty  pitch  periods,  and  the  processing 
only  concerned  itself  with  one  period  of  the  overall  waveform. 

It  should  be  pointed  out  that  the  test  signals  used  in  these  tests 
were  not  actual  speech  waveforms.  The  input  waveform  was  modeled  to 
closely  resemble  the  speech  waveform  during  voiced  sections.  There 
were  no  comparison  tests  made  on  test  signals  modeled  as  unvoiced 
speech.   In  the  next  section  the  results  of  these  tests  for  the  differ- 
ent systems  implemented  will  be  discussed. 

5.3  Results  of  Test  Signal  Processing 

In  order  that  the  results  of  the  test  signal  processing  be  pre- 
sented as  clearly  as  possible,  many  figures  will  be  employed  in  this 
section.  There  are  three  basic  types  of  presentation  used  in  this 
section.   First,  a  time  domain  presentation  is  used  to  show  how  the 
various  systems  process  the  input  waveforms,  and  from  these  figures, 
the  input  and  output  waveforms  can  be  compared  visually.  Second,  a 
frequency  domain  presentation  is  used.  The  test  signal  was  formulated 
to  simulate  a  signal  that  had  been  sampled  at  10  kilohertz.  A  short 
time  spectrum  was  computed  by  resampling  the  input  waveform  at  5  kilo- 
hertz,  multiplying  by  a  Hamming  Window,  and  then,  the  spectrum  was 
computed  using  a  Fast  Fourier  Transform  (FFT)  program.  The  bandwidth 
of  the  test  signal  was  well  below  the  Nyquist  Frequency  for  the  new 
sampling  rate.  The  frequency  domain  presentation  displays  logarith- 
mically the  square  of  the  magnitude  of  the  Fourier  Transform  of  the 
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slgnals,  and  this  was  computed  to  eliminate  the  effects  of  any  phase 
errors  that  were  introduced  by  the  filtering.  Finally,  an  error  func- 
tion is  shown  in  some  cases,  and  this  error  function  can  be  expressed 
in  the  following  equation: 

E  =  20  log{|X(eJ'w)|2  }     -  20  log{  |  Y(eJ'w)  |2  }     (5.4) 
where:   X(eJw)  is  the  Fourier  Transform  of  the  input  signal 

Y(eJw)  is  the  Fourier  Transform  of  the  output  signal. 

Notice  that  the  error  function  defined  above  can  be  thought  of  in 
another  manner.  Consider  the  linear  system  shown  in  figure  5-7  below: 


Figure  5-7 


Clearly, 


y(n)  =  x(n)  '*  h(n) 


(5.5) 


Expressing  equation  (5.5)  in  the  frequency  domain 


Y(e>)   =  X(e>)  H(eJW) 


(5.6) 


Now,      H(ejw)  =  Y(e^) 

X(eJw) 


(5.7) 
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Taking  the  magnitude  and  squaring  both  sides  gives 


H(eJw) 


2  _ 


Y(eJw) 


X(eJw) 


2  =  |v(eJ>)|2 
|X(eJ'w)|2 


(5.8) 


Now,  taking  the  logarithm  of  both  sides  gives 


log |H(eJw) |2  =   I 


og 


[XTeWl2 


logflYCeJ^)!2}  -  |og{  |X(eJw)|2}     (5.9) 

Notice  that  this  is  almost  the  same  form  as  the  error  function.  l-KeJ03) 

is  often  referred  to  as  the  system  transfer  function,  and  the  error 
function  may  be  represented  as: 


E  =   log 


1 


HCeJw)l2 


(5.10) 


Therefore,  the  error  function  could  be  thought  of  as  being  related  to 
the  inverse  system  transfer  function. 

Before  proceeding  the  following  discussion  will  be  covered  in  or- 
der to  show  how  the  error  function  will  be  used  to  calculate  system 
performance.  Consider  the  system  proposed  in  figure  5-7.  Suppose  that 
the  following  result  was  desired: 


y(n)  =  a  x(n) 


(5.11) 


Expressing  this  in  the  frequency  domain  and  because  the  system  is 
I inear: 


Y(eJw) 


a  X(eJw) 


(5.12) 


Now, 


XCeJw) 
YCeJ") 

1_ 
a 

log 

|x(eJ'w)|2 
[Y(eM|2 

= 

log 

_R2_ 

E  =   log{|X(eJw)[2}  -  log{ |Y(eJw) |2}  =   log 


E  =   log  (1)  -  2  log{|a|}  =  -  2  log{|a|} 


1 

a|2 
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(5.13) 

(5.14) 

(5.15) 

(5.16) 


The  above  results  will  be  used  in  later  parts  of  this  section  for 
comparative  analysis  between  the  systems  implemented. 

The  first  test  described  in  section  5.2  was  named  the  input  pass 
test.  This  test  was  performed  with  the  three  basic  systems  developed: 
Shields',  the  adaptive  overload,  and  the  adaptive.   For  all  systems 
five  coefficients  were  used,  and  these  coefficients  were  taken  from 
the  Hamming  Window  function. 

Figure  5-8  shows  the  time  domain  input  and  output  waveforms  for 
Shields'  system.  Several  things  can  be  noticed  in  the  output  waveform, 
First,  there  is  a  definite  amplitude  modulation  throughout  the  wave- 
form. Second,  there  are  some  areas  where  an  "overload  bump"  occurs 
in  the  output.   Figure  5-9  shows  the  same  presentation  for  the  adap- 
tive overload  system.   It  can  be  noticed  that  the  amplitude  modulation 
does  not  exist,  but  there  are  still  some  problems  with  overloading. 
In  figure  5-10,  the  adaptive  system  results  are  shown:   input  and  out- 
put waveforms  are  essentially  identical. 
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A  similar  comparison  can  be  made  by  viewing  the  spectra  of  the 
input  and  output  waveforms.   Figure  5-11  shows  the  spectrum  of  the 
input  waveform,  while  figures  5-12,  5-13,  and  5-14  show  the  output 
spectra  of  Shields'  system,  the  adaptive  overload  system,  and  the 
adaptive  system  respectively.   Figure  5-15  presents  the  error  function 
for  each  of  the  systems. 

For  this  test  the  system  should  be  performing  as  an  identity  sys- 
tem. This  says  that  the  value  for  a_  in  equation  (5.11)  should  be 
equal  to  unity.   Substituting  this  into  equation  (5.16): 

E   =   -2  log  (1)   =  0  (5.17) 

This  result  is  approached  in  all  three  systems,  but  the  adaptive  system 
comes  the  closest  to  matching  the  ideal  error  function  for  this  test. 

From  all  three  presentations,  the  adaptive  system  clearly  performs 
much  better  than  the  other  two  systems  on  the  input  pass  test. 

The  second  test  performed  was  named  the  input  reject  test.  The 
purpose  of  this  test  was  to  determine  how  well  the  systems  would  re- 
ject or  attenuate  an  undesired  speaker.  This  corresponds  to  the  separ- 
ational  aspects  of  speech  enhancement.  The  manner  in  which  this  test 
was  conducted  was  simple.   A  waveform  was  generated  with  different 
pitch  periods,  and  the  pitch  marks  from  that  signal  along  with  its 
pitch  table  were  used  to  filter  the  original  waveform. 

The  second  waveform  had  a  pitch  period  that  also  varied  sinusoid- 
ally  about  some  average  value  that  was  identical  to  the  original 


-68- 


■g<3  ni  Q3avn6s  BaniiNovw 


3 
CD 


-69- 


N 
h- 

cr 

Ll) 

X 

o 


>- 
o 

UJ 

o 
cr 


to 


ZJ 
C 


>- 

CO 


"O 


CO 


3 
l_ 

+- 
U 
<D 
Cl 

co 


ZJ 
Q- 


Z> 

O 


CN 


O 
CN 


LA 


in 


in 


■go  ni  3yvn6s  aaniiNovw 


m 

CO 

i_ 


■70- 


o 

C\l 


LT\ 


in 


o 


I 


'aa  ni  QByvnos  3aniiN9vw 


N 


cr 

UJ 

X 

_j 

in 

o 

CM 

^ 

^ 

>- 

o 

?r 

UJ 

D 

O 

LU 

cc 

u_ 

o 

• 

CM 

W 

in 
<o 
Q_ 


CL 

C 


E 
0) 

+- 
in 
>- 

CO 

"O 

ro 
O 


0 
> 

O 

CD 
> 

+- 

O- 

ro 

TO 

< 

I 

E 
a 

L_ 
+- 

o 

CD 
CD- 
CO 


3 
CL 


I 

CD 

1_ 
13 
(J) 


-71- 


M 
UJ 

o 


^ 


> 
o 


o 

LU 


CO 
0) 


CO 


z> 

CL 

C 


E 
0 

4- 
tn 

CO 

CD 
> 


CL 

ro 

"a 


i 


ZJ 
L. 
4- 
U 

a) 

CL 

CO 


Z> 
Q. 


3 

o 


in 


in 


in 


*aa  ni  a3yvn6s  3amiN9Vw 


i 

CD 
1_ 
3 
CD 


-72- 


-73- 


-74- 


o 


LTi 


in 


'aa  ni  aadvnOs  aaniiNsvw 


-75- 
waveform.   The  6  of  the  second  waveform  was  changed  in  order  to  give 
the  different  set  of  pitch  periods.   It  should  be  pointed  out  that 
this  procedure  allows  for  several  of  the  corresponding  periods  to  have 
the  same  pitch  period.   This  is  synonymous  to  a  crossing  of  pitch  con- 
tours in  a  pitch  period  versus  time  presentation.   Ideally,  the  output 
from  the  systems  should  be  zero,  but  this  result  could  not  be  achieved. 

Figures  5-16,  5-17,  and  5-18  show  the  time  domain  results  of  the 
input  reject  test  for  Shields',  the  adaptive  overload,  and  the  adaptive 
systems  respectively.   From  these  presentations  it  is  extremely  diffi- 
cult to  establish  which  system  is  doing  a  better  job  of  attenuating 
the  input  waveform.   Each  system  has  certain  areas  that  it  does  well 
in,  and  areas  that  it  does  not  do  well  in. 

The  frequency  domain  presentation  does  provide  a  better  picture 
of  how  the  systems  performed  on  this  test.  The  input  signal  spectrum 
has  not  been  included  in  this  series  of  figures  because  it  would  be 
the  same  spectrum  as  the  one  shown  in  figure  5-11.  The  output  spectra 
for  the  Shields',  the  adaptive  overload,  and  the  adaptive  system  are 
presented  in  figures  5-19,  5-20,  and  5-21  respectively.   In  all  systems 
the  general  level  of  attenuation  is  on  the  order  of  five  decibels. 

The  system  shown  in  figure  5-7  may  again  be  used  to  measure  the 
performance.   In  equation  (5.11),  if  the  value  of  a_  is  less  than  one, 
this  corresponds  to  an  attenuation  of  the  input  signal.  Consider 
equation  (5.16)  again.   If  the  value  of  a_  is  less  than  unity,  then, 
the  logarithm  of  a_  wi I  I  be  negative  making  the  value  of  E  positive. 
As  a_  approaches  zero,  the  value  of  E  approaches  infinity.  Therefore, 
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in  evaluating  the  system  performance  for  this  test,  the  magnitude  of 
the  error  function  is  proportional  to  the  amount  of  attenuation  applied, 
The  error  functions  are  given  in  figure  5-22,  and  these  error  functions 
show  several  characteristics.  The  average  level  of  the  three  error 
functions  is  approximately  the  same.  Therefore,  it  may  be  stated  that 
overall  the  adaptive  systems  seem  to  do  no  worse  than  the  system 
proposed  by  Shields  on  the  input  reject  test. 

The  next  test  concerned  the  performance  of  the  systems  in  the 
presence  of  white  noise.   The  white  noise  was  generated  by  an  analog 
noise  source,  lowpass  filtered  at  4.7  ki  lohertz,  and  then,  sampled  at 
10  ki lohertz.  The  gain  of  the  noise  was  selected  so  that  the  signal- 
to-noise  ratio  would  be  much  greater  than  one.  The  noise  waveform  was 
then  added  to  the  signal  as  described  in  figure  5-6  to  produce  the  in- 
put waveform  for  the  enhancement  systems. 

The  conclusions  from  Shields  thesis  were  that  the  system  did  not 
perform  well  in  the  presence  of  white  noise.  The  observation  that  the 
comb  filter  changed  the  wideband  noise  into  noise  that  was  highly 
harmonic  is  exemplified  by  figure  5-23.  The  output  waveform  in  figure 
5-23  is  quite  different  from  the  input  signal.   The  same  statement  may 
be  made  concerning  the  output  of  the  adaptive  overload  system,  which 
is  shown  in  figure  5-24.   The  output  of  the  adaptive  system  as  shown 
in  figure  5-25  is  quite  different  from  the  outputs  of  the  other  two 
systems.  The  input  and  output  waveforms  are  almost  identical,  and 
this  is  quite  encouraging.   The  adaptive  system  passes  the  signal  and 
noise  combination  almost  exactly.   In  the  adaptive  overload  and  in 
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Shields1  system,  the  output  v/aveforrn  is  distorted  immensely  due  to 
the  white  noise.   On  the  other  hand,  the  adaptive  system's  output  is 
noisy,  but  the  distortion  of  the  other  two  systems  is  not  present. 
Figure  5-26  shows  the  noiseless  input  signal  and  the  output  of  the 
adaptive  system  with  seven  coefficients  in  the  filter.   It  can  be  seen 
that  the  effects  of  increasing  the  number  of  coefficients  in  this  case 
is  negl igible. 

The  spectrum  of  the  noise  used  for  this  test  is  shown  in  figure 
5-27,  and  the  spectrum  of  the  input  signal  added  to  the  noise  is  shown 
in  figure  5-28.   Figure  5-29  is  the  spectrum  of  the  output  from  Shields' 
system  for  the  case  of  the  signal  with  additive  white  noise.   It  can 
be  seen  that  the  spectrum  is  composed  of  distinct  harmonic  bands  that 
were  not  as  prevalent  in  the  other  cases.  The  areas  outside  the 
harmonic  bands  have  been  lost  in  the  filtering  process,  therefore, 
causing  distortion  in  the  output.   The  output  spectrum  for  the  adaptive 
overload  system,  figure  5-30,  resembles  the  input  more  closely  than 
Shields'  system,  but  there  is  some  energy  getting  lost  as  the  level  of 
the  spectrum  is  less  than  that  of  the  input's  spectrum.  The  adaptive 
system's  output  spectrum,  figure  5-31,  has  only  minor  differences  from 
the  spectrum  of  the  input.   Figure  5-32  shows  the  error  functions  of 
the  outputs  in  the  same  manner  that  they  have  been  presented  in  the 
past. 

There  is  little  doubt  that  the  adaptive  system  performs  better 
than  the'  other  two  systems  from  the  standpoint  of  this  particular  test. 
The  adaptive  system  provides  an  output  signal  that  may  be  noisy,  but 
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the  other  two  methods  provide  an  output  signal  that  is  both  noisy  and 
distorted. 

The  final  test  was  conducted  for  the  purposes  of  continuity. 
This  test,  as  pointed  out  before,  is  a  combination  of  tests  one  and 
two.   The  system  proposed  by  Shields  is  presented  in  a  time  domain 
representation  in  figure  5-33.   The  adaptive  overload  and  adaptive 
systems  are  shown  in  figures  5-34  and  5-35  respectively.  The  spectra 
of  the  three  systems  are  shown  in  figures  5-36,  5-37,  and  5-38  for 
Shields'  system,  the  adaptive  overload,  and  the  adaptive  systems 
respectively.   Upon  careful  examination  of  these  spectra,  the  level  of 
the  spectrum  from  Shields'  system  is  generally  lower  at  most  frequencies 
than  either  of  the  adaptive  systems  or  the  spectrum  of  the  desired 
speaker  shown  in  figure  5-11.  This  point  seems  to  agree  with  the 
phenomenon  observed  earlier  in  the  input  pass  test.  The  amplitude  of 
the  output  from  Shields'  system  in  the  input  pass  test  was  being  ampli- 
tude modulated,  therefore,  the  waveform  would  have  to  contain  less 
energy  at  some  frequencies  than  an  unmodulated  waveform  of  the  same 
type.   The  error  functions  shown  earlier  have  been  omitted  for  this 
test  because  they  are  almost  identical  and  not  very  benefical  for 
compari  sons. 

In  summary,  this  series  of  tests  was  performed  in  order  to  provide 
some  other  means  of  system  evaluation.   It  should  be  stated  that  the 
fact  the  adaptive  filter  did  not  suppress  the  noise  waveform  better 
than  the  system  proposed  by  Shields  is  not  surprising.  The  main  ad- 
vantage of  the  adaptive  system  is  fidelity,  not  rejection.  Some  of 
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the  phenomena  and  errors  observed  from  these  test  signals  were  help- 
ful In  Implementing  the  systems  that  worked  with  actual  speech.  V/hen 
the  listening  tests  are  performed,  a  more  accurate  comparison  can  pro- 
bably be  presented.  However,  before  these  tests  are  performed,  other 
than  informal  listening  evaluations  and  speech  spectrograms,  these 
results  are  all  that  are  available  for  system  comparison. 


-1  IT- 
CHAPTER  VI 

COMPUTER  I MPLEMENTAT I  ON 

6.1   Description  of  Computer  System 

The  computer  work  for  this  thesis  was  done  on  a  PDP-11/45  digital 
computer.   This  computer  is  a  highly  sophisticated  and  powerful  16-bit 
word  machine  with  32  K  of  core  memory  and  dual  cartridge  disk  drives 
capable  of  holding  1.2  million  words  each. 

Other  peripheral  equipment  included  the  Lab  Peripheral  System,  or 
LPS,  a  storage  oscilloscope,  a  Hewlett  Packard  7004  B  plotter,  and  a 
VT05  Alphanumeric  Display  Terminal.  The  system  was  also  equiped  with 
analog  lowpass  filters,  attenuators,  and  audio  amplifiers. 

The  LPS  is  a  modular,  real-time  subsystem  that  houses  a  12-bit 
analog-to-digital  converter  (A/D) ,  a  programmable  real  time  clock,  and 
a  display  controller  which  includes  two  12-bit  digital-to-analog  con- 
verters (D/A).   With  the  LPS  the  task  of  input  and  output  to  the  com- 
puter was  easily  handled. 

The  LPS  real  time  clock  was  programmed  to  sample  an  input  signal 
at  10  ki  lohertz  after  a  lowpass  filter  with  its  cutoff  frequency  at 
4.7  ki lohertz.   The  computer  was  set  to  sample  two  channels  simultan- 
eously, and  the  speech  signal  along  with  the  glottal  accelerometer 
signal  were  put  onto  the  disk  in  the  form  of  a  data  file.  The  samples 
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were  stored  in  an  interleaving  format  so  that  the  time  reference 

between  the  two  signals  would  be  preserved.   The  program  that  conduct- 
ed the  sampling  was  set  to  sample  up  to  3.2  seconds  of  speech  with  its 
corresponding  accelerometer  signal.   This  amount  of  time  was  chosen  so 
that  most  sentences  could  be  accommodated,  but  with  minor  modifications, 
the  program  could  be  changed  to  allow  any  length  file.  The  program 
that  conducted  the  sampling  also  took  advantage  of  the  direct  memory 
access  (DMA)  option  of  the  A/D  converter,  which  allowed  the  conversions 
to  be  stored  in  memory  at  the  maximum  rates  without  processor  inter- 
vention. 

The  LPS  also  provided  a  capable  display  control.  Along  with  the 
knobs  and  switches  provided,  a  display  program  was  written  that  was 
capable  of  displaying  and  editing  waveforms.  This  feature  was  very 
beneficial  in  correcting  pitch  markings  and  in  viewing  processed  and 
unprocessed  waveforms.   X  and  Y  cursors  with  knob  controls  were  in- 
cluded with  an  LED  readout  of  position  in  order  to  locate  specific 
points  in  the  waveform. 

The  routine  used  for  audio  output  was  driven  from  the  D/A  convert- 
er in  the  LPS.   The  signal  from  the  D/A  converter  was  lowpass  filtered 
at  4.7  kilohertz  and  then  amplified.   At  this  point  the  output  signal 
could  be  recorded  on  magnetic  recording  tape  for  future  use. 

6.2  Shields'  System 

The  filtering  methods  discussed  by  Shields  were  implemented  on  the 
PDP-11/45  using  the  standard  RT-1 1  Fortran.   It  was  decided  that  an 
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implementation  in  fortran  would  be  a  slower  version,  but  after  algo- 
rithms were  perfected,  the  execution  speed  could  be  improved.  Flow- 
charts for  the  systems  described  in  this  section  are  included  in  the 
appendix.  Since  the  A/D  and  D/A  converters  used  12-bits,  the  input 
and  output  data  values  were  in  integer  form  although  the  internal 
multiplications  were  done  in  floating  point  format. 

There  were  two  different  systems  implemented  for  the  comb  filter- 
ing techniques,  and  the  difference  in  these  two  methods  was  the  manner 
in  which  the  unvoiced  segments  were  handled.  A  brief  statement  about 
these  two  methods  seems  necessary  before  the  systems  are  described. 
The  first  method,  the  attenuated  input  method,  stopped  the  comb  fil- 
tering when  an  unvoiced  segment  was  encountered.  The  program  then 
began  to  attenuate  the  input  by  some  specified  constant  until  the 
next  voiced  section  appeared.   The  second  method,  the  "inertial 
filter",  ignored  the  fact  that  an  unvoiced  segment  had  been  encount- 
ered and  continued  to  filter  the  input  using  the  parameters  of  the 
last  known  pitch  period  value.  The  name  inertial  filter  resulted 
from  the  physical  aspects  of  the  processing.   The  filter  could  be 
thought  of  as  moving  through  the  input  samples  with  some  velocity, 
and  without  any  external  forces,  the  filter  and  its  motion  would 
remain  unchanged  as  it  moved  through  an  unvoiced  section.  With  the 
exception  of  the  manner  in  which  unvoiced  speech  sections  were 
handled,  the  two  systems  were  basically  the  same. 

The  filter  coefficients  could  be  derived  from  several  sources. 
A  subroutine  was  written  in  order  to  let  the  source  for  the 
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coefficients  and  the  number  of  coefficients  be  determined  as  the  pro- 
gram commenced  execution.  The  choices  of  sources  included  the  four 
window  functions  that  were  used  in  the  original  system  proposed  by 
Shields.   The  window  functions  were  the  rectangular,  the  Hanning,  the 
Hamming,  and  the  Blackman.  A  second  choice  was  the  coefficients  from 
a  lowpass  filter  designed  by  the  Parks-McClel I  an  Algorithm  that  was 
available  on  the  computer.  The  final  choice  was  one  that  allowed  the 
user  to  type  in  any  coefficients  desired. 

After  the  coefficients  had  been  determined,  the  input  buffer  was 
initialized  with  the  first  section  of  data.  As  mentioned  in  the  pre- 
vious section,  the  waveform  had  been  sampled  using  two  channels.  The 
signal  was  stored  on  one  channel,  and  the  pitch  marks  were  placed  on 
the  second  channel.   This  type  of  implementation  could  be  considered 
wasteful  from  the  standpoint  of  storage,  but  from  the  aspects  of  the 
display,  the  waveform  along  with  pitch  marks  could  be  viewed  very 
easily.   As  the  filtering  commenced  the  center  coefficient  checked 
each  value  on  the  second  channel  for  the  pitch  marks,  and  when  a  mark 
was  encountered,  the  filter  was  changed.  A  pitch  table  file  contained 
the  distance  between  pitch  marks  and  was  used  to  furnish  a  new  pitch 
period  value  after  a  mark  was  encountered.   If  the  value  of  the  pitch 
period  was  greater  than  20  milliseconds,  the  section  was  treated  as 
an  unvoiced  section.   The  unvoiced  speech  procedures  described  above 
were  then  implemented  depending  on  the  system. 

The  pitch  table  was  also  marked  to  denote  areas  that  were  silent. 
These  silent  areas  were  detected  as  mentioned  in  Chapter  III  from  an 


'   •"  -121- 

energy  measurement  over  a  segment  of  data.  These  si  lence  marks  were 
used  only  to  speed  up  the  processing  at  the  beginning  and  ending  of  a 
sentence.   If  a  value  from  the  pitch  table  indicated  a  silent  area,  the 
output  of  the  filter  was  set  to  zero,  and  no  multiplications  were  per- 
formed.  Procedures  for  handling  the  silent  areas  that  appeared  inside 
the  sentence  limits  were  not  thoroughly  investigated,  and  this  could 
be  a  topic  for  future  research. 

From  all  practical  aspects  and  to  the  best  of  the  author's  know- 
ledge, these  systems  were  identical  to  the  ones  described  by  Shields. 
The  execution  time  for  these  systems  was  on  the  order  of  seventy-five 
times  real  time. 

6.3  Notch  Fi Iters 


Another  method  that  was  investigated  for  the  speech  enhancement 

problem  was  a  notch  filter  implementation.  Shields  suggested  that  a 

notch  filter  be  used  with  the  pitch  marks  of  the  unwanted  speaker  to 

33 
reject  the  unwanted  speaker. 

In  the  comb  filter  implementation,  the  impulse  response  of  a  low- 
pass  function  was  used.  Using  the  principle  of  duality  a  notch  filter 
implementation  yields  a  frequency  response  for  one  of  the  possible 
filters  as  shown  in  figure  6-1. 

The  high  pass  prototype  filter  was  designed  by  means  of  the  Parks- 
McClellan  Algorithm.  There  were  several  parameters  that  had  to  be 
varied  in  order  to  decide  what  type  of  prototype  filter  was  needed. 
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Figure  6-1 


These  parameters  were: 

1 .  Fi  Iter  order 

2.  Stop  Band  Width 

3.  Transition  Width 


The  filter  order  of  the  prototype  high  pass  filter  was  limited  to 
less  than  twelfth  order.   This  value  was  determined  from  the  implemen- 
tation.  If  a  number  of  zeros  equal  to  the  pitch  period  were  inserted 
between  coefficients,  then  the  effective  order  of  the  notch  filter 
grows  quite  rapidly.   For  example,  if  an  eleventh  order  prototype 
filter  were  chosen  and  the  pitch  period  were  100  hertz,  then,  the 
order  of  the  notch  filter  would  be  1100th  order.  A  filter  length  of 
this  magnitude  is  on  the  borderline  of  the  quasi -period! c  assumption 
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made  for  the  speech  waveform.   For  the  most  part,  the  high  pass  proto- 
type filter  order  was  between  seventh  and  ninth  order. 

The  second  consideration  was  the  width  of  the  stop  band  in  the 
prototype  filter.  This  parameter  can  be  linked  directly  to  the 
width  of  the  spectrum  of  the  speech  waveform  at  the  various  harmonics 
of  the  fundamental.   If  the  stop  band  is  too  narrow,  then  the  energy 
at  each  harmonic  from  the  unwanted  speaker  will  not  be  sufficiently 
suppressed.   On  the  other  hand,  if  the  stop  band  is  too  wide,  the  notch 
filter  may  suppress  the  desired  speaker  and  cause  distortion.  The 
above  discussion  relates  back  to  the  basic  tradeoff  for  speech  enhance- 
ment. Again,  this  tradeoff  is  the  desired  speaker  distortion  versus 
the  undesired  speaker  attenuation. 

The  final  primary  consideration,  the  transition  width,  was  related 
to  the  first  two  considerations  already  mentioned.  The  transition 
bandwidth  can  not  be  made  sufficiently  low  without  increasing  the  or- 
der of  the  filter.  As  the  transition  width  is  decreased  without  chang- 
ing the  filter  order,  the  ripples  in  the  pass  and  slop  bands  increase. 
Therefore,  the  limit  for  the  transition  width  is  related  to  the  maxi- 
mum ripple  that  can  be  tolerated. 

From  the  prototype  filters  that  were  designed  the  transition 
widths  were  relatively  large  due  to  the  limitations  in  the  filter  order. 

After  the  coefficients  from  the  prototype  high  pass  filter  had 
been  calculated,  these  coefficients  could  be  used  in  the  same  programs 
as  used  in  the  comb  filtering  implementation. 
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6.4  Adaptive  System 

The  adaptive  filtering  systems  were  implemented  much  in  the  same 
manner  as  those  systems  mentioned  in  the  previous  section. 

Since  the  methods  used  to  implement  the  system  that  corrected  for 
the  overload  problem  encompassed  the  methods  used  for  the  system  that 
did  not,  only  the  former  system's  description  will  be  covered. 

The  initialization  was  the  same  as  in  the  previous  description, 
but  an  array  was  used  to  hold  the  information  on  the  spacings  be- 
tween the  coefficients  since  this  system  did  not  use  uniformly  spaced 
coefficients.   Initially  all  spacings  were  set  equal  to  the  first 
voiced  pitch  period.   As  the  filter  moved  through  the  input  data,  and 
pitch  marks  were  detected,  the  spacing  array  was  continually  updated 
to  reflect  the  correct  spacing  for  each  coefficient.   If  an  unvoiced 
area  were  reached,  the  scheme  was  to  perform  the  inverse  operation 
described  above.   The  coefficient  that  moved  into  the  unvoiced  area 
would  retain  the  last  voiced  pitch  period  spacing,  and  the  filter 
continued  in  this  manner  until  all  coefficients  were  clear  of  the 
voiced  area.   In  unvoiced  areas  the  input  was  attenuated  as  in  the 
previous  system. 

Another  array  was  used  to  correct  for  the  overload  problem.  This 
array  contained  the  initial  starting  locations  in  the  buffer  after  the 
filter  had  been  changed.   If  a  coefficient  were  about  to  enter  an  area 
where  overloading  would  occur,  that  coefficient  was  set  to  zero,  and 
the  rest  of  the  coefficients  were  rescaled  so  that  the  output  would 
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have  a  constant  level.   The  array  that  contained  the  starting  locations 
was  used  to  detect  when  an  overload  situation  existed. 

The  rest  of  the  system  involved  only  bookkeeping  operations  that 
were  needed  to  handle  the  rules  used  by  the  system.  The  system  was  a 
little  slower  than  the  one  proposed  by  Shields  and  operated  on  the 
order  of  ninety  times  real  time. 

In  summary  it  should  be  pointed  out  that  the  fortran  implementa- 
tions were  not  designed  for  speed.  Their  basic  design  philosophy 
was  one  of  user  interaction  with  the  programs  to  ensure  correct  pro- 
cessing.  If  speed  had  been  a  factor  in  this  implementation,  several 
items  could  have  been  changed.   First,  the  input  waveform  could  be 
composed  of  a  single  channel,  and  since  the  data  samples  were  only 
12-bits  long,  the  higher  order  bits  could  be  used  for  the  pitch  mark 
information.   With  this  change  in  effect  longer  segments  of  data  could 
be  stored  in  core  memory,  eliminating  several  input/output  operations. 
The  filter  implementation  could  take  into  account  the  even  symmetry 
of  the  FIR  filter  and  eliminate  one  half  of  the  multiplies  per  output 
point.   The  addition  of  assembly  language  subroutines  to  do  the  calcu- 
lational  aspects  of  the  program  may  also  speed  up  the  programming. 
With  the  implementation  of  these  changes,  the  system  begins  to  lose 
the  user  interaction  capabilities  that  have  proven  to  be  very  helpful 
in  this  work. 
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CHAPTER  VI  I 

SPEECH  WAVEFORM  RESULTS 

In  order  to  determine  whether  or  not  the  systems  would  perform  as 
planned  on  real  speech  waveforms,  several  sentences  were  processed  by 
the  various  systems.  A  true  evaluation  of  the  performance  of  the 
three  systems  involves  listening  to  the  processed  outputs  when  speech 
signals  are  used.   These  systems  will  be  tested  in  extensive  listening 
tests  at  a  later  time,  but  from  informal  listening  results  and  spectro- 
graph^ analysis  will  be  discussed  in  this  chapter.   It  should  be 
pointed  out  that  it  is  very  difficult  to  describe  the  results  of  the 
various  systems  with  speech  waveform  inputs.  The  descriptions  and 
comparisons  of  the  systems  will  be  made  with  characteristic  words  that 
are  not  quantitative  in  the  least. 

The  initial  parameters  used  in  processing  the  speech  waveforms 

originated  from  the  information  provided  in  Shields'  thesis.   Shields 

concluded  that  the  optimal  value  for  the  number  of  coefficients  was 

seven.   This  number  of  coefficients  provided  the  best  compromise  for 

the  tradeoff  between  desired  speaker  distortion  and  undesired  speaker 

separation.   He  stated  that  the  Blackman  Window  provided  the  best  set 

of  coefficients  from  the  window  functions  used,  but  he  also  stated 

that  the  differences  in  the  results  between  the  Blackman  and  Hamming 

34 
Windows  were  slight.   With  these  facts  in  mind,  these  parameters 
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could  remain  fixed  while  the  enhancement  systems  were  varied. 

Before  the  results  of  the  tests  similar  to  those  in  Chapter  V  are 
discussed,  a  few  phenomena  characteristic  to  the  individual  systems 
will  be  introduced. 

7.1   Comb  and  Adaptive  Filter  Results  with  Speech  Waveforms 

There  were  several  phenomena  observed  in  the  different  systems 
that  were  not  noticed  in  the  test  signal  section.   In  the  attenuated 
input  method  of  Shields'  system,  an  attenuation  constant  that  was  too 
small  caused  the  output  waveform  to  appear  chopped.  This  chopping 
sensation  resulted  from  unvoiced  sections  that  were  being  attenuated 
by  a  great  amount,  and  this  caused  a  drastic  change  in  the  amplitude 
between  a  voiced  and  an  unvoiced  section.  At  attenuation  constant  of 
0.3  seemed  to  be  the  smallest  value  for  the  constant  that  allowed  the 
chopping  to  be  imperceptible. 

The  inertial  filter  implementation  of  Shields'  system  did  not  pro- 
duce the  chopping  sensation,  but  a  soft  reverberation  was  introduced 
which  was  audible  in  quiet.   This  reverberation  could  be  directly  attrl 
buted  to  the  fact  that  unvoiced  sections  were  processed  by  the  comb 
filter  with  constant  parameters.   Both  the  chopping  sensation  and  the 
soft  reverberation  were  phenomena  that  could  be  attributed  to  the 
unvoiced  speech  segments.  Although  the  time  intervals  for  unvoiced 
sections  in  normal  speech  are  small,  the  phenomena  occurring  in  these 
unvoiced  sections  seemed  to  carry  over  into  the  voiced  sections  and 
disguise  the  actual  results  for  the  voiced  sections. 
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The  same  types  of  tests  performed  on  the  test  signals  in  Chapter  V 
were  implemented  on  the  actual  speech  waveforms.  These  tests  were  eva- 
luated by  informal  listening  and  by  spectrograms.   Spectrograms  will  be 
included  as  a  final  comparison  on  the  methods  that  were  implemented. 
In  these  tests  the  two  systems  implemented  for  Shields'  method  were 
equivalent  since  the  desired  speaker  consisted  of  a  totally  voiced 
waveform. 

For  the  input  pass  test  all  systems  performed  equally  well  for  the 
most  part.   The  most  noticeable  characteristic  of  Shields'  system  was 
a  slight  hoarseness  in  several  words.   The  adaptive  systems  were  char- 
acterized by  a  slight  reverberation.  All  output  waveforms  had  very 
similar  spectrograms  with  only  very  minor  differences  from  the  spectro- 
gram of  the  input  sentence. 

The  input  reject  test  seemed  to  show  that  the  adaptive  systems 
provided  about  the  same  attenuation  of  the  undesired  speaker,  but  the 
outputs  from  the  adaptive  systems  had  more  reverberation  than  those 
from  Shields'  system. 

In  the  white  noise  test  Shields'  system  seemed  to  perform  a  little 

35 
better.  The  "reedy"  sound  mentioned  in  Shields'  thesis    was  present 

in  the  output.   In  the  adaptive  systems  the  output  had  a  "buzzing" 
sound  that  seemed  to  be  more  annoying  than  the  distortion  in  the  output 
from  the  comb  filter.   Again  the  spectrograms  indicate  no  major  differ- 
ences for  this  test. 

The  combined  signal  test  has  been  presented  with  spectrograms,  and 
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the  combined  signal  spectrogram  is  shown  in  figure  7-1  (a).   The  two 
waveforms  were  taken  from  two  different  male  speakers  with  approximately 
the  same  average  pitch.   The  two  signals  were  aligned  in  time  so  that 
both  speakers  would  be  talking  at  the  same  time.  This  alignment  would 
show  how  well  the  filtering  performed  when  both  speakers  were  talking. 
Figures  7-1  (b)  and  (c)  show  the  spectrograms  of  the  two  waveforms 
before  they  were  added,  and  the  spectrogram  in  figure  7-1  (b)  has  been 
designated  as  the  desired  speaker.   Figure  7-1  (d)  is  the  spectrogram 
from  the  output  waveform  of  Shields'  system  with  seven  coefficients, 
and  a  Hamming  Window  as  the  lowpass  prototype.   Finally,  figure  7-1  (e) 
is  a  spectrogram  of  the  adaptive  system's  output  waveform.  The  output 
waveform  in  this  figure  was  produced  with  seven  coefficients  and  a 
Hamming  Window  as  the  lowpass  prototype  also.   Both  figures  7-1  (d)  and 
7-1  (e)  are  very  similar  with  each  one  having  only  minor  differences 
from  the  other.  This  type  of  result  was  normally  encountered  in  other 
sets  of  sentences  that  were  examined. 

As  mentioned  earlier  in  this  section,  the  adaptive  systems  pro- 
duced a  characteristic  buzz.   The  origin  of  this  problem  was  investi- 
gated by  examining  the  input  and  output  waveforms  in  the  areas  where 
the  buzzing  sound  was  most  prominent.  These  waveforms  are  shown  in 
figure  7-2.   In  the  output  waveform  it  can  be  observed  that  the  major 
peak  in  the  input  waveform  is  being  clipped  and  distorted  by  the 
adaptive  system.  This  phenomenon  was  not  observed  in  Chapter  V  when 
working  with  test  signals.   In  the  case  of  the  test  signals,  the  pitch 
period  was  allowed  to  change  very  rapidly  in  some  areas,  and  the 
adaptive  filter  performed  satisfactorily.   It  should  be  pointed  out 
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(a) 
NPUT    WAVEFORM 


(b) 

OUTPUT   WAVEFORM 
ADAPTIVE  FILTER  DISTORTION 


Figure  7-2 

that  in  the  case  of  the  test  signals,  the  impulse  response  of  the 
vocal  tract  did  not  change.   Therefore,  the  distortion  that  is  shown 
in  figure  7-2  (b)  may  be  due  to  the  increased  variation  of  the  impulse 
response  in  short  intervals  of  time.  This  phenomenon,  fortunately, 
does  not  occur  frequently  in  normal  speech,  but  it  may  occur  often 
enough  to  be  an  annoying  problem. 


7.2  Notch  Fi Iter  Results  with  Speech  Waveforms 

The  notch  filters  were  implemented  to  reject  the  unwanted  speaker, 
Several  sentences  were  processed,  and  the  overall  opinion  of  the  notch 
filter  was  that  it  did  not  improve  the  quality  of  the  output  waveforms. 
The  notch  filter  was  implemented  by  itself  and  in  series  with  a  comb 
filter.  The  notch  filter  or  series  combination,  at  best,  performed 
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only  as  well  as  the  single  comb  filter  by  itself.   In  most  cases  the 
results  of  the  notch  filter  or  series  combination  of  notch  and  comb 
were  worse  than  those  of  the  single  comb. 

The  distortion  of  the  desired  speaker  was  higher  in  the  notch  fil- 
ter implementation.   This  increase  in  distortion  was  probably  due  to 
the  wide  transistion  width  that  resulted  from  the  design  specifications, 
The  series  combination  also  had  the  problem  of  "harmonic  overlap". 
The  "harmonic  overlap"  problem  may  be  explained  in  the  fol lowi ng  manner : 
If  the  teeth  of  the  comb  filter  are  set  to  pass  the  desired  speaker 
and  the  notches  of  the  notch  f i  Iter  are  set  to  attenuate  the  unwanted 
speaker,  there  may  be  some  areas  in  the  frequency  domain  when  the 
teeth  and  notches  will  overlap.  This  problem  can  occur  even  when  the 
pitch  of  the  two  voices  are  very  different.   The  overall  result  of 
this  combination  of  notch  and  comb  has  to  be  lower  in  performance  when 
the  "harmonic  overlap"  occurs.  Because  of  the  results  from  several 
processed  sentences  the  notch  filter  implementation  was  not  pursued 
further. 

In  this  chapter  the  descriptions  of  the  speech  enhancement  systems 
that  were  implemented  were  discussed.  This  chapter  was  included  for 
the  purposes  of  showing  that  the  systems  did  in  fact  provide  some  type 
of  speech  enhancement.   These  systems  will  be  tested  in  another  effort, 
and  at  that  time,  some  more  analytical  results  may  be  revealed.  A 
demonstration  tape  for  all  of  the  tests  described  in  this  section  along 
with  other  sentences  was  made  and  may  be  found  in  the  Digital  Signal 
Processing  Group  Library,  M.  I.  T.. 
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CHAPTER  VI  I  I 

RESULTS  AND  OBSERVATIONS 

The  next  two  chapters  will  summarize  the  results  and  conclusions 
of  this  thesis.  Chapter  VIM  will  deal  with  the  overall  results  while 
Chapter  IX  will  cover  the  opinions  and  general  conclusions. 

The  general  goals  of  this  thesis  were  fulfilled.  The  system 
that  was  described  by  Shields  has  been  fully  implemented  on  the  com- 
puter and  was  tested  thoroughly  for  correct  behavior.   As  an  alterna- 
tive to  the  comb  filtering  methods,  an  adaptive  system,  has  been 
proposed,  implemented,  and  tested.  A  limited  amount  of  comparison  was 
performed  with  test  signals  and  actual  speech  waveforms. 

The  adaptive  methods  can  be  summarized  in  the  following  manner: 
The  method  originated  from  the  time-varying  structure  of  voice  pitch. 
The  filter  was  modified  to  conform  to  the  variations  in  pitch  in  or- 
der to  pass  the  desired  waveform  with  as  little  distortion  as  possible. 
It  was  shown  in  Chapter  V  that  the  adaptive  system  was  able  to  pass 
a  time-varying  speech-like  waveform  when  the  impulse  response  of  the 
vocal  tract  was  constant.   From  this  viewpoint,  it  was  shown  that  the 
adaptive  methods  surpassed  the  methods  formulated  by  Shields.  On  the 
other  hand,  it  was  stated  in  Chapter  VII  that  the  adaptive  systems  had 
some  problems  when  the  impulse  response  varied  rapidly  in  short  inter- 
vals of  time.   It  should  also  be  pointed  out  that  the  adaptive  system 
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was  more  complex  and  more  time-consuming  in  the  computer  implementa- 
tion than  the  system  proposed  by  Shields. 

The  adaptive  system  stressed  fidelity  of  the  desired  speaker,  and 
apparently  does  no  worse  on  the  rejection  of  the  unwanted  waveform. 

Finally,  it  was  observed  that  the  notch  filter  and  the  notch  and 
comb  filters  operating  in  series  did  not  provide  any  substantial  im- 
provements in  the  overall  quality  of  the  output  waveform.   In  fact, 
both  of  these  methods  showed  substantial  degradation  in  some  cases. 
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CHAPTER  IX 

CONCLUSIONS 

In  concluding  this  discussion  the  following  items  should  be 
stressed:   First,  the  method  proposed  in  this  thesis,  the  adaptive 
filter,  works  well  in  some  cases  in  the  speech  enhancement  problem. 
In  other  cases,  the  method  has  its  limitations.  The  good  and  bad 
points  of  the  adaptive  system  were  summarized  in  Chapter  VIM.  There 
is  little  doubt  that  the  adaptive  filter  reduces  the  level  of  the  un- 
desired  speaker  in  all  cases.   The  desired  waveform  was  sometimes  dis- 
torted and  occasionally  unintelligible. 

Second,  the  method  used  for  pitch  detection  as  developed  in  Chap- 
ter III  appears  to  be  a  very  good  approach  for  the  pitch  detection  pro- 
blem. The  general  pitch  detection  problem  was  not  considered  in  this 
thesis,  but  the  entire  field  of  automatic  pitch  detection  algorithms 
has  been  the  topic  of  much  research  for  many  years.  This  method  may 
possibly  be  improved  by  first  taking  the  derivative  of  the  glottal 
waveform,  and  then,  using  the  peak  picking  algorithm  on  the  resulting 
waveform.   There  is  still  one  major  drawback  in  this  method:  The 
transition  areas  between  voiced  and  unvoiced  speech  are  not  easi  ly 
marked  whether  automatic  or  manual  processing  is  used.  There  is  a 
general  lack  of  structure  in  these  areas,  and  a  decision  has  to  be 
made  as  to  whether  or  not  a  mark  is  needed  and  where  the  mark  should 
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be  placed.   This  method  of  pitch  detection  also  proved  to  be  proficient 
in  the  automatic  detection  of  silent  areas.  Although  the  maximum  bene- 
fit of  knowing  where  these  were  located  was  not  taken  advantage  of,  the 
silent  area  detection  helped  speed  up  the  processing  at  the  beginning 
and  ending  of  the  sentences.   Overall,  this  method  was  generally  accur- 
ate and  swift.  The  feasibility  of  the  method  may  be  in  question  due 
to  the  fact  that  the  glottal  waveform  is  not  as  accessible  as  the 
speech  waveform. 

Third,  it  can  be  concluded  that  the  notch  filter  used  to  reject 
the  unwanted  speaker  does  not  provide  substantial  overall  improvement 
when  both  attenuation  of  undesired  signals  and  faithful  reproduction 
of  the  desired  signals  are  considered.  This  conclusion  is  not  sur- 
prising, because  the  notch  filter  is  basically  the  dual  of  the  comb 
filter.  Based  on  these  observations,  the  notch  filter  should  not 
be  used  because  it  introduces  the  serious  distortion  of  the  desired 
waveform.   The  conclusion  is  that  it  seems  more  feasible  to  pass  the 
desired  speaker  with  as  little  distortion  as  possible  than  to  atten- 
uate the  unwanted  speaker  to  a  higher  degree. 

A  major  observation  that  should  be  examined  is  the  question  of 
how  effective  these  methods  are  in  speech  enhancement.   It  was  con- 
cluded that  for  normally  mixed  speech  waveforms,  there  is  a  limit  to 
the  amount  of  enhancement  that  can  be  achieved  with  either  of  the 
filtering  methods  discussed.   It  is  believed  that  even  optimally 
designed  filters  for  specific  speakers  or  even  for  specific  sentences, 
would  only  provide  limited  improvement  because  of  inherent  overlap 
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in  the  signal  spectra. 

The  remaining  areas  open  for  improvement  in  this  type  of  approach 
to  the  problem  seem  to  be  in  the  area  of  polishing  the  adaptive  system 
so  that  it  approaches  this  limit.  The  difficulty  in  this  area  lies 
with  the  fact  that  much  more  information  concerning  the  speech  wave- 
form will  have  to  be  revealed  before  the  techniques  can  be  improved. 
Until  that  time,  alternative  methods  should  be  examined  for  the  pro- 
blem of  speech  enhancement. 
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1 

Shields,  Vaden  C.  Jr.,  "Separation  of  Added  Speech  Signals  by 

Digital  Comb  Filtering",  M.I.T.  Master's  Thesis,  M.I.T.,  September, 

1970. 

2  Flanagan,  James  L.,  Speech  Analysis,  Synthesis,  and  Perception, 
Springer-Ver lag,  New  York,  2nd  edition,  1972,  p. 11. 

3  Oppenheim,  A.  V.  and  Schafer,  R.  W. ,   Digital  Signal  Processing, 
Prentice-Hall,  Inc.,  Englewood  Cliffs,  N.J.,   1975,  p. 512. 

4  Fant,  G. ,  Acoustic  Theory  of  Speech  Production,  s 'Gravenhage, 
Mouton  and  Company,  1960,  p. 17. 

5  Shaffer,  H.  L. ,  "Information  Rate  Necessary  to  Trasmit  Pitch 
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Ibid.,  pp. 25-26. 

25 

Parsons,  Thomas,  W.,  "Automatic  Separation  of  Simultaneous 

Speech  of  Two  Talkers",  Paper  Z- 16,  88th  Meeting  of  Acoustical 
Society  of  America,  November  8,  1974,  pp.  1-6. 
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in  Unorthodoxy",  Proceedings  of  the  Institute  of  Electrical  and 
Electronics  Engineers,  Volume  58,  No.  5,  May,  1970,  pp.  707-712. 
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APPENDIX 

A.  1   I ntroduction 

This  appendix  contains  the  actual  computer  programs  that  were 
used  in  implementing  the  systems  described  in  this  thesis.  They  were 
included  as  an  appendix  for  several  reasons.   First,  if  these  systems 
need  to  be  implemented  again,  these  programs  will  provide  a  reasonable 
starting  point.   Second,  if  these  systems  are  to  be  compared  with  other 
systems,  there  will  be  no  question  as  to  whether  or  not  the  systems 
described  in  this  thesis  have  been  correctly  implemented.   Finally, 
these  programs  were  placed  in  the  thesis  instead  of  a  separate  note- 
book in  order  that  the  two  sections  would  not  become  separated. 

These  programs  were  written  in  RT-11  Fortran,  and  all  non  stand- 
ard fortran  routines  will  be  included  for  completeness.   It  should  also 
be  mentioned  that  some  of  the  features  of  RT-11  Fortran  may  not  trans- 
late directly  into  statements  usable  with  another  Fortran  compiler, 
but  for  the  most  part,  the  RT-11  Fortran  is  compatible  with  other 
versions. 

The  PDP-11  structure  stresses  the  use  of  modularity  in  program- 
ming, and  this  concept  has  been  used  to  a  great  extent  in  these 
implementations.   For  each  module  or  program,  a  brief  flowchart  is 
included  with  the  computer  code.  This  will  probably  be  a  helpful  fea- 
ture if  the  programs  in  this  appendix  are  to  be  dissected. 
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Figure  A-1  shows  a  general  layout  of  the  overall  systems  with  the 
modules  shown  as  blocks  and  the  computer  file  names  that  are  used  as 
inputs  and  outputs  to  these  blocks.   Figure  A-1  is  used  in  the  pro- 
cessing of  waveforms  exactly  as  shown  if  there  are  no  pitch  errors  to 
be  corrected.   However,  if  pitch  errors  occur,  then  these  errors  have 
to  be  corrected  as  shown  in  figure  A-2. 

These  computer  programs  were  not  written  with  speed  of  execution 
as  a  primary  concern,  but  with  user  interaction  capability  and  compre- 
hension for  new  users  as  the  primary  concerns. 
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COMPUTER  PROGRAM  INDEX 


Page 

1.  Computer  System  Layout  146 

2.  Correction  of  Pitch  Marks  147 

3.  ADDS  Program  149 

4.  SEPAR  Program  152 

5.  PITV02  Program  158 

6.  PEKS  Program  164 

7.  ZEROS  Program  169 

8.  ADDNS  Program  172 

9.  REVAMP  Program  176 

10.  Shields'  System  Documentation  180 

11.  SHV03  Program  189 

12.  SHV04  Program  197 
13-  COEF  Program  205 
14.  WINDOW  Program  209 
15-  Fl LIN  Program  211 

16.  CALC  Program  213 

17.  Adaptive  System  Documentation  215 

18.  N0NV01  Program  219 

19.  N0NV02  Program  230 

20.  BEG  Program  240 


-145- 


COMPUTER  PROGRAM  INDEX 
(Conti  nued) 


Page 

21.  REBEG  Program  242 

22.  CALC2  Program  244 
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COMPUTER  SYSTEM    LAYOUT 


SPEECH 


PITCH 


A/D 


Sample  waveform 


SIGNAL. XXX 


SEPAR 


S  IGNAL.SSS 


PITCH. SSS 


P I TVQ2 


PI  TCH  .  DAT 


u  u 


ADDNS 


SIGNAL  DAT  PI  TCH. TAB 


Li 


FILTERING   SYSTEM 


OUTPUT.DAT 


FIGURE  A-1 
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CORRECTION  OF  PIT  CH   MARKS 


PITCH,  SSS 


PITCH-DAT 


* I 


ADDNS 


PITCH.MAR 


Create  backup 
file 


i 


P  ITCH. MAS 


Correct   marks 
with    editor 


PITCH. MAS 


SEPAR 


PITC  H  T  AB  PITCH. DAS 


NON.OOO 


Li 


REVAMP 

P 

ITfH.TAR 

PITCH.DAT 

Discard    file 


FIGURE    A-2 
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ADDS 


Initial  ization 


Accept  f  ile  names 


no 


yes 


3QQ 


Quit 


20 


30 


fill  buffer   1 


f  i  II  buffer    2 


40 


Add   buffer  ^ 

Scale 
Resul  t— ^buffer  3 


50 


Write  buffer  3ondisk 


FIGURE    A-3 
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C 
C 

c 
c 
c 

c 
c 
c: 
c 


title:  adds. for  750122 
this  program  adds  two  waveforms  togtthlr  om  the 

ODD    CHANNELS,     ft (M 0    Dl\/irLj:     THE    SUM    BY    2    TO 
PREVEKl     OVEKFLO^. 

INITIALIZATION 

If.lFPtK    DA  lBuF(2b60)  ,  CAjBUF  (2560)  • BAt«UF ( 2560 ) 

C  Al  A    IV/AR*  JVAR«KVAR/2,2»2/ 

WRITE     < 7« 1 ) 

FORMATt"  i»iNAME  OH  FILE  1  <ril.  MAM.  ExT>  t  ) 

k  R I T  E  (  7 ,  2 ) 

FORPAT( •  • ♦ •  • ) 

CALI   ASSIGN  (?t •KKlSSIGNAL.XXX' ,-15) 

URITE(7,3) 

FO-^ATC     '1'NAME    OF    ETlE    2    <FIl  MAlvi.  EXT>«  ) 

WRITE     (7«?) 

CALL  ASSIGN  (3, • KKl ; SIGNAL. SSS» ,-15) 

W  R I T  E  (  7 1  4  ) 

FOR*'AT(«  •••NAME  Of-  CUIPuT  FILL  <H  ILNAM.  EXT>  •  ) 

k  R 1 T  E ( 7,2) 

CALI      ASSIGN    U»  •KKl'SIGNAL.OAT'  ,-15) 
EEFTNE    FILE    ?     (2bl»256,U,IVAR) 
fEFT'E    FliE    3     (2t)lt256»U»  JVAR) 


2" 


3n 
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CEFTNE  FILE"  4  (2bl,256«U,KVAR) 

C  PROCESSING  BEGINS  HERE. 

1"      IF(IVAR#EQ«2F2)  ^0  10  30  0 

C  FILL  PUFFFRS. 

CO  20  lA=ltin 

REAn(2*lVA«)     (UATBUp(IB    +     <  ( T A~3  )     *    256) > « IB=1 « 256 ) 

f  omtinue 
CO    ^n    IC=1,10 

REAn(3»jVAR)  (CATBUpdD  +  <(IC-1)  *  256) ) , Id=1 » 256 ) 
CONTINUE 
C  ACD  T|C  TWO  TOGETHER. 
CO  40  IE=1»2559»2 

BATBUF(IE)=(CATBUF(IF)  +  C «TBUH ( IL ) ) /? 
IZ=IE+1 

batbuf( 1z) -q 
40     comtinue 
c  write  the  output 

r  o  Ro  if  =  1,10 

WRlTf      (tf'KVAR)     (BATEUFdG    +     <<1F-1>     *    256)  )  t  IG=1 1  2t>6) 
COMTINUe 

f-0    TO     10 

ENDFILE  ?     I  TURN  THlCGS  OFF 

ENDFILE  3 

ENDFILE  H 

■  STOP 
END 


50 


if'p 


SEPAR 
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Initialization 

1' 

Accept   file 
names 

yes 


Fill  the  buffer 


Write  onto  disk 
file  1-odd  values 
file  2  even  values 


Quit 


FIGURE  A-4 
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C 
C 

I 

c 

c 

l 

c 

( 

c 

c 

c 

c 

c 

c 

(• 

c 

c 

c 


THi  e:  sepa^.fok 


750121 


THIS  PROGRAM  IS  USFU  10  sfPaRaTl  THE  TWO  CHANNELS  OF  A 
TILE  /jND  WKJLTE  lp,tM  INTO  INDIVIDUAL  FILES  WITH  NO  DATA 

r r,  ruE  even  channels. 

IF  SWITCH(0)  I?  UP  THE  PROGRAM  ALLOWS  NAMED  FILES  TO  BE 

TYPE[  Iv,   IF  SWITCH  (1)  IS  Up  I  HF  OUTPUT  ?  FILE  IS  IM 

FVEN  FORMAT.   IF  SWITCH  (1)  IS  DOWN  THE  OUTPUT  2  FILE  IS 
IM  0{jC  FORMAT. 

SUBROUTINES  USED 

ISSWI  (SFF  REVAMP  pROGP/f  tor  EXPLANATION.) 

INI  TI  A  LIZ  A  1'ION 

INTEGER    DArBl.iFlPbfcO)  ,     CAyBUF(  2560  )  «     BArBUF(2.560) 

CATA  IVAR,  JVARtKV/»R/?,2«2/ 

UR1TL ( 7,1 ) 


URlTE( 7,2) 

FORMATC  •••SWITCH(O)  UOwN  -  DEFAULT  FILES  USED») 

W  R I T  E  (  7 ,  3  ) 

FORMAT (•  »,»TNPUI  -  SIGNAL. XXX*) 


IP 


]  3 


12 
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U  R I TE  <  7  ,  i\) 

FORP"AT(«     »»»pUTpUTS    -SIG|\jAl.  SSS       PITCH.  SSSM 

U  R I T  E  <  7 ,  5  ) 

FORmaK*  »,«SWITCH(o)  UP  -  I/O  FILLS  ^AY  BE  CHOSEN*) 

W  K I T  F (  7,6) 

F(HrM(«  i  ♦  •sWITCH(i)  DOWN-OUTPUT  2  IN  oPD  CONFIG. «) 

U  R  I  T  F  (  7  ,  7  ) 

FOF^aTC1  »,»FWITCH(i)  UP  -  OUTPUT  2  IN  EVEN  CONfIG.*) 

PAUSE  'PAUSE* 

IF  (  ISSWI(0  )  .EQ.D  CO  TO  8 

CALL  ASSIGN  ( 2 ♦• RKl : SIGNAL. XXX • ) 

CALI  ASSIGN  (3t 'KK1:SIGNAL«SSS') 

CALL  ASSIGN  (4,  •KKljPIl  'CH.SSS*  ) 

C:0  JO  15 

in,  R I T  E  (  7  1  9  ) 

FORMATS  «i«MAME  OF  INPUT  FILE  <F ILNaN. EXT>» ) 

URITE< 7, in) 

FORMAT ( •  •  «  »  •  ) 

CALL  ASSIGN  ( 2  ••  KK1  I  SIGNAL  •  SSS  ♦,  -14') 

U  R I T  E  (  7 , 1 1  ) 

FORMATC     »«»MAME    °F    CUl'PUT    FTLt     1     <F  ILL'A^.  LXT>»  ) 

U .  R  1  T  F      (  7  t 1 U  ) 

CALI      ASSIGN     (3, •KK1;SIGNAL«0Q1V ,-14) 

UKITF     (7*12) 

F-ORMAT(«  •♦•NAME  OF  OU1PUT  FTLF  2  <F  ILL'aH.  EXT  >•  ) 

URITF{ 7,lo) 
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CALI  ASSIGN  (<4»  » KKl 'SIGNAL. 002 •  ,-14) 
lb      CEFjnE  FILE  ?  (2bl»256tUiIVAR) 
( EFTNL  FILL  7  l2bl,256»U,  JVAR) 
CEFlNt  FILL  ti  (2bl,256,U,KVAR) 

c  are  wf  Thru? 

37      IF( IVAR.GE.252)  G0  K  300 
C  FILL  THE  PUFFER. 

CO  ?0  IA  =  1,10 

READ(2»IVAR)     (DATFUF(IB    +    ((T^-])    *    256)  )  .,  IB=1 » 256  ) 

20  COMTIMUF. 

C  FILL  1\1E    ARRAYS  -  CHAMNELC+)  =  UDP  POINTS; 

C  CHANNEL (14)  =  EVEN  POINTS. 

CO  30  IC=1 ,2559,2 

CATPIJF(  IC)=DATBUF  (  Ic  ) 

CATRUF ( IC    +    1 )    =    o 
50  CONTINUE 

TO    qn    ID=1,2559,2 

IF     (ISSWI(1).EQ*0)     [TTBUF(Jn)     =    DATBUF(IP    +    1) 

IF    (ISSwl(l)  .EQ.O)    BATBUFdD    +    1)     =    0 

IF(ISSWl(i),rQ.l)     t/TPUF(IU)     =    0 

IF     (ISSwI(1)..EU.l)     tATBUF(JD    +     1)     =    DA  I  Bl'F  (  I  D+l  ) 

40     continue 

c;  WRITE  THr  OUTPUT  ONTO  THt  SPEcIHEO  FILLS. 
CO  r^ 0  IE  =  1,10 
■  u.RlTE(3t  JvAR)  (CAlbLr(iF  +  ((IE-1)  *  256 )  )  « IF  =  1 ,  25b ) 

50      CONTINUE 
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L.  o   r-  n    I  g  =   1,10 
URITEU'KVaR)     (BATEsLFCIH    ♦     ((IG-1)     *    256>  )  » IH=1 ,25b) 

60  comttnue 

bO    TO     1? 
300  EI\iOFlLE    2     !  TURN     H  1 1  C  S    OFF. 

ENDFILE    3 
tN^FlL.L    4 

SI  OP 
F  N  0 


PITV02 
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300  f 


yes 


Quit 


50 


Initial  ization 


Find  first  zero 
crossing 


Calculate  pitch 


GOTO  230 


no 


© 


Set  ZFLAG 


@ 


FIGURE   A-5 
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P1TV02  (CONT.) 
0  © 


yes 


Find    zero 
crossing 


Calculate  pitch 


Bump    pointer 


230 

Bump    IR 

GO   TO    100 

FIGURE  A-5(CONT.) 
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( 

c 
c 
c 
c 
c 
c 
c 
c 
(." 
c 
c 
c 

c 
c 
c 
c 
c 
c 
c 
c 
C 


TITLE:   P1TVH2    7bC?2U 

SUBROUTINES  USED: 
PEKS 

7ERCS 

PARAMFTERS  USED  A*'D  PFFli .IT10N: 

jFLAG:   THIS  is  USED  Id  i.ENOTE  a jvj  OVERFLOW  IM  THE  BUFFEk. 
IT  OCCIJKS  WHEN  A  PEAK  IS  DETECTED  AT  THE  END  OF  a 
RIiFFrR,  AMD  THE  CORRESPONDING  ZERO  CROSSING 
HAS'  MOT  BEEN  FOUjv-D.   WHtw  jF[AG  is  SET  (=1) 
THE  PROGRAM  LOOKS  FCR  A  ZERO  CROSSING  IMMEDIATELY 
AFTER  a  BUFFER  S|  TFT  HAS  BEEM  PERFORMED. 

7FLAC:   THIS  FLAG  I^>  USE[  TO  DENOTE  A  SILENT  ARE  A 

AMP  WHEN  SET  (=1)  A  PITCH  PERIOD  IN  THE  PITCH 
TABLE  IS  GIVEN  t-     NEGATIVE  SIGN  TO  DENOTE  SILENCE. 

TIM  T  I  AL  I  7  AT  JON 

INTEGER    Da  lBUF(?bG0)  »  A(  128^  >•  B<  1280  ),  BL'F<  2560  ) 
•     INTFGEK    SUM(5) ,ZFLAG,     STaR I.     TH,     ENDl 
COMMON    DATtUIC 
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C  A T  A  B  JF  f  T  «  J , ST  AKT/2 5 6 0 * 0  » 3  1 1  ♦  1  / 

I,  AT  A  IVaR, JVaRiKVaR, iRtZFLAG, JF LAG/2, ?,2*1 » 0» 0/ 

CALl  ASSIGN  (2» •HK1;PITCH.SSS* ) 

CALL  ASSIGN  (3i »RK1:PITCH. » AB» ) 

CALL  ASSIGN  <i*i  •KK1;PITCH.UAT»  ) 

CEFIME  FlLt  •>     (2bl  ,256»U,  IVAP) 

CEFINE  FILE  ?  ( 1000 , 2 . U, JVAR) 

pEFTME  FILE  n  ( 2bl ,256 » U, KVAR ) 
10      IF  ( IVAr,e«,252)  Gc  TO  300      jPo  WF  QUIT? 

1=1        IMP,  RESEl  THE  POINTER. 
(.'  r  ILL  THE  BUFFER 

LO  ?p  11=1 , in 

READ(2*IVAR)  (DATBUp(J  t  (dl-l)  *  25& ) ) ♦ J=l« 256) 
2P      CONTINUE 
C  FIND  THE  PEAKS  ANP  ^F^O  cROSSlNbS, 

(ALI   PEKS  ( A,Bi IV AH, SUN) 

CAL!  ZEROS  (A, 6) 
C  CALCUI  Al F  pitch 
C  TS  ThL  0\/EHFlOUi  Fl  AG  SF  To 

TF  (JFLAG.EQ.O)  00  10  100 
C  ThllS  SECTION  PIck^  UP  Th[  ZlKO  CROSSING  FpH  AM  OVE.RFl.OW* 

L=l         !YFS,  FIN[  THE  FIRS!  ZERO  CROSSING, 
bn      jf  (  A(  J)  ,Ne»  0  )  ^°  *>  °  7  0 

C    =  J  +  i 

IF( J.Gr.lpsO)  00  TO  250 

GO  TO  bp 


-1  Go- 
To    emoi=j         ! found  it. 

lu  =  E^Dl-STAKT      !C/\LCULATt  THE  PITCH  PERIOD. 

START=£MD1         !HLSET  THl  START  POINTER. 

TlVDFy  =  ?*J  !CALCULAT£  poitiun. 

EUF( INDEX)  =  4UQ5        i  MARK  THE  PITCH  EPOCH, 

fcU  =  lMDEX/2^6    !T|  TS  SECTION  DETERMINES  THE 

IF(  (KU*?56)  .EQ. INDEX)  KU  =  KM  -  1  !REC«  NO.  THAT 

IREC  =  HU  +  KVAR  ! CORRESPONDS  TO  THE  SPEECH  FH-E, 

IF  ( IREC.pr.P5i)  IREC  =  251 

IF  (ZFLAG.EQ.O)  GC  TC  BO  !lS  THIS  A  SILENT  AREA? 

r  =  -1  *  M  |YEb,  CHANGE  PUCH  PERIOD  TO  A  |\IEG.  VALUE, 

7FLAG  =  0  ICIEAP  Th£  SILENf  FLAG. 
80      |,  RlTE<  3t  JVAR)  h,  IR^C 

I  =  J 

cFL-AG  =  0  ICLEAR  THL  OVlFFLOW  FLaG. 

CO  TO  230 
100     IF(I.GT#12B0)  G°  T0  250  !NO,  A^E  WE  THRU? 

TH  r  SU|^  {iR)      JO^TAIN  THRESHOLD. 

If   (TH.fC.U)  CO  TC  120 

?FLAG  =  1 

CO  TO  ?0  0 
120     IF(B( I) .LT.TH)  GO  T(  200  !MOt  IS  THE  PLAK  VALUE<  TH? 

,„  =  !  "NO,  GO  CHECK  ZERO  CROSSINGS, 

l^r     }F  (A(J).nE.O)  r-0  To  180         !ZLKO  CROSSING? 

u=J+3  !No«  BUMP  THE  COUNTED* 

IF  (J.LE.12BP)  GO  TC  IdO        ! THRU  BUFFER? 
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STAPT  =  START.i  260  «V(_StST0RL  start  location. 


]  t'r 


3^0 


200 
23  n 


2  bo 


J  LAG=1 
GO    TO    ?50 
E  l\i  0 1     =  J 

p=endi-start 

IN3FX=?    *    J 


!  Si  T     r  Ht.    UVfc  RK.OW    FLAG, 
!Go    RLAD     fHE    BUFFER, 

!YLr,  zeRu  crossing  occurred. 

JC/.I  CULATt  THt  PITCH. 

!C/|CULATt  POSITION  FOR  MARK, 
r  U  F  (  I  N  U  F  X  )  =  4u?b  !r/-RK  1  Ht  ^HcH  EPOCH. 
KU  =  IIMdEX/?56    UjICULaU  RECORD  NUMBER 
If   ( (KJ*256) .EQ.I^DeX)  KU  =  KU  -  1 
IREC  =  KU  +  KVAR  !Oh  SPEECH  WAVEFORM. 
IK  (If Ec.Gf.P51)  IREC  =  251 

IF  (ZFLaG.EQ.O)  GO  If  190  ! IS  THIS  A  SILENT  AREA? 
|<  =  -J  i     M         !YlS,  DEMOTE  WITH  NEG.  VALUE. 
7 FLAG  =  0 

URITF<  3«  J\/AR)  M,  IKLC 

START  =  Er.'Ui       IKISET  THE  START  POINTER. 
I  =  I  +  2n        IReSEI  I  ^Y  MINIMUM  PITCH. 
GO  TO  250 
1  =  1  +  1 

IR  =  1/256 

IF  ( ( IK  *  2  56) *NEa  l )  iR  =  IR  +  1 
GO  TO  10  0 

TF  (JFLAG.EQ.O)  STAFT=STAR r-1280 
C  WRITE  OUT  THE  BUFFER  OMc  THE  F1LF  ON  THE  DISK. 
I  0  2  6  0  I A  =  1 , 1 0 
|,R1TF(  4»KVAR)  (RUF(IE  +  U*A  -1 ) *256 ) ) « ID=1 « 256) 
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2fc  p  fOMTTrnJf 

C    CLEAR    THF    GJFFEK,     BUF,     FCR     rHE    NEXT    CALL. 

[ 0    ?70     ]C-\ , ?5b0 

f-  U  F  (  T  C  )  =  U 
27c     CGMTINUE 

r  o   to   l o 

30 r  FMDFlt E    2 

ENDFILE    3 

FlMDFTLE    '4 

si  op 
t  nd 
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PEKS 
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C 
C 

c 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

L 

c 
(• 
c 
c 
r 
c 
c 
c 
c 
{ 
c 
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1  Hi  f  :   PrKS.FOR  7bollO 

THIS  SUBROUTINE  FTNDS  Tl,L  PLAKS  OF  a  WaVEF^RM  THaT  HaVf. 
BEEN  STOPTD  ]M  a  BUFFER  (PATBuF).   'HP  RETURN  CONTAINS  AM 
ARRAY  "Fh  Wljh  rHF  VALUE  OF  THE  pFAK  IF  Op.iE  OCCURRED  AMU 

was  greater  THaN  ZERO,  if  no  plak  occurred  at  that 
point.  Then  a  zero  is  returned, 

PROGRAM  PARA|«ETLR?: 

THE  EI\E*GY  SUMMATION  IN  lO-LOOP  3  HAS  a  ScALe  FACTOR  THAT 
IS  USFD  TO  PREVENT  OVERFLOW.   (In  THIS  CASE  100.  IS  USED 
FOR  TIL  SCALE  FACTOR.) 

THE  ARRAY  "SjjMm  CONTAINS  THL  THRESHOLDS  THAT 
Ir'EPE  LALTULATEn  FROM  THL  F^LRGY  PlF  ASURL  wE|\jT .   THFSE 
THRESHOLDS  O-VL    KETUfcNED  TO  THE.  |WlAlN  PROGRa/  FOR  USE  IN 
PETERFIMTNG  THE  CORRECT  tFAKS, 

THL  TFRESHOLC  SETTING  Tl/J  WAS  LINKED  10  a  PARTICULAR 
MEASl-RE^FNI  CF  LNFRGy  WAS  DETERMINED  PY  EXAMINATION  OF 
SEVERAL  WAVEFORMS  Al\lP  PLcTTilMG  ENERGY  VERSUS   TKUb 
GLOTTAL  PEA«S.   THIS  AREf  OF  THt  PROGRAM  MIGHT  BF  IM- 
PROVE F  3Y  A  MORt  nETAlLE[  STUDY  OF  THE  DATA  FOR  THE 
OPTIMAL  THRESHOLD  VALUES. 
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c: 

t:  PARAMETER  "JrUG"  Ji>  USE[  TO  DEMOTE  a  NEGATIVE 

C  SLOPE.   IFLAG  =  1  -  NEGATIVE  SLUPf 


TFLAG  =  0  -  POSIT  JVt  SL»Pr 


SUBROUTINE  PrK'o     ( E » f  , I VAK, SUM) 

C       INlTlALIZAIION 

INTEGER  FA'BUF(?560) , E( 12B0) , F ( 1280 ) , SUM(5) 

COMMON  FATBUF 

CO  60  I A  =  1,5 
C  PERFORM  AN  ENERGY  SUMMATjON  OVE*  256  POINT?  OF  THE 
c:  GLOTTAL  WAVEFORM. 

S-Uviy  =  o.O 

r;0  3  ib  =  1,511*2 

SUvIX  =  SUfX4  (  (  (FATBUI  (((lA-l)*51?)  +  iB)-20«f8)/300)**2) 

3  C.OiMTINUe 

C  THIS  PA*T  OF  PROG,  ASSIG|v.S  A  TRESHOLD  VALUE  BASED  ON 

r  tpf  ei\erc-y  summation  computed  above. 

IF(SUMX,NE»0)     bO    TO    3  0     !  SUiviX  =  0    MEANS    SILENT    AREA 

7HR    =    0 

GO    TO    55 
]  n  IF(SUMX    -    12P00)     15,40*^0 

15  t  F ( S  U  M  X    -    b  5  n  o )     20,35*35 

2Ti  IF(SUMX    -    1QO0)     25,30*30 

25  •   THR    -    ??0n 

( 0    TO    55 
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30  tH*   =   ?eon 

f 0    to    55 
3b  iH*    ~    5C00 

GO    TO     55 
4  0  IF     (  S U vlX    -    1  <■  0 0 0  )     4 b ,  5 0  »  50 

M-b  THR    =    StfOO 

GO    to    55 
bf  THR    =    3800 

bb  SU^I     (IA)     =    thR 

br  CONTINUE 

C    FIND    THE    PEAKS    hOP    THE    SeG^ENI     Of    DATA. 

10    10  0    KH  =  1»1280 

V  (KTT)=0 
inn     rONTlNUE 

K  =  l 

K2  =  1 
130     IF(FATBUF(»U?)  -  FATBUF(K)J  l*f0,150»l&0 
IMC      IF  LAG  =  1 

GO  TO  170 
1 bn     K  =  K  *  2 

GO  TO  UO    . 
16  0      1FLAG  =  0 
170     LO  pio  K=3,255y,2 

F  =  K  -  2 
.  IF  (FAl'BUF(K)  .GT.FATF'UF(P)  »  GO  TO  20  0 

IF  (IFLAG.EQ.l)  GO  TO  205 
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20(1 
?05 
210 


IF  (FATbUF(K).LT.2C46)  GO  10  205 
F(<?)  =  FATBUF(K-2)      IFUUND  THL  PFAk 
IFLAG  =  1 
GO  TO  2o5 

i  F  L  ,a  0  =  0 

K2  =  K?  +  i 

continje 
re:  turn 

END 


ZEROS 
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(enter  ) 

»p 

Initialization 

if 

Zero  array 
C 

< 

' 

Set     up   IFLAG 
for  initial  slope 

V 

Mark  all    zero 
crossings  with 
negative  slopes 

> 

(pet 

r 

URN  ) 
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THi  r:   Zeko^.fok 


750110 


('  THIS  SUBROUTINE  FTMbS  ml  ZERO  CROSS jroS  OF  f\    WAVt- 

C  FORM  STOPED  jN  A  PUFFER,   jf-  A  ^EKO  CHoSSyMO  OCCURS  AT 

L  SOME  FOlMTi  THaT  P01M1  Is  ^aRKED  WITH  A  ONE  IN  THE  ARRAY 

C  11  c » .   IF  NOT,  A  zrRO  IS  LEFT  IN  The  ARRAY. 

r 

C  PARAMETERS  USED! 

C  "IFL/,G"  TS  JSED  IM  THIS  SUBROUTINE  TO  DLNOTh  AM 

C  A  VAliE  OF  THF  UATA  THAI  jS  NEGATTVF. 

0  I.E.       IFUAG  r  1   DATA  VALUE  IS  NEGATIVE. 


i  n 


TFLaG  =  0   OATA  VALUE  IS  POSITIVE. 

SUBROUTINE  ZEROS  (C,D 

INITIALIZATION 

INTEGEK  CA  rBl»F{  2560  )  ,C(12BO)  ,D(  1280) 

COMMON  CATBUF 

ro  in  1 1  i  =  if i2ao 
c  (  itt)  =  n 
CONTINUE 
N  =  l 

K2  =  1 
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2^      IF ( rAl HUh ( r)  -  ?0MC)  30,40*50 
3  0       I F  L  A  G  =  0 

60  TO  bo 
C  T'ANK  A  ZFRO  CROSSTMG. 
fO      C(M2)  =  1 

h    -    N  +  2 

1\2  =  N?  +  i 

IF  (N»Ge,2560)  GO  Tc  20  0 

GO  TO  25 
50       iFLAG  =  1 

60      IF  (CATbUF(N)  *  20«+6)  /0,1n,po 
70      IF  (IFLAG.EQ.O)  GO  TO  ?5 

C ( M?)  =  1 

rFLAG  =  0 
75      |\  =  N  +  2 

l\2  =  N2  +  ] 
•  IF  (N.GE.2560)  GO  TL  ?00 

GO  TO  fee 
80      IF  (IFLAG.EQ.l)  GO  TO  75 

( :(  Mp)  =  i 

I F  L  A  G  =  1 
GO  TO  7^ 
20 p     RETURN 
F  N  D 
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title:  addns.fdk 


7  5  0  1 20 


THIS  fR0p,Pa^  ADDS  TwO  FIlFS  TOGETHER  WITHOUT  A  STALE 
FACTOR.   TT  IS  ytrtif    TO  Lt  USED  ON  Am  ODD  -  EVEN  FORMAT. 
(I.E.  OMr  FILE  Oqp  FORMAT*  THE  OTHER  EVEN).   IF  SWITCH(U) 

is  up  thr  program  allows  the  in^ut  amd  "output  files  to 
pf  napes,  jf  switch  <oj  is  dowivj,  them  the  program 

ASSUMES  THE  CEFAi'L1  NAMES  THAT  ARE  USED  IN  SERIFS  WITH 
THE  CCMPLETE  SHIELD^  SYSif1^. 

SURROLTINES  USED 

1SSWI  (SEF  REVAMP  FOR  EXpL.ANANTlON.  ) 

INITIALIZATION 

INTEGER  DATBUF(2560) , CATBUF(2560 ) ♦ BAjBUFC 2560 ) 

TATA  I^AR.  JV/\R»KVAR/2»2«2/ 

URITE( 7,1) 

FOR-MAT(«     tttpROGRAM    CONTROLS    -    -    -    SyHCH    (0)') 

URITE     (7«2) 

FORMAT (•  • ••SWITCH* C)  DOWN  -  DEFAULT  FILES  USED» ) 

URlTF  (7«3) 

FORMAT! •  »»»JNPUT'-  SIGNAL. SSS») 

•  UR1TF( 7,q ) 
FORMAT («  »t»0UTPU7  -  SIGNAL. PAT • ) 
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WRITE(7,5) 

5  FORMAT(«  i»tcwiTCH(o)  UP  -  I/O  FILLS  MAY  BE  CHOSEN') 
I  AJSF  ' P  A 1 1 S  E  * 

IF  (  I  SSI,:  I  (  0)  .EH-1*  GC  TO  6 
CALL  ASSIGN  ( 2 f' «K1 'SIGNAL. SSS* ) 
CALL  ASSIGN  (3,»KKl«PITCH.UAT») 
CAL!  ASSIGN  (iff  »KK1:SIGNAL»DATM 
p  0  TO  1 o 

6  WRITE (7,7) 

7  formaTc    h»name   cf    input    FILE    <FILNaM.EXT>»  ) 
WRITE (7,8) 

8  FORMAT < '     • t •     » ) 

CALI   ASSIGN  (2,  •RKl  ;MGNAL.  SSS»  ,-14) 
CALL  ASSIGN  (3,»RKl:PITCH.DATt  J 

WRITE  <7f9) 

9  FORPATf  •♦•f'AHF  OF  OUTPUT  FILL  <h"ILNAM.  EXT>»  ) 

WRITE ( 7,6) 

(ALI.  ASSIGN  ( if  ,•  KKl :  SIGNAL.  DAT ', -14 ) 

10  DEFINE  F 1 1  F  ?  (2bl,256f U, IVAR) 
DEFINE  FILE  *  ( 2bl ,256f Uf JVAR ) 
DEFINE  FILE  M  I  ?b3  ii^b6f  U,  KVAR  > 

(  PROCESSING  BEGINS  HERE, 

13      iF( IVAK,EQ.2^2)  GO  TO  30  0 

C  FILL  BUFFERS. 

no  ?n  iA=itio 

PEAP(?»IVAK)     (DATBUf(IB    +     ((IA-1)     *    256 ) ) , IB=1 f 256 ) 
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20  CONTINUE 

do  * n   i c = 1 » i r 

REAP(3'JVAR)     (C/\TBUp(ID    +     ((IC-3)     *    256) ) » ID=1»256> 
30  CONTINUE 

c    ADD    II-!.    TwO    TOGETHER. 

CO  HO  U=l  ♦  P?5^»2 

FAtRUF(IE)=D/»TBUh(It)  +  CATBUF(IE) 

BATpUf ( IE+l )     =    CATBLF < IE+1 » 
ifO  CONTINUE 

C    WRITE     7HF    OUTPUT 

CO    ^0     IF     =    1,10 

URlTF  U'KtfAP)  (BATfcUFClG  +  ((IF-D  *  256) ) t IG=1 »2&6) 
hn      CONTINUE 

r-o  to  11 

300     ENDFILE  2  ITl.'RN  THIf.GS  OFF* 
ENDFILE  3 
EMDFILE  4 
STOP 
E  N  D 


REVAMP 
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lllir:   REVAMP. FUR 
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subroutines  used: 

isswi  -  this  subroutine  is  writ i en  in  assembly 

language  and  allows  fop  the  display  register  switches  tu 

pE  USED  AS  PROGRAM  CONTROLS. 

FOR  SUITCH  ( 1 )  UP 

INPU1  AMP  OUTPU"  FILES  h/ Y  BE  SPECIFIED. 

FOP  St  ITCH  (0)  DOU'N 

INPUT  -  SIGNAL. DAS 

OUTPUTS  -  1.   SIGNAL. UAT   2.   PITCH. "TAB. 

Th'IS  FROpRA^  IS  USED  TO  REVAMP  A  SIGNAL  FILL  THAT  HAS 
BEEN  EDlTFD.   THE  OuP  VAlUES  OF  THE  SIGNAL  ARE  PRESERVED 

in  this  transformation,  tpe  eve'M  values  or  hitch  marks  are 

ZEROES  JNLESS  THEY  APL  Cl  TO  i|0/5.   AT  THE  SAME  TIME 
THE  DISTANCE  BETWEEN  PITCH  MARKS  IS  MEASURED  AND 
STORED  IN  PITCH. TAB  FOR  USAGE  IN  OTHER  PROGRAMS. 

INITIALIZATION 

INTEGER  DA  I'BUF  (  256C  )  ,BATBUH  2560)  t  START  »A(  256) 

CATA  STARTiA/2»256  *  0/ 

LATA  IVAR, JVAR»KVAR/2»2i2/ 
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C.    PEVIEU    PREVIOUS    PTTCH.TAe    FILE 

CALI      ASSIGN    (2i »KK1 -PITCH. TAB* ) 
DEFINE    FILE    ?     (1000,2»U*JVAR) 

1  REAP     (2»  J\/AR,Ei\iD    =    ^)     M,     IKEC 
IF     (M.LT.O)     A(1REC)     =1 

GO    TO    1 

2  ENDFILE    2 

C  MOW  RElMlTlALlZt. 

kRJTF  (7«3) 

3  FORMAT  {'  >,*SWTTCH(0)  DOWN  -  DEFAULT  FILES  USED.') 

U  R I T  E  (  7  » '4  ) 
1       FO^MATC  »,'TNPUT  -  SIGNAL. PAS;  OUTPUT  -  SIGNAL.DAIM 
UUTr  (7.5) 

5  FORMAT  (•  tvtsWlTCH(0)  Up  -  I/O  FILES  P"AY  BE  CHoSLimM 
PAUSE  'PAUSE1 

IF  (  ISSWI(O)  .EQ.D  CO  TO  6 
CALI  ASSIGN  (2» ,KK1:SIGNAL.DAS' ) 
CALL  ASSIGN  ( 3 i * KK1 ; SIGNAL . DAT ' ) 
CALL  ASSIGN  (4i *HK1 :PI rCH. I AP» ) 
GO  TO  10 

6  IvRlTF  (7*7) 

7  FOr<'»AT     (i     t,i|y|APIl     0|      INPUT    FTLL    <F- ILNAM,  EXT>«  ) 

k  R 1 T  E     (  7 «  8  ) 
P  FORMAT    (••,»•) 

CALI  ASSIGN  (2»  • KK1 : SIGNAL. DAS'  ,-m) 
k R I T  E  <7»9) 
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9  FORMAT  (•  »,»NAr*E.  Cf  OUTPU!  Flip  <KILNAM.  EXT>«  ) 

I,,  R  1  T  F"  (  7  »  P  ) 

CALL  ASSIGN  (  3  »•  RK1 •  SIGNAL.  DAT  ',  -l«f ) 

CALF  ASSIGN  ( tf  »  •  KK1 ;  PITCH.  '  AB»  ) 

10  ft  FINE  FILE  ?  ( 251 ,256 , U, I VAR ) 
CEFlNE  F  iLt:  ?  (2bl,256iU«  JVAR) 
DEFINE  FILE  4  ( 1000 ,? , U» KVAR) 
v-VAp  r  2 

C    /\RE     V  E"     TlifvU? 

11  IF     (I VAR, EG. 252)     GO    TO    30  0 

C  FILL  THE  BUFFER  WITH  DAT/  FROM  INPUT. 

L  0  PC  1  =  1,10 

REAn<2»lVAR)  (DATBUF(J  *  (  (  1-1)  *256)  )  «  J=D  256) 
?r      CONTINUE 

C  TRANSFER  THE  ODD  SAMPLES  AS  ThEY  ARE,   TRANSFER  THE  FVEN 
C  SAMPLFS  AS  ^ERO  UNLESS  T|-FY  EQUAL  M095.   ALSO  THE 
C  PITCH. TAP  FILE  IS  BEING  WRITTEN  OUT  IN  THIS  LOOP. 

LO  FO  H  =  1,2559,2 

FATni.'Ft  II)  -    DATHL'F  (ID 

IF     (DATBUFi  1 1  +  1)  .GE.4075)     t>0    TO    25 

FATrUF ( II    +     1)     =    0 

C-  0    TO    50 
25  FATRUF ( II    +    1 J     =    4055 

[v     =     (  (  Il  +  l)     -    START )/2 

,.U    =     (  U     +     D/256 

1F(  (  JU*?5c,).rQ.  (  13+1  )  )     JU    =    JU    -1 


mi 


50 
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,KtC  =  jU  +  jVAP  K./LCULATL  REC.  NO. 

C  CHECK  FOR  SILENJ  />RLAS 

If-  (M.LI  .  ?no  )  ^n  TO  no 

jF  (  A<  JREC)  .NE.D  &C  TO  40 

p  =  M  *  -1  !AK[A  IS  A  SILENT  ONE, 

URITE  U'KVAR)  P^  URFC 

START  =  II  +  1 

CONTINUE 
C  WRITE  THF  OUTPUT  INTO  SIGNAL. DAI  FILE 

LO  70  12  =  1,10 

L,RITF(  3t  JwAR)  (BATBUFII3  +  ((  12-1  )  *25&  )),  13  =  1 »  25&) 
70       fOMTIMUE 
C  RFSE1  START 

START  =  START  -  256c 

GO  TO  1] 

ENDFILE  2  ! TURN  THINGS  OFF. 

ENDFILE  3 

eNdfile  q 

STOP 

[NO 


30  p 


-180- 


ITTLL;   SHIELD. OOc 
DOCUMENTATION  FOR  SHV03  AfviD  SHV04 


AHFTRAC 1 


these  pkogr/>v^  impli  pfut  ihf  systems  designed  and 
mwulated  py  vaden  shields  in  an  mit  master's  thesis  fkom 
1970.  the  purpose  of  thest  systems  was  to  separate  a  signal 

FRCP  NOlSE  (USUAL«-Y  A^OTHth  SPE/\KtR).  THLRE  ARE  TWO  SYSTEMS 
IMPLEMEM  Ep,  AND  THESE  UlMfR  IN  THE  MANNER  Of-  UNyOlLEU 
SPFECH  PROCESSING.  SHVOj  USES  IMF  ATTENUATED  INPUT  METHOD 
WHjlE  SEVDu  USfS  IHF  lNLRTiAL  FILTERING  METHOD. 

SUBROUTINES  USED 


1.  WINDOW  (!\'TYPL,K)  .  THIS  SUBROUTINE  ALLOWS  FOR  I  HE 
USfP  TO  PRFSCRIPE  THE  WINDOW  FUNCrION  DESIRED  IN  THE  PROGRAM 
AlVT  TO  pESlGNATE  IHF  VALUE  cF  K  FOR  THE  FlLTEK  LENGTH, 


2.  COEF  (NTYPFtK, A,L )  -  THIS  SUBROUTINE  CALCULATES  I  HE 
COEFFICIENTS  Of  ThE  FILTER  PRESCRIBED  IN  THE  WINDOW 
SUPROUTINE  AND  SyORrS  tHESl  VALUES  IN  AN  ARRAY  "A"  THAT  LAN 
Hf-  OF  MAXIMUM  DIMENSION  01  15. 

OTHEH  PROGRAMS  I ' S (-  p  I1^  CONJUNCTION  WITH  SHIELD, 
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1.   RFSA^p  -  THIS  PROGRAM  SA«iP|_hS  THE   AnAlOG   WAvEFuRn 

a"»  iu  khz.  on  r r c ' J  °F  ?  channels  (4  and  14)  and  places  its 

UUTPUT  CN  FILL  SIGNAL. XXX  ,    THE   CHANNELS   ARE   STORED   IN 

INTERLEAVING  FORMAT. 

CHANNFL  <4)  -  pPu  SAMpi  E  VALUES.   (ODD  FORMAT) 
CHANNFL  (10)  -  e^FN  SArPLE  VALUES.   (EVEN  FORMAT) 


2.   SfPM  -  I'HIS  PKOGFaM  SEPARA'fcS  THE  TWO  CHANNELS   ON 
ONE  FILE  Ann  pl  AC^S  EACH  Of  a  SEPARATE  FILE. 

infui  -  signal. xxx. 

OUTPUTS  -  CHANNEL  (H)  -  SIGNAL. SSS 
CHANNEL  U4)  -  PITCH, SSS 


BOTH  APE  ASSIGNED  IHE  FORMaT   °F  A   CHANNEL   (4)  ONLY 

Mir.    (I.E.    ODD   SAMPLES   -   DAT A  J  EVEN   SAMPLES  -  0  BY 

DEFAULT.)  TF  DESjREP  THE  SecOND  OUTPUT  CAN  pE  PUT  JNTO  EVEN 

H!P  AT. 


3.  =>JTV0?  -  THIS  PROGRAM  DMFC1S  PITCH  FROM  'HE 
PITCH. SSS  FILE.  TH£Rt  /  RE  TWO  OUTPUT  FILES  FROM  THIS 
PROGRAM.  FIRST,  PITCH. I AE  IS  A  TABULAR  FILE  OF  THE  VALUES 
OF   TH[    PITCH   PLRTOu.   If  IS  IS  USFD  AS  AN  IMPyT  TO  SHIELD. 

second,  pitch, pat  is  a  fill   ihat  ts  used  ro  mark  piich 

PtRlODS   IN  A  WAVt'FPRh  FILt.   P1TLH.DAT  HAS  TWO  VALUES  ONLY* 
(Ll  OR  4095).   THE  VALUE  4l',r  IS  USED  10  Ufc  NOTE   AN   ENnPOiNT 
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OF   A   FlTcH   PERIOD.    THIS   IS   ijSKn  AS  AM  INPuT  TO  ADDNS. 
PITCH.   DAT  IS  Im  EVEN  FORMAT. 

'4.     apdns  -     this     prc.cram     adds     together     two     fills* 

signal. XXX  AND  PlTcH. UAT.  SINCE  UNt  FILL*  SlbNAL. XXX*  IS  IN 
Oir,  FO[^I  ftNO  PITCH. DAT  IS  IM  EVEN,  » Hf  CHa^ELS  aRL  KLPT 
SEPARATE  BUT  ^UT  ifoTO  ONI  p  TL  E  CALLED  SIGNAL. HAT, 

inputs  to  shield 

1.   SlRiMAI  .uAT  -  FROM  AUDNS 
?.   F ITCH. TAP  -  F  POi^  PI  rV02 

OUTPUTS  FROM  SF<1FlD 

1.   OUTPUT. uAI 

Al  RORITFM 


1.  Tmf  PROgKa1""  FTKS1  FROVIDLS  THE  USER  WITH  THE  CHOICE 
OF  TYPE  OF  WlMDOW  AMD  THE  f-ARAMETLR  K.  (2K  +  1  COEFFICIENTS 
AHE  USEr  IN  TH[  FILTER) • 


2.  the  ^RORK^r  then  calculates  thl  coefficients  ^ok 

TIF  FILTER  AMD  STORES  THFf";  JM  ARRnY  "A". 
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5.        T|tr  MtXT  SFCTTOM  /I  LOWS  i-OP  THfT   INI  1  I ALIZATION   OF 

THF  LOGICAl  UMUS  T f ' A  '  WILL  RL  USLD  IN  THIS  IMPLEMENTATION, 

A.  LOGICAL  Uf„TT  ?  -  USED  FOK  THE  SIGNAL  FILE. 

B.  LOGICAL  UNJTT  3  -  USEU  FOK  THt  PITCH  FILE- 

C.  LOGICAL  UNIT  H     -  USED  FOK  THF  OUTPUT  FILL. 


NOTE:  THL  INPUT  aNU  Ol.IPUT  SIGNAL  HILE-S,  OM  LOGICAL 
UNITS  g  AND  i\  A«E  Dl\/IDtL  INTO  KECOHqS  OF  LENGTH  256,  I  HE 
PITCH  FILE  IS  DIVIDED  INTO  RECORDS  nF  LENGTH  IWO. 


4.1HE  PUFFER  IS  1'HtN  SLT  Up  INITIALLY  BY  LEAVING  IK 
KEcOROS  FULL  Op  7EKOS  (IR  IS  COMPUTED  IN  THE  FOLLOWING 
MANNER:  {l/o  Op  Ihf  WINDOW  LENGTH)  I IM  SAMPLES  *  1/^56 
(HFCORL/SA^pLEs)  +  1  IF  F  THE  NUMBER  OF  SAMPLES  IS  NO  I  A 
lvU|  TIPLf  OF  ?5f,.)  THEN  THL  SIGNaL  IS  IaKlN  fKOM  THE  FILE  AND 
PI  fCEO  INTO  THE  rUFFER.  THl  NUMBER  OF  RECORDS  RFAD  IS 
VAPIABLF  DEPENDING  ON  THE  flRST  VALUE  OF  PITCH  AND  ON  K. 


5.   IN  THIS  PROGRAM  IfF  VALUt  I  IS  USED   AS   A  BUFFER 

POINTER.    THE   VALUE   1   FOlNTS  '0  THF  CENTER  OF  THE  FlL'EK 

THpT  IS  BEING  USED.   i  IS  NOW  S£T  UP  10  BE  THL  FIRST  SAMPLE 
IN  THE  CATA  SECTION  OF  THE  F-UFFER. 


6.        THF     VALlJL     OF     riTCj      IS    READ       AND       THROUGH       THE       SAME 
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SERIES   OF   MEaSUKEkENTS   /,S   FOUimH   I|\!   THE  I  M  T  T  I  AL  IZAI  T  I  ON 

PROCESS,  IMF  ALGORITHM  CHECKS  TO  SEE  IF  THE   t XTRLNE   VaLUeS 

NEEDED  IN  THE  plFFERENCt  L&Ua'TION  APf  WITHIN  THE  BUFFER, 

7.   IF  THfY  /\RF  (MOT,  THE  BUFFER  IS  SHIFTED;   THE   VALUE 
OF  I  REPOSITIONED  AND  MtV  i  aTa  IS  READ  IN  FORM  THE  FILf. 


3#  MOW  THE  CALCULATIONS  CAN  BE  PERFONED.  I  HE 
UlFFERtfCE  EQUATION  *S  I  [  ,'  L  F^|E,\Tt.U  Ar'D  An  OUTPUT  POINT  Y  IS 
CALCULATED  AND  WRITTEN  ONTO  THE  OUTPUT  FILE. 


9.  TfF  PROCt-SS  CONTINUES  UN»IL  1  h'E  PITCH  FLAG  IS  SET 
HY  READING  THE  MftKK  IN  TEL  INPUT  FILE.  A  NEW  VaLUf  OF  PIlCH 
PERIOD  IS  THEM  READ,  AND  T|F  FILTER  IS  MODIFIED  FOR  THIS 
CHANGE, 


10.  A  SLIDE  OPERATION  OCCURS  WHEN  THE  VALUE  OF  I  IS 
TOO  SFALL  TO  ACCU,v]wODATE  li  F  CURRENT  FILTER.  THE  BUFpER  IS 
REPOSITIONED  AND  1  PFSFT  Sc  THAT  ALL  INPUT  VALUES  NEEDED  BY 
IHF  FILTER  ARE  PRESFNTLY  If  THE  BUFFER, 


11.  THE  PROCESSING  C  (:T  TINUES  UNTIL  A  ?ILENI  AREA  IS 
DETERMINED,  A^d  THEN  AFTER  FILING  THE  OUTPUT  FILE  WITH  ZFKOS 
THE  PROGRAM  HAlTs. 


SHV03 
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Initialization 


Accept   CON, 

Attenuation 
Constant 


CALL  WINDOW 


CALL   COEF 


20 


READ    M     from 
PITCH. TAB 


50      r 


Take 

Measurements 


SQL 


Fill   buffer 

initially 


100       | 


Clear   PFLAG 


T 


FIGURE    A-10 
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SHV03 


Clear    ZFLAG 


READ    M   from 
PITCH. TAB 


no    < 

'   M    <    0  ? 

\ves 

V 

'' 

Set     ZFLAG 

130 


£ 


no 


Clear    LARGE 


yes 


20  j 


Set      LARGE 


Take 
Measurements 


GO  T0 100 


yes 


GOTO  100 


Assign 
(!25toKXX 


1)0- 


GO  TO  200 


READ    in  a  new 
buffer 


end   of 
f  i  le 


129  j 


© 


1=1 


SOTO  300 


Shift  buffer 


£ 


Reset    I 


01 

FIGURE    A-10     (CON'T.) 
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SHV03 
© 


© 


Does 


DO    ^buffer 


need   \ 

yes 

170 

Slide   buffer 

u 

Reset  I 

- 

190 


200 


T 


Assign  140   to 
KXX 


Set    up  location 
for   output 


no 


yes 


Implement 
difference  equation 
result  — »-Y 


yes 


Add     zeros  at 
beginning 


no 


no 


OUTPUT  =  CON* 
INPUT 


OUTPUT:  2048 


OUTPUT;  Y 


FIGURE    A-10     (CON'T.) 
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SHV03 


270 

Output  for    2  nd 
chan-   2048 


no 


yes 


WRITE  OUTPUT 
cnto    disk 


End 


GO  TO  KXXC125, 
140) 


400. 


Quit 


300 


Add     zeros   to 
end   of  buffer 


no 


yes 


GO    TO   159 


GO  TO     129 


FIGURE    A-10    (CON'T.) 
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C 
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title:  shv03     author:  r.  fka7Ier 

adapted  fkcm  thesis  py  v.  shields 

this  frdgRam  uses  the  method  of  shields  to  perform 
the  separation  of  a  speaker  a|\d  noise  by  means 

OF  A  DIGITAL  C 0 fl B  Fit  TEK.   IHIS  PROGRAM  USES  THE 
ATTENLAlION  TETHOn  FOR  U| VOICED  SECTIONS, 

SUBROUTINES  USED! 
WINDOfc  (MTYPFiK) 
TOFF  (NTYPEiKt AtL) 

F ARAhf Tlps  used: 
sprays : 

1.  DATBUF  -  CONTAIKS  THE  INPUT  SAMPLE  VALUES. 

2.  A  -  CONTAINS  THE  COEFFICIENTS, 

3.  OUTPUT  -  TEMPORARY  BUFi-EK  FOR  OUTPUT  VALUES. 
VAPIAFLES! 

1.   SFLAG  -  USED  TO  IN1TIALI7E  OUTPUT  EILE 
SO  THrr  STEP  DISCONTINUITY  DOES  MOT 
OCCUR. 

?.   7FLAG  -  USED  TO  DENOT  BEG.  OF  SENTENCE, 

3.   PFLaG  -USED  10  |  f NOTE  NEW  PITCH  PERIOD. 

ft,   LAKGE  -  I'SLD  TO  rENOlE  UNVOICED  SECTIONS. 

5.   CON  -  IHF  ATTH-LATION  CONSTANT  F0«  UNVOICED 
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i  m  t  i  ^  l  i  ?  a  r  I  n  r ; 


C  SECTIONS. 

r 
c 
c 

LlviF^S10^    AM  5) 

INTEGER  DATBl'F(*K>96)  , START* OUTPUT  (H5&)  *ZFLAG»SFLAG 

LOGICAL  PFLAR»LAKGE 

[ATA  DAtEiji-  »  SFL AG*  OUTPUT/ l+U 96*0  ,1«.256*0/ 

L/l A  IVAR. JVAR* KVAR/2*2«2/ 
C  CALL  THE  FOLLOWING  SUBROUTINES  HOW  INITIALIZATION, 

W  R I T  E  (  7 , 1 ) 
3        FO*PAT(*  •« 'VALUE  Of  AITEN.  CONSTANT?  <1»<  F3.2>») 

ACCEPT  2,  CON 
2       FORMAT (F3. 2) 

CALI  ASSIGN(?,  •RKi:siGNAL.DAT»  )  jOuTPuT  FROM  ADDNS,, 

CALI  ASSIGN(3» 'PITCH. TAB* «9)    fOUTPUT  FROM  PITDET. 

CALI  ASSIGN(«I«  'OUTPUT*  «6)        'OUTPUT  OF  THE  COMB 

•FILTER. 

rEFTNE  FILE  ?  ( 251 , £56* U, KVAR ) 

EEFINE  FILE  3  ( 1 000  *  2 « Ut I VAR) 
C  WINDOW  DETERMINATION. 

(AL!  WlrPOW  (NTYPEtK) 
C  CALCULATE  COEFFICIENTS. 

CALI  CUf.  F  (  NTYPF«K.i  A»L) 
C  SPECIFY  THE  LENGTH  uF  Th[  OUTPUI  FILE. 

I  I  "1 T  T  =  253  +  (  K  *  i)     +  1 
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LLFTrt  FILE  u  (LIMlT»256tU«  JVAK) 
C  FILL  THE  PUFFER  INITIAI  L  Y 

t  PEAU  IHl  FMsT  Va!  UL  HRO|v  THE  PITCH  TABLE  FUR 
C  INITIALIZATION  KuPPOSES. 
2  0       \LI\D     (5»IVAR)  M 
C  TOUHLL  THIS  VALUE  TU  COMPENSATE  FOR  TWO  CHAT!, 

(-  =IABS(2**|) 
C  TS  M  TOO  LARGE? 

)F  (!w.GE,tK)0)  00  10  ?0 

r  now  we    "Iay   ^et  up  buffer,    take  •vie:asureivients. 

bO      C=L-j 

I V  =  (  (  C  /  2 )  *  v )  +  1 

lR=TV/256 

IF  ( IR*25fi.NF.IV)  1  (  = I K  +  1 

IP=IR*256*2 

I=(IP  f  1)  -  (K  *  K) 
C  FILL  -Wit    PUFFER  FROM  WHERE  THE    POINTER  STARTS. 

CO  80  II=(2*TR)+l«lfc 

READ(2,KVAR)     (DATBUf(J    +     ((II-l)     *256>>»    J=1,2S6) 
8G  COiMTINUe 

IVAR  =  2  !K[SEI     PUCH.TfB    FILE. 

C    MOW    BEGIN    THE    PROTEASING. 

iop    fflao  =  .false,  ithis  flag  dft.  when  we  change  pitoh. 

2  F  L  A  G    =    0 

REAr)    (  5  •  I  \/  a  R  )     M 

jF     (P'.Ll  ,0)     7FLA0    =    3 
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ly  =  lABS(?*V|) 

C  TS  M  TOO  LARGE? 

IF  (Pt.UT.ifOO)  GO  TO  j.30  iM>i400  COKKES.  TO  PITCH<50. 
C  YFS,  TURN  OFp  pi,  Tf.K. 
120     L  ARGE  =  .TKUF  . 

I?1-"     IF(PFLAG)  bQ  TU  100      !IS  PFLAG  SET? 
C  PFLAG  IS  NOl  SET. 
C  ARE  Wf  rHROUGH  WITH  THE  [UFFER? 

IF  (T.GJ.U09A)  GO  It  127 
I  NO,  GC  ReTU*N  an  OUTPUT  POINT 

ASSIGN  125  TO  KXX 

CO  TO  200 
t  YFS,  kE  ARE  THRU  PUFFER,  READ  JN  A  NEW  BUFFER, 
1?7     CO  126  jj  =  1,16 

IF (KVAK,GE.2K2)  00  10  300 

FEAP(2»KVAR)  (DATBUfcICB  +  ((JJ-1)  *  2S6) ) ,ICB=l»236) 

12P     CONTINUE 
129     ]  =  ] 

GO  TO  125 
C  |w  WAS  NOT  TOO  LARPEi  PROCESS  THE  SIGNAL. 
160     LARGE  =  .FALSE. 
C  TAKE  MEASUREMENTS. 

lV=((C/2)        *    M)    +    3  »IV  =  1/2(LENGTH    OF    THE    wilMDUW) 

lR=IV/256 

IF     (IR*256.NF.IV)     lR=IR    +    1 

lP=Tp*256 
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C  IS  n.l     PTTCH  FLAG  SET? 

l^f     IF(PFLAG)  bo  TO  10(j 

c  ro,  rcrs  the  buhfer  meed  to  be  shiftfd? 

150     If-  ((I  +  I\f)  ,LT,1Q9b)     t.0  TU  160 
C  YFS»  SHJrT  THE  BUFFER 

ihoi  n=i 

hVAPrK^AR- ( ?  *  IK) 

CO  1  5fi  JJ=1  ,56 

IF  (KVAR.Gti.P52)  GO  TO  300      !  ARE  WF  THRU? 

REAn  (2*KVAR)  (DATBUF(J  +  ( ( JJ-1 ) *256) ) « J=l » 256) 
1 ^r  COMTINUE 

C  RESET  THF  POINTER  I. 
159     START  =  (1b  -(2  *  )|))  *  2^6 

J  =  IHOiD  -  START 

GO  TO  190 
C  IS  I  103  SMALL? 
16P     IF(I.GE.IV)  GO  TO  I50 
C  YFS,  SLIDE  MiE  RUFFLR. 
3  70     uHOLD  =  I 

cVEX  r  IV  -  T 

uVEY  =  JVFX/P5b 

IF((JVEY  *  256) .NE.oVEX)  JVEY  =  JVEY  +  1 

KVAn  r  (KVAR  -  16)  -  JVEY 

CO  175  IPy  =  It  16 

REAp(?'KVAK)     (uAlFUH     IDX    +     ( ( IPY-1 ) *256> ) « IDX=1 1 2S6) 
175  COMTIWUE 


194- 


C  RESET.  1 

T  =  JHOLD  +  (256  *  JVEY) 
IVr     ASSIGN  l,4n  ]p  KXX 
C  COMPUTATION  AND  OUTPUT  SECTION 
20 C     1ND  =  I 
21 0     IP ( INP.LE.25F)  GO  1 c  220 

1IM0  =  1IYD  -  P56 

GO    TO    2]0 
22r  IF    (LARGE)     Go    10    26>0  !  LARGE    SET? 

C    MO*     CALCULATE     a^'D    OUTPU1     POINT. 

Y-0.0 

r.o    P^C    MM=1,|     -    1 

IP0IIMT  =  ]     -     (  (K    +    1    -    MM)     *       M) 

Y  =  Y     +     (A(WM)     *     PATtLF  (  IPOI'MT)  ) 
2<+r  COiMTlNUE 

C     IS     THIS    THE     FIRST    OUTPUT     TO    BE    ^KITIFM. 

IF"     (SFLAG.EQ.O)     GO    if     255 
C    YES,     COM«T     a |  l 0 ia'    DRASTIC    STEP    INPUT. 

EC    ?!sq     IU    r    1  ,  (  IND-^)  ,2 

CUIPIIT     (  Iu)     =    20^P 

CUT  PUT  (  IU  +■  1)  =  0 
2^'4     CONTINUE 

SFLAG  =  0  IFOREVEP. 
C  WRITE  1  HT  OUTPUT  POINT  * 
25?    .  OUTPUT  (  IMD)  =  Y 

GO  TO  ^70 
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26  0     IF(7FLAg,EQ.P)  GO  U  265 

OUTPUT! IND)  =  2018 

( 0  TO  ?70 
26?     CUTPUT(IMD)  =  (TUT  *  (UATRUF(I)  -  204SJ)  +  2048 
2  70     CUTpi.lT  (  IN[)  +  1  )  =  c 
C  ARE  WE  READY  TO  DUMP  OUTf-L'T? 

IF( IND.NE.255)  00  To  200 
C  YES » WE  ARE! 

WRITE (4*  JVAR,EI\in=^00)  (0UT^UT(KKK)»KKK  =  l,2  56) 
2>;0      1  =  1  +  1 

IF(PATBUF( 1) .EQ.409£)  PFLAG=J   (CHANGE  PITCH? 

1  =  1  +  1 
2Vn     (0  TO  KXX,  (125, 110) 

C  THIS  SECTION  IS  u?EU  1<^  F  ILL  DA1BUF  WITH  zEKOS  WHEN 
t;  CALLED. 
300     LO  310  It_AX  =  JJ»16        IDtTERI^lNL  HO*<  MANY 

DO  305  JLAX=1 ,25b 

DATpUF(jLAX  +  ( ( IUAx-1 ) *256 ) ) =0        IRECORDS  NEEU 
10  PE 
305     CONTINUE  IAUDFI)  WITH  ZEKOS. 

33  n     CONTINUE 

IF ( I  ARGE )  BQ  TU  129 

GO  TO  159 
C  YES,  UUIT 
100    •  F  NDF3 LE  2 

f  NDFILE  3 
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EIMDFILE    i* 

STOP 

END 
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C 
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c 
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HUf:  SIU/OI+      AUTHOR;  R,  FM7ILR 

ADAPTFU  From  THESIS  BY  V.  SHIELDS 

TITS  |  R  0  n  F?  A u'  USES  ThE"  ME1HOD  OF  SHIELDS  TO  PERFORM 

THE  SEPARATION  OF  A  SHEAKFR  A^D  NOISE  BY  |v;EAnS 

OF  A  DlGTTAL  COMb  FILTER.   I  HIS  PPt'tiR/vM  USES  TUL 

/\TTLMAlTON  METHOD  FOR  UNVOICED  SECTIONS, 

SUBROUTINES  USED! 
WINDOW  (p'!TY^LtK) 
COEF  (NTYPE»K« AtL) 

parameters  used: 
arrays: 

1.  n/mbuf  -  contains  the  input  sample  values. 

?.  a  -  contains  the  coefficients. 

3.   OUTPUT  -  TEMPORARY  Bu^ER  !-  OR  OUTPUT  VALUES, 
VARlAFLEs: 

1.   SFL/\G  -USED  TO  INITIALIZE  OUTPUT  FILE 
SO  THAT  ^TtP  DISCONTINUITY  DOES  MOT 

OCCuR. 

?.   ZFLAG  -  USED  TO  DENOTE  BE,  OF  SENTENCE. 

•  3.   PFLaG  -  USED  TO  r  ENTOE  NEW  PUCH  ptPTOD. 

14.   LARGE  -  USED  TO  DENOTE  UNVOICED  SECTIONS. 
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INITIALIZA II ON 


i  luirPSIoM  A(15) 

INTEGEH  DAl'BUF(t096)  ,  START,  OUT  PUT  (256)  »ZFLAG,SFLAG 

LOGICAL  PFL AG, LARGE 

LATA  DAtBuK , SFLAG  *  OLTPUT/ fU96*0 , 1 , 256*0/ 

r  'AT  A  IV/\R,  J\//»R,KVAR/?,2,2/ 
C  CALL  THE  FOLLOWING  SUBROUTINES  FOR  INITIALISATION, 

CALI  ASSIGN(P, »RKi:siGNAL.UATt ) jOuTPyT  FROM  ADDNS, 

C.Ali   ASSIG|\I(?,  •PITCH. TAB*  ,9)     'OUTPUT  FROM  PITV02. 

CALL  ASSIGN* " ♦ • OUTPLT1 ,6)        ! OUTPUT  OF  THE  COMB 

! FILTER. 

CEFINE  FILE  ?  (  251  »  2?^  »  U  ,  KvAR  ) 

REFINE  FILE"  ?  ( 100  0  ,  2  ,  U»  IVAR) 
C  WINDOU  DETERMINATION. 

CALL  WIND0W< MTfPfc ,K) 
C  CALCULATE  COEFFICIENTS, 

CALI  CDf  F(NTYPE«K«'At|  ) 
C  SPECIFY  THE  LENGTH  UF  TH[.  OUTpUl  FILp. 

L I  "I  I T  =  251  +  ( K  *  2 )  +  1 

EEFINE  FILE  q  (  LIMI-,  ,  256,  U«  JVAR  ) 
C  EJLL  THE  PUFf-ER  INITIALLY 

C  PEAL)  THE  FI«ST  VALUE  HROp  THE  PITCH  TABLE  FOR 
C  INITIALISATION  Hu-PGSES. 
20       h E  An  <  3  t I v A  r )  I i 
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C    DOUBU      THIS    VAL^E    TU    COMf-ENSATE    FOK    TWO    CHANNELS. 

|v  =  I  a  B  S  (  2  *  W  ) 

PHOLO  -  M  'THIS  WILL  RE  uStD  FOR  UNVOICED, 
C  TS  ft    TOO  LAKqE? 

IE  (P'.Oif.  .'400  )  GO  H  20 
C  rOW  WE  V'/\Y  St  T  Up  nUE'EEH,  TAKE  MEASUREMENTS. 
bO      C=L-1 

iv=(  (CV2)  *  r)  +  i 

lR=IV/?56 

IF  (IR*256.NF.1V)  lR  =  IR  +  1 

IF  =  tr*  256*2" 

I  =  (  I P  +  1  )  -  (  K  *  I") 
C  TTLL  THE  -UFFF  R  FROM  WHEI  E  THE  POINTER  STaRTS. 

r  0   p o    II  =  (2*TR)+l«l6 

READ(2»KVAK)     (uATE-ui(J    +     ((IT-])     *?56 )  )  %  J  =  l »  256  ) 
AC  CONTINUE 

IVAR=2  !KfcSET     PITCH. TaB    FILE. 


C  row  BEGIN  THf  PROCESSirG. 

100     FFLAG  =  .False.   !HjS  Fl_AG  OE I .  WHEN  WE  CHANGE 
PITCH. 

7 FLAG  =  0 

fhoi  d  -  ft 

RE  An     (3»IVAR)     ft 
•     IE      flw.LT.O)     PFLAG    =    1 
lu  =  lARS(2*vl) 
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c;  TS  M  too  large? 

IF  (  r.GEI.  4  00  )  M  =  IM  GLD 
C  TAKE  I1  EASUREl^ENTS. 

lV=((C/2)   %  M)  +  1      |IV=1/2(LENGTH  OF  THE  WtNDU' 

lR=TV/256 

IF  ( IR*256.MF.IV)  IR  =  IR  +  1 

IP=IR*256 

c  is  thf  pitch  flag  set? 

140     IF (PFLAG)  ^0  TO  100 

C  rOi  DOES  THE  BUFFER  NEED  TO  BE  SHIFTED? 
lbO     IT   ((I  +  JM)  .  I  .1  .MO^t)  GO  TU  160 
C  YESi  SHIFT  HIE  PUFFtP 

IH0LD=I 

KVAP  =  K\/AR-(2  *  JK> 

FO  150  JJ=3 ,16 

IF  (KVAp.GL.?5<=;)  GO  10  300       'ARE  WE  THRU? 

READ  (2»KVAR)  (DATBUF(J  -f  ( ( JJ-1 ) *256> > t J=l«256) 
15fi     CONTINUE 
G  RESET  I 
159     START=(16-(2*IR)  }  *  ?56 

I=lHOLD-START 

GO  TO  ?0  0 
(   IS  I  I  00  SMALL? 
16  0     IF(  I.GE.Iv)  <~Q  TO  2oO 
C  YFS,  SLIDE  THE  BUFFER. 
170     oHO|.P  =  I 
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uvEx  =  iv  -  t 

wVEY  =  jVpy/?56 

IFUJVEY  *  256)  .NE.jVEX)  JVEY  r  JVEY  +  1 

KVAR  =  (KVAR  -  16)  -  JVEY 

CO  175  JDy  =  lil6 

RLAn(2»KVAR)     <UMFut(     iD>;    +     (  (  IDY-1 )  *256  )  )  *  lDX=l  t 2b6) 
175  rOVTTNUt 

C    RESET     I 

T    =    JHOLD    +     (256    *    JVEY) 
C    COMPUIAI ton    /\rjD    OUTPUT     SECTION 

20  0  I NO    =    1 

210  IP ( TMD.LE .25c)     GO    lC    220 

I(MD    =    I|SID    -    P5o 

GO    TO    210 
C    r'O,     C/sLCULATE     AMD    OUTPUT    POINT. 
2Rp  IF     (ZFLAG.EQ.l)     GO    TO    253 

Y  =  0  .  0 

r  0    ?'|0    MM=1«I.  -1 

IP0I!\|T  =  I    -     (  (K    +    1    -    MM)     *       N|) 

Y  =  r    +     (A<MP1)     *     (U£TLl'F(  IpOINT)     -     2048)> 
?4r              fOMTlNJf- 

c  is  this  the  tirst  output  to  be  written? 

?b3  IF     (SFLAG.EQ.O)     GO    TO    255 

C    YFSi     COM»T     OllqkJ    PRASTlC     STEP    TKANSlENT. 

CO    ?5M     IU    =    5 , ( INP    -2) i     2 

CUTPUT     ( IU)     =    2  04^ 
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curpuT   ( iu  +  i)   =  u 

25i\  CONTINUE 

Sf  LAG  =  0  ! FOREVER 
C  WRITE  THE  OUTPUT  P0IM1 . 
255  IF(7FLAG.EP.D)     OUTPUT  (  IIMD)     =    IFIX(Y    +    .5)     +    20<+fl 

If-  (ZFLAG.EQ.l)  OUTPUT  (INO)  =  2048 
?7o     IF  (DATdUF(I  +  1  >  •  G»E  • ,+  07^ )     OUTPUTCINd  +  1)  =  40°5 

IF     (DATbUF(I    +    1> •L7.4O75)     OUTPUKIND    +    1)     =    0 
C    APE     WE     *F/\DY    TO    DUMP    OUTpUT? 

TF  (  IND.ME  .  *55)     GO    To    230 
C    YF.StWE     ApF! 

URITE< 4» jvAR,Ei,n=MOo)     (0UTPUT(KKK)»KKK=lt256) 

280  1=1    +    1 

IF(nATBUF(  1)  .EQ.4095)     PFLAG=1        iCHAISIGE    PITCH? 

1  =  1    +    1 
2Vr  gO   TO    UiO 

C  THIS  SECTION  ADOS  ?LPUS  TO  DATB^F  WHEN  CALLkD. 
300     CO  310  ILAX=JJ»16 

r 0  305  jLA/ri ,256 

CATBUF'tjLAX  +  (  (  ILAx-l)*25b)  )=2048 
305     COMTINUE 
310     CONTINUE 

GO  TO  159 
G  YES»  GUIT 
HOP     EMDFILE  2 

ENDFILE  3 
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ENDFTLE    4 

sior 

END 


COEF 
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(enter  ) 

1 

1 

Initialization 

PARKS 
MCCLELLAN 


(return) 


FIGURE    A-11 
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T  I T  |  E:   COEF.FUP  74i?06 


C  ThIS  SUBROUTINE  T^KlS  THf  VALuES  THAT  wfcRf  HEAD  RY  WTNDUw 

C  AND  CALCULAIFS  rue  COLFFlClLNTS  FOH  THE  APPROPRIATE 

C  WINDOU  lf-'AT  WAS  REQUESTER,   IT  RETURNS  FHESt  COtFFlCIENIS 

c  to  the    "Iain  program  in  array  "a". 

C"    DEFIMTIONS: 

C  1.        A    -    ArkAy    C0IMTAI|\iI|\jG    CUEFF  If  IL.PJTS. 

2,  ntype  -  the  value  that  corresponds  to  the 

window  OP  filtek  olsirld. 
s.  k  -  the  value  That  determines  the  length 

OF  THE  p.'PUlSf  FvFSPONSL  nF  THE  WINDOW 

i 
ANO  THE  NUMBt-R  OF  COEFFICIENTS  CALCULATED. 

'4.   L  =  2  *  K  +  2 


SUBROUTINE  COEF  ( NTyFE* Kt A* L) 
DIMENSION  A (15) 

r  o  in  i  =  i,i5 

Ad)  =  0.0 

ir    continue 

I ]  =  5. l M 1D9 

SCALE  =  0 . P 

L  =  ?  *  K  +  9 

(0  TO  ( i00» 200*3^0  u+OOt 600* 650) ,  NTYPE 
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C    P/\MUir.  G    WINDOW    SEI  FCTLT 

10  0  CO    110    I    :    1,L-1 

A(I)     =    .5    -     ,5*(C0S( (2*PI*1)/L) ) 

SCAI  F    =    SCALF    +    Ml) 
13  n  CONTINUE 

GO    TO    50  0 
(;    M^JING    WINDOW    SEI.  ECTLD. 

20  C  [0     ?]0     I     r     1,1.-1 

Ml)     =    .54     -     (  .46*  (COS  (  (2*^1*1  )/L)  )  ) 
SCALE    =    SCALF    +    A(  I) 

21  r  {  ONTINUe 

GO    TO    50  0 
C    PLACKPAM    wIMDOw    SELECTED 
300  (0    310     1     =    1,1-1 

Ml)     =    .^2-{.5*<C0S( <2*PI*1)/L> ) )     +     (  .OA*(COS<  (4*P1* 
1    D/L)  )  ) 

SCALE  =  SCALF  +  A(  1  ) 
31o     CONTINUE 

(-0    TO  50C 
("  PECTaNGJLAr  WINDOW  SFLTL1FO. 
4np     CO  «H0  I  =  1,L-1 

A( I)  =  1.0 

SCAI  F  =  SCALF  +  A(  I) 
41 r      CONTINUE 
C  CO^.P-UTE  THE  SCALE  FACTOR, 
50  0     LO  520  I  =  1,15 
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M 1 >  =  A(I) /SCALE 
5?C     CONTINUE 
CO  TO  70  0 

c  ^r:  arbitrary  set  of  coefficients  may  be  Typed  in, 

600     CAL|   ASSIGN  (5t fKKl: » i-1) 

DEFINE  FILE  *  ( 2»25fc»U« NVAN) 

r\VAp  =  i 

REAP  (5iN\/AR)  (A(lY),lY  =  1»?*K+1) 
FNDFILE  5 

GO  TO  70  0 
C  PfPKS-MCCLELLAN  COEFFICIENTS  ARt-  READ  IN  FROM  FILE. 
650     CALL  ASSIGN  ( 5 « • KK1 : « , -1 ) 

hEAn  (  b)  |\!h'n  ,iMFG 

MM  =  (NFlL  +  l)/2 

REAp  (5)  ( A(IY) i 1Y=] ,NN) 

ENDFILE  5 

CO  660  J=  1,MN 

A(NFIL  +  1  -J)  =  A  (J) 
bi   0     CONTINUE 
700     RETURN 

END 
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WINDOW 


(enter  ) 

" 

Initialization 

i 

r 

Accept  filter 
type 

I 


Accept  number 
of  coef  f  icie  nts 


(return) 


F  IGURE    A-12 
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i  Hi  r:      wTMnow.Fup. 


7'4  1??0 


C    THIS    SUBROUTINE    ftl LOWS    FOl'R    POSSIBLE    WINDOWS, 
C    f\    PARKS-WCCLFLLAN    Fll   TER,     OR    (\N    ARB.     FILTER    TO 

c  pr  selected  and  a  value  cf  k  (Proportional  10  the 

C    |  ENGTF    OF    THF     WIISIPOVJ)     TO    TE    SELECTED    FOR    THE     r,AlN 
I"    PROGRAM    TO    USE. 

c 

SUBROUTINE  WINDOW  (|TYPE,K) 
U R  I T  E  (  7  ♦  1  u  ) 
10      fORmaTC  »««TYPE  IN  NUMBER  COKKpSP.  TO  WINDOW*  ♦/»  IX, 

1  iDESIRED:  HAMN=1i  HaFM.=2»  BLACK. =3,  RtCT.=M»  ,/, 

2  »ARplTKARY  =  5,  PAR*- S-MCCLELLAN  =  6  (II)*,/) 
READ  (7,20)  MTYPE 

20      FORmaT(H) 

WRlTEt 7,30 ) 
3n      FORMAT('  ♦ t • THE  VAI-L.E  OF  K=?  (2K  +  1=  LENGTH)') 

R  E  A  n  (7,40)  K 
MO      F0R|"'AT(I2) 

RE  TURN 

E  ND 
C 
C 
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F  I  L  I  N 


(enter  ) 


Initializati  on 


Set    up  array 

S    wi  th  next 
pitch   value 


Compute  f  ilter 
length 


(return) 


FIGURE    A-13 
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C 
C 

c 
c 

C  THIS 
C  PFGIf, 
C  AREA 
C 

c 


Till  r :   FiLiN.i-  ok 


75U315 


IT 

C  I  SUM 


SUBROUTINE  ArTb  10  irillALl7t  1 HE  FILTER  f\l     THE 
|\I\|G  OF  THE  PROGHAM  OR  AT  I  HE  ONSET  op  A  VOICED 

after   af  uimvoiceu  arf  a. 
initialization 

SU3RbUrjNE  FILIN  ( K, iPlT, S» ISUM) 
INTEGEK  S(Jit) 

CO    10     *C     -    1,2+K 

S(<C)     -    IAC5S    (1PIT    h     2) 

COMTINUE 

=  the  length  of  the  filter* 

TSJr  =  U     *  K  *  THll       IRLMFW  ISUM. 

re  turn 

f  I  i  D 
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CALC 


(enter  ) 


Implementation 

of  difference 
equation 


(return) 


FIGURE   A-14 
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C 
C 
C 

r 
c 
c 
c 
c 
c 
c 
c: 
c 
c 
c 
c 
c 
c 
c 


TITLE:   CMC.  FOP  75Q313 

THIS  SUBROUTINE  Ic  USLD  t-lTH  THt  MAIN  PROGRAM  NONVOl 
/\NO  IS  JSEP  TO  CALCULATE  AN  OUTPUT  POINT, 
PARAP  [  TERS: 
ARRAYS : 

1.  A  '  CONTAINS  THL  COEFFICIENTS 

2.  DATbUF  .  CONTAINS  THE  INPUT  DATA 

2.   S  -  CONTaH'mS  THE  SPECINGS  BETWEEN  THE 
COEFFICIENTS, 
VARlAELES: 

1-  t  "  pointed  °r  the  current  processing  in  datbuh 
?•  k  -  Value  that  lfteRMInes  the  no.  of  coef. 
3.  iout  -  a  value  keturneo  to  the  n />  i m  program, 

THF.  OUTPUT  OF  THE  FlLlEH  EpCM  ONE  POINT. 

SUBROUTINE  C.*LC  ( A t DAT6UF« S, I , K, IOuT) 

DIMENSION  A ( 1 5 ) 

INTEGER    DA>'BUF     (M-09fc.)  iS(  14->  tOUl 
COMPUTATION    OF     ' Hf    OUTPUl     POINT 

Y  =    0.0     UNITIALIZE     THL    SU"1    VAUiE. 

Y  -    A(1)     *    DATBUF(I)  !PtRFOK|w    FIRST    MULT. 
IN Dry    =    I 

[0    ?0     IY     =    ?,(2*K)     +    \        IPLRFOHM    KEST    OF    MULT, 
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rNoirx  =   imuev  -  sd\   -i) 

Y  =  Y  +  (  /\(  IV)  +  DATlUF(  iNOty)  ) 
2P      (OMTINUE 

IOJT  =  IP  IX  (Y  +  .CJ)      IRUUMD  THE.  OUTPUT  VaLUL 

hi  TURN 

FMD 
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I  I  II   L;         ADAPI  .HOC 
DOcU^FNTATlUM    i  OR    NONV01     AMP    NONVn? 


AI-'STPACl 


ThFSE  HROGkai*S  J  f 'l|  I  E>|E|\|T  THE  ApfFTjVL  SySTLpoS 
FORMULaTEJ  1M  THIc  THESIS.  MOMVOl  CORRESPONDS  To  I  HE 
ADAPTIVE   OVERLOAD   SYSTF.f   WHILE   NONV02   IS   THE   ADAPTIVE 

SYSTEM. 

SUpROUTlNES  used: 

A.   Mnr-'VOj 


1.  mllNDOW  (NTYPl»K)  -  THIS  SUBROUTINE  ALLOWS  FOP  I  HE 
USpR  TC  PRESCRIBE  'HE  IINDOW  FUNCTION,  HAR«S_MCCLLLLaN 
KILTER,  0*  aN  ArRITKaRy  I  H TEN  DESIRED  IN  THE  PROGRAM  A^P  TO 
DfcSl^NATE  THE  VALUE  OF  K  FoR  1  HE  FILTER  LENGTH, 


2.  CoFF  (NTYpF,K,A»L)  -  THIS  SUBROUTINE  CALCULATFS  I  HE 
COFFFIcIE'^TS  OF  THE  FjlTLK  PRESCRIBED  IN  THE  WINDOW 
SUpROOTlNE  ANO  STURFS  THEM  IN  AN  ARRAY  "A"  THnT  LAN  RE  OF 
PIA>  IMDM  DIMENSION  OF  1?. 


6.  FALL   (  A.nAlPUF  ,S  ,  I,  K.IDDT) 


THIS    SUBROUTINE 
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KtTURNS  CNF  VAlUF  or  t  Ml"  OUTPUT  PtK  CALL.   jj  IMPLEMENTS  I  HE 

AU/iPTIVf  FTI  T  LF  DIFFERENCE  FOUA1ION. 


'+.  FlLiN  (KtMEWPl,StISUM)  -  ThIS  SUpRpUTiNE 
INITIALIZES  1'HE  FTLTEH  AT  THE  BtGINNING  OF  HIE  PROGRAM  ANu 
Al  THE  CNSFT  Op  A  y/nicru  ARFA  FOLLOWING  AN  L'NvOICED  ARt-A. 
THF  Fl|lEp  SpAC^MPS  APf  All  SET  EQUAL  TO  THE  NEW  pIlCH 
PtRlOO  VALUE  IINilTiAl  LY. 

B.   MONV02 


1.        rfTf'pOW     (iMTYpLtK)     -     (SEE    EXPLANATION    ABOVE). 


2.   COFF  (NTYPPt*,A*L)  -  (  SEt  EXPLANATION  ABOVE). 


3.   FTI.  IN  (K,NrWHI«S,  ISUH)  -  (SEE  EXPLANATION  ABOVE). 


H.      PALc?    (1  iPATLir  »Si  I  »  * »  IOl'T,  ISTaRT  )      -  THIS 

SUpROUTlNE  IS  BASICALLY  THE  SAME  AS  THF  CALC  SUBROUTlNf  WITH 

THF  EXCEPTION  THfT  TT  PROVIDES  FOK  THf  OVERLOAD  PROBLtr*  AND 
C( | PECTS  Fpp  IT . 


5.   3( -r,  (pSTARTtS.IiK)  -   THIS   SUBROUTINE   INITIALISES 
tV'F   ISTA^T   ARRAY   WHILE   IS   USED   IN  THE  L'LTECTION  of  ihe 
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OVFPLOAC  PROBLEM. 


6.   RFBEG  (ISTARTflfK)  -  THIS  SUBROUTINE  RE- INITI aLI^FS 
l|T  ISTAR1  APKAY  AFTER  A  BUFFER  SHlFI  hAS  BFEn  pERFORMEO. 

(HlTp  PROGRAMS  USfD  WITH  » |  F  S  L  SYSTEMS 

(SFE  THt  SHItLn.D^C  EXPLANATION;  ALL  PROGRAMS  ARE  THE  SAML.) 

INPUTS  TO  adative  system 

1.  SjRNAt.DAT  -  FROM  ADDNS 

2.  PjTCH.TAB  -  FROM  |  J  TV  02 

OUTPUT    f-ROr    1H[     ADAPTIVE    SySTEM 
1.        OUTPUT.DAT 

Al RORlThM 


i.     the   algokithm  is  jasjcally  the  samf  as  thf  shielus* 

SYSTEM.   ONLY  THF  EVCEPTlOfS  ARE  LISTED  I|\l  THIS  DESCRIPTION. 


2.  IHE  PITCH. JAB  EIL[  Is  IN1TIAU.Y  STORED  IN  AN  ARKAY 
"P"  SO  EASIER  ACCESS  Tc  THt  PITCH  PERlOp  VALUES  CAN  BE 
OBTA INEP. 


3.   AM  MRAY  "S"  TS  USED  TO   HOI  U   THE   SPACING   VALUES 
BtTWEEN   T  MT   r OEFF T C A ENTS .    THE   VARIABLE   ISUM  IS  USEO  TO 
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i  I  I  OIL'  ihE  I  EW&TH  Or  I  UK  f  \\   !  K . 


<+.  IN  MONVO?  AN  AKRA'j  "IStAKT"  IS  uSf  P  TO  HOLD  I  HE 
J I '  T  T I A  l  STARTING  VALUES  \(R  A  PARTICULAR  FILTER  RASED  ON  A 
HITCH  UK  J  Op  VaUJE.  IT  IS  USED  JN  DETECTING  a^1  OVfRLUaD 
CONDITION, 


b.  vhe  output  puint  is  computed  by  thf  difference 

Ef'iAJlcr>       iiSED       FUR        I  HE       APAPTlyt       FILTER,  IN       N0NV02     "HE 

uvrpLO/u  problem  is  cuPKfiirn. 


6.   Tf'F  BUFFER  IS  SHlf TFU  IN  A  SlplLlAR  MANNER   To   I  HE 
SHIELDS'  PPOGK^r. 


7.   DI'RINp  UMVriCED  SECTIONS  THE  INPUT  IS  ATTENUATED  BY 
THE  VAR TABLE i  CO|\|,  TO  OBTAIN  THE  OUTPUT. 


8.   WHFN  THE   LOGICAL   VARIABLE,   DONE,   IS   TRUE,   I  HE 
PROGRAM  PERFOIS  a  SpRiTS  Of  CLOSING  STfPS,  A'MD  THEM  HAL.TS. 


t  r  p 
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C 
C 

c 

c 
r 
c 
c 
c 
c 
c 
( 
c 
c 
c 
c 
c 
c 
c 

3 

( 

c 
c 
c 

L 
C 


TITLE*   NONv^l.rUR 


75U31 1 


THIS  fROGRf'1  IMPLEMENTS  THE  NONUNJFORMLY  SPACED 
ADAPTIVE  FILTER  WTTHOUT  CCC?KECT1^!^;  RoR  1HE  UVLRl  QAD 
PROBLEMS. 

SUBROLTINES  USED: 

3  .   WINDOW  ( IV'TYPE  »K) 

?.   COEF  ( MTYPEtKi AiL) 

3.   FILIN  (K,NEWPI,S,ISUM) 

U,       CALL  ( A,DATBUF,StI«K,IOUT) 

A  R  R  A  y  S  : 

1.  A  -  CONTAINS  COEFFICIENTS. 

2.  HA1RUF  -  INPUT  VALUES. 

3.  OUTPUT  -  TEMPORARY  OUTPUT  ARRAY. 
il.  P    -    USED  TU  STCRF  PITCH  PEHjOUS. 

5.  S  -  USED  TO  STORE  SPACINGS  BETWEEN  COEF. 

VARIABLE? ! 

1,  DOME  .  LOGICAI  l  t  NOTES  END  OF  SENTENCE. 

2.  ATTEN  -  PENOIFS  UNVOICtD  ARFAS. 
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C        3.   J  "  POINTER  1W  I  ITCH  PtRTOP  ARRAY  "P". 

C       i*a   TSUh  -  LENGTH  C^  THE  FILTER. 

C 

C       INITIALIZATION 

CI  PENSION  A  (15) 

INTEGER  DArBHF(^096)  ,  OUTPUT  (256)  ,P(  1000)  ,S(m) 

LOGICAL  DOME*  ATI  EN 

LATA  DATBUK,  OUTPUT  ,  p  ,  SM0  9&*0  »  ^56*0  ,  100  0*0  %  11*0/ 

CATA  IVAR,  J\/AR»KVAR/2,  2»2/ 

CONE  =  .FALSE. 

ATTEM  =  .EAL^F. 

CALL  ASSIGN  (  2  »•  RK1  •  SIGNAL.  DAT*  ) 

CALt  ASSIGN  (3* •RK1;PITCH.1AR» ) 

CALL  ASSIGN  ( if  »♦  RKl :  OUTPUT.  DAT*  ) 

CEFINE  FILE  2  (  2bl  «£56  .  U,  I  v  AR ) 

CEFlNE  FILE  3  U  000  ,  2  »  Ui  JV/»R) 

CO^  =  .3   I  ATTENUATION  CONSTANT 
C  CHOOSE  WINDOW  TYpr 

CALL     WINDOW    (NTYPEtK) 
C    CALCULATE    THE    C^EEFIC iL hi S 

(A  LI      COeE     (1MTYPE«K»A»L) 
C    SPECIFY    THE    |EmGTH    Uf    THf     OUTPU1     FILE. 

i l^TT    =    251    +     IK*?)     +    1 

[EFTNE    FILE    "     (  L  1  iv  1  I  »  ?  d6  1  U  *  K  V  AH  ) 
C    FILL    THE    PITCH    ARPAY. 

CO    100     I A    =    1,100  0 
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REAp  (  3»  JvAR,  END  =  1 1  n  )  P(JA) 

IPO     CONTINUE 

llr     COMTIN'JE 

u  =  1  !J  IS  THE  HOIMER  IN  THF  PITCH  PtROD  ARRAy, 
C  THIS  PART  DETERMINES  IT  i,F  HAVE  A  SILENT  AREA  AT 
C  THE  [ [ GIHNI MG  Of  THE  F ILL. 
12C     IF  (P(U).GE.O)  GO  Tc  200  !  SILENT  AREA? 

I"  =  IABS(P(J)  *  2)  !YES,  SLT  THE  OUTPUT  =  0. 

CO  150  KD  =  J  ,  <M-1)  ,2 

KA  r  KD 
130     IF  (KA.LE.ii55)  GO  To  140 

hA  =  KA  -  ?5^ 

(-C  TO  130 
l«*p      CUT  PUT  (KA)  r  20H6 

CUTpUKKA  +  ])  =  0 

IF  (KA.NE.ii55)G0  TC  150  j  IS  IT  TIME  TO  WRI1F  A  pEC? 

UR1TF(  ttKVAR)  (OUTPLiT(KB)  «KB=1«256)    !YES. 
150     CONTINUE 
C  PUMP  THE  COUNTER. 

c   =   J   +   1 

(0  TO  1?0 

c  pitci-    PtPioo  was  iuoi   a  Silent  one. 

200     IF(P< J) .GE»400)  GO  10  50  0  0 

C  VOICEC*  srT  yP  f'HF  FILTEf.  SPACING. 

ASSIGN  350  TP  KXX 

CALI   FILIN  (KtP(U) iSi I  SUM) 
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C  n  ASURE^FK'TS  SECTION 

1R  r  I SUM  /2C6 

IF  ((IK  *  256).NL.ISUM)  IR  =  IK  +  1 
G  SFT  UP  FOR  INITIAL  KFADING  INTO  BUFFER, 

IV AP  =  KVAK  -  XR 

CO  ^00  KE  =  1  ,16 

READ(2»IVAR)     (UATBUF(KF    +     (  (KE-1)*256)  )  tKF    =    1«?56) 
30  p  f.OMTlNUE 

c  sft  uf  iMF  Buffer  pointer 

IF  (KA.EQ.25?)  KA  =  -1    } RLSEI  THE  POINTER 

! TU  THe  LAST  OUTPUT. 

I  =  (IK  *  ?5£)  +  <Ka  +  2) 
C  PC  THE  DUTPJT  CALCULATIOKS 
350     CALL  CALC  ( A,DATBUF , S» I t K« *OUT ) 
C  COMPUTE  THL  ADDRESS  FUR  IFF  OUTPUT  POINI. 
36C     IND    =  1 
365     IF  (IND.LE.255)  g0  1  °  370  ! DO  WF  HAVE  A  GOOD  ADDRESS 

IUD  =  InD  -  P56  |MO  Go  KECHECK, 

C-0  TO  3&5 
37p     iF(ATTErj)  GO  TU  57b  !IS  THIS  A  VOICED  SEGMENT? 

CUTPUT(IND)  .=  10UT  | YES*  GtT  CAl CULATED  VALUE. 

PC  TO  3p0 
C  ATTENl  ATP  THE  INPUT. 

375     CU1  PUT  (  jrjrj)  =(CON*(  DflTBUF  (  I )  -2  0*+8 )  )  +20^8 
360   •  CUTpijT(  IIMD+1)  -  0        [SECOND  CHANNEL  =  0. 
P  IS  IT  TIME  U  DUvip  I  HE  OUTPUT  BJFFt  H? 
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IF     (IN0.NE.255)     GO    to    400 

URITEC*»KVAR,EI\ID=3000)     ( OUTPUT ( KB) » KB=1 »256)     !YrS. 
40  0  1=1+1 

C"    TS    THE  RE    A    PITCH    I^RK? 

IF     (DATBUF( I) •LT.4075)     GO     10    1J>?0 
bnp  f ,  E  ^  p  T    =    P(J)     *    ?     !YESi     GET    A    MEW    VALlJF    OF    PITCH. 

c    =  J  *  3  IBUMP  THE  PITCH  PERIOD  POINTER. 
C  ARE  K-  ENTERING  A  SILENT  AREA? 

)F  (NEWPl.l T.O)  GO  TO  1600 

c  mo,  AhE  we  entering  aim  unvoiced  area? 

IF  (NEWpi.GT.fOO)  &C  TO  1BO0 
C  NO,  IhE^rrOKf  we  rRL  ENTERING  A  VOICED  ARFA. 
C  CHANGE  the  filter. 

ASSIGN  35n  T^  KXX 
C  TS  THE  FILTER  BEING  INITIALIZED  AFH R  AN  UNVOICED  AREA? 

IF  (ATTEN)  GO  TO  80o 

GO  TO  900  !THE  FILTER  IS  NUT  BEING  RE- INIT 1 AL IZED. 
BOO     ATlFT1  =  .FALSE.   ICLEAR  THl  FLAG. 

(ALL  FILIM  (KtNEWPIvStISUm 

GO  TO  1220 
90 C     hAA  r  2    *  K  • 
IOP0    IF  (KAA.EQ.l.)  CO    10  jlOO 

S(KAA)  =  S(KAA  -  3)  'SHIFT  THE  SPACI|\!GS 

^AA  r  KAA  -1 

•  go  to  icnn 

npn    Ml)  =  NEWPI  \*{'l    TfF  NEW  PITCH  PERIOD  SPACING 
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•    ISJ!"!    =    0     ICO^TUTL    ThE     LEfviGIH    OF     THE    FILTER 

10    1?0n    Kl    -    1«2*K 

IS  Jr.    =    JSUM    +    S(K1  ) 

yypv       continue 
i??r       1=1+1 

(  JS  II  TIME  10  SHIFT  THE  BUFFER? 

IF  (  T.LT  .40  9?  )  PO  T(  140  0 
("  YE  Si  SHIFT  I  PL  BUFFER. 

IR  =  I  SUM  /  J>56 

IF  ( (1**256) .NE.ISUp)  IR  =  IR  +  1 

IVAR  -  IVAr  _  IR 

CO  130  0  k?  =  It  16 

IF  (IVAR.GE.P52)  60  TO  290 U 

READ  (2»IVAR)  (DATBUF(K3  +  ( ( K2-1 ) *25& » > »K3  =  1,25b) 

1300    CONTINUE 

(  RFSEl  I 

1310    I  =  (I«  *  25b)  +  3 

1400    f:0  TO  KXX,  (  350  i  360  , 1 840  ) 

C  THIS  /iREA  RESERVEF  FOR  SILENT  A*EA  LATER. 

If  00    KEWPI  =  IABS(NLHHI)  IGET  RID  Of-  MINUS  SIGN 

COMF  =  .TRUE.  !St-T  DONE  FLAG. 

fO  TO  1800 
C  UNVOICE3  ARE/5 
18P0    1=1+3  IBLMP  THE  POINTER. 

IF  (NEWpI.LE.ISUM)  GO  10  20  On 
1P20    K7  =  1  !Sfc7  UP  A  COUNTER. 
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1 P  ?  5         K  A  6    =    0 

c  /\RE  we    through  shifting? 

1830    IF  (KAB.EQ, ( (2*K>-K7) )G0  TO  1835 
KNEW  =  2    *    K  -  KAB 

i<o^p  =  \<Nru   -  i 

S(KMEW)  -    S(KOLP)         I NU» CONTINUE 

KAB  =  KAB  +  1 

GO  TO  l83n 
1P35    K8  =  1 
1837    ASSIGN  10-4  0  to  KXX   . 

GO    TO    5b0     !GO    COMPUTE     AIM    OUTPUT, 
1840  Kfi     =    K8    +     2 

IF  (K8.lT.Sm)  GC  10  1837      !  ARE  WE  THRU? 

K7  =  K7  +  1  !YtS,  BUMP  K7 

C  ENTER  TUB  ATTENUATION  PROCEDURE 
18^0    ATTpN  =  .TRUF.  !SET  THE  ATTfNUATIUN  FLAG. 

IF  (POMr)  GO  TO  ISOC      !A*E  WE  ALMOST  THROUGH? 

ASSIGN  360  TO  KXX  ![Vr,  THIS  IS  ONLY  jsU    UNVOICED  ARLA. 

GO  TO  3f0 
190  0    I H  D  =  I 

C    ro    Wf     HAVE    A    SiJii.'vRLr    INl? 
1LI?0  IF     (  IND.LE.2^5)     GO    1 0    195  0 

JND  =  IND  -  ?5b  !NO«  TRY  AGAIN. 

GO  TO  1920 
19«S0  .  CUTPUT(IND)  =  ^OHP  ! SILENT  AREA,  THEREFORE 

! SLT  OUTPUT  =  2048  (OR  o)  . 
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CUTPUT  (  IMP  4  1)  -  i. 

II   (  IN0,NE!.255)  GO  10  1970  !Tir"E  ro  DUMP  OUTPUT  BUH? 

URlTFlC+'KVARiENn  =  3000)  (  OUTPUT*  KB)  ♦  KB  =  1 ♦  256) 

INO  =  -3  IRtSEI  POINTER 

1970    IIMD  =  IMP  +  ? 

GO  TO  l92n         ICCMT.  UMUI  END  OF  FILE. 
("  TPIS  t\hl(\    WILL  HE  DEALING  WITH  UNVOICED  SEGMENTS  WITH 
C  ISUM>|\EaIPI. 
2000    GO  TO  3P0  0 
29Q0    CO  P9I0  Il/^X  =  K2»lfc 

I  0  ?90b  JL^X  =  3.  »25< 

LATRUF(JLAX  +  ( ( ILAx-l)*25b)  )  =  2048 
2Vp5    COMTINUe 
29^0    COMTINUE 

GO  TO  1 3 1 n 
30  0  0    EIMDFILE  2 

tkOFTLE  3 

FNDFILE  4 

STDP 

END 


N0NV02 
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1  n  i  t i a  li  zat  ion 

V 

CALL  WINDOW 

'! 

CALL   COEF 

' 

> 

Fill    pitch 
array 

yes 


Fill    out  put 
file  with  zeros 


}\_yes 


Assign  350  to 
KXX 


Qui  t 


CAL  L  FILIN 


300 


Fill    buffer 
i  nitiul ly 


i 


CALL     BEG 


FIGURE  A-15 
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N0NV02 

350 Jr_ 


CALL  CALC2 
Result— >IOUT 


360 


Prepare    for 
output  address 


37  0 


no 

s'  U  n v  o  i  c  ed^x^^ 

yes 

T 

y 

' 

\.    S  eg  m  e  n  t  ?/>~ 

"i 

OUTPUT=IOUT 

OUTPUT  .CON  # 
1  NPUT 

Output   for    2nd 
channel  =  2048 


yes 


Write   output 
onto    disk 


End 


400 


no 


GO  TO      1220 


Quit 


yes 


1600        j 


1800 


Assign 

350    to 

KXX 

© 


Set    DONE 
FLAG 


Move     filter 
out  of  voiced 
area 


<D 


FIGURE  A -15  (CON'T.) 
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N0NV02 


no 


900 


1850 


Shi  f  t  spucings 


Set    ATTEN 


Clear  AT  TEN 


i 


yes 


SI  ide  but  f er 


hcl 


?     ~>  'IV 

_ r. 


Assign    360  to 
KXX 


Reset     I 


t   1 


920 


GO  TO     360 


Fill    output  file 
with    zeros 


820 


CALL  F  I  LIN 


Quit 


CALL   BEG 


1220 


Bump      I 


Shift   buffer 
Reset      I 

+ 

CALL 
REBEG 

1400 


GO    TO    KXX 
(350,360,1840) 


FIGURE  A-15(CON'T.) 


-230- 


c 
c 
c 
c 
c 

L 
C. 
.C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

L 

c 
c 
c 
c 
c 
c: 


t i t i  e:      NoiMvnp.ruK 


7  50  32  't 


THIS  ffOrPAl  IMPLFMENTS  THE  NONUNlf-ORPiLY  SPACED 
ADAPTIVE  FILTER  WITH  JUL  CORRECUON  FOR  THE  OVERLOAD 
PROBLFM. 

SUBROUTINES  USED: 

1.  WINDOW  (NTYPEiK) 

P.  COEF  ( MTYPL*K« A,L) 

3.  FILIN  (K,NEWPIiS»ISUM) 

H  .  P  F  0  (PSTART*S,I,K) 

5.  CALC2  (A,DATBUF,SiI»K»10UTf ISTART) 

fc.  REBEG  ( ISTART* I,K) 


arrays: 

1  . 

7  . 
3. 


A  -  CONTAIMS  COEFFICIENTS. 

p  -  used  for  rescaled  loef.  during  overload. 

DATrUF  -  CONTAINS  lNPLM . 
H.        OUTPUT  -  TLrPOK/PY  OUTPUT  APRAY. 

e.  p  -  contains  pitch  periods. 

£.       S  -  CONATINS  SPACING  FUR  FILTER. 


VARIAPLES  : 
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c 

1 

c 

? 

c 

3 

c 

M 

c 

C 

I 

DOMF  -   LOGICAL  FFNOTFS  END  OF  SEN1ENCL. 
ATTEM  -  LOGICAL  DENOTES  UNVOICED  AKFAS. 
J  -  POINTER  IN  THE  PITCH  I'tRIOD  A?*HaY  "P". 
ISUM  -  LFNGTH  Op  THE  FILTER, 

INITIALIZATION 

[ IMFNSION  A ( 1 5 ) .  B(15) 

JNTEGEK  DA  I  BI'F(4UCH, )  tOUTRUI  (?5fe)  ,P(100n)  ,S(  14) 

INTEGER  1ST APT  (15) 

L  OGJCAL  DONE,  ATI  EN 

DATA  DATBUK,pUTPUT»FiS/f0  9S*0»256*0,lOOO*Oil«f*0/ 

DATA  I  START,  6/5  5*0,15*0.0/ 

CATA  IVAR, JVAR»KVAR/2»2»2/ 

FATA  DO^F,  ATTEN/. FALSE., .FALSE,/ 

CALl  ASSIGN  (2, ♦ RKl 'SIGNAL. DAT* ) 

CALF  ASSIGN  (  3  *  •  KKl :  Pi  TCH.  I'  AB»  ) 

CALL  ASSIGN  ( if ,' RK1 : OUTPUT. DAT* ) 

DEFINE  FILE  ?    (  251  »256»  U,  I VAR  ) 

LEFINE  FILL  3  ( 1 U00  ,  ? , U , JV ar ) 

COM  r  .3 
C  CHOOSE  WINDOW  TTpr 

CALL  WINDOW  f Ml YPt  »h  ) 
C  CALCULATE  THF  COEFFICIENTS 

(AH  CUFF  (NTYPEiK,A,L) 
C  SPECIFY  THF  LENGTH  UF  TH|_  OUTPUI  FILE. 

L  1*1  IT  =  25-t  +  (K*2)  -f  1 
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CEFINE  FILL  '!  ( LIMIT* 256» U» KVAR) 

C  F ILL  THE.  PI rcH  ARRAY. 

DO  10  0  1 A  =  1,10  0  0 

HE  An  ( 3* JVARtEMD  =  110)  P(IA) 
inn     CONTINUE 
110     CONTINUE 

d  =  3       |J  IS  THE  P01NTEK  IN  THE  PITCH  PLROD  ARRAY. 
C  THIS  fAKT  DETERMINES  IF  t.F  HAVE!  A  SILENT  AREA  Al 
C  THE  BEGINNING  Of  THE  F I L  (  . 
120     IF  (P(J).GE.n)  GO  Tc  20.0         'SILENT  AREA? 

T  =  IABS(P(J)  *  2)  !YES,  SET  THF  OUTPUT  =  0. 

CO  1  50  KD  =  i  ,  (M-l  )  ,  ? 

KA  =  KO 
130     IF  (KA.LE .255)  R0  U  1^0 

KA  =  KA  -  25ft 

GO  TO  130 

mo    output  <KA)  =  ^ol^n 

OUTPUT (KA  +  1  )  r  0 

IF  (K.A.ME.25S)  GO  Tc  150  US  II  TIME  TO  WRITE  A  RfL,? 

URITEC+'KVAR)  I OUTPuTt KB) » KB=1 t 256)  JYLS. 
150     CONTINUE 
("  PUMP  THE  COUNTER. 

C   =   J   +   1 

(-  0  TO  1?0 
C  PITCH  PERIOD  WAS  ^0 1  A  Si LENT  ONE. 
20  0     lFI'MJl.C-f.qno)  GO  t  C  30  0  0 
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c  voicn  »  set  up  i'hf  filtli-  spacing. 

ASSIGN  350  TO  KXX 

(.ALL  FILIN  (K,P(  J)  »Si  ISUM) 
L  MEASUREMENTS  SECTTOM 

IR  =  I  SLIM  /2r^G 

If-     ((IK    *    256)  .NE.ISl'M)     IR    =    IH    +    1 
C    SET    UF     FOR     J  ri  IT  I  At     READING     IN10    BUFFER. 

IVAR     =    KVAK     -     IP 

[0    30  0    K E    =    1,16 

RtAn<2f  IVAR)     (DATBUf:  (KF.    +     U KE-1 ) *256) > » KF    =    1«256> 

30o    co-mtinue 

C  PPT  UF   H'F  BUFFER  PUTNTEf, 

IF  (KA.EG.25B)  KA  =  -1   iRLSPl  THE  POINTER 

!0F  THE  LAST  OUTPUT. 

I  =  (IK  *  256)  +  (h„  +  2) 
C  INITIALIZE  rhE  ISTART  ARRAY 

CAM  BEG  r 1START«S,I,K) 
C  PC  ThL  OUTPUT  CALrULATlOf.S 

350     CALL  CALC?  (  A  ♦  DATBUf.  ,  S  »  1 1  K  ♦  IOUT  ,  1ST  ART  ) 
C  COMPUTE  THE  ADDRESS  FUR  THE  OUTPUT  H0IN1. 
3C-r     IND  =  I 
3'6fi     IF  ( IND.Lp.255)  GO  TO  370  !  DO  Wp  HAVp  A  GOOD  ADDRESS 

IND  =  IND  -  P56  JISIO  GO  RECHECK. 

(-0    TO  365 
370     JF(ATTEM)  GO  TO  37b  J  IS  THIS  A  VOICED  SEGMENT? 

CUTPUT(IND)  =  IOUT        !YLS,  OFT  CALCULATED  VAlUE. 
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GO  TO  360 
C  ATTENl  AIT  THE  INPMT. 

375     CUTPUT(IND)  =  (COTJ  *  (  DATBiJF  (  I  )  -20'I0  )  )  +  2018 
380     CUTpUK  IND+1)  =  0        (SECOND  CHANNEL  =  0, 
C  TS  IT  TIME  TO  DUMP  I  ML  Ol.TPUT  BUFFER? 

IF  ( IMD.NE.255)  GO  1  0  400 

UR1TE(^»KVAR,END=3000)  ( OU I  PUT ( KB) t KB  =  1 «256)  !YFS. 
400     1  =  1  +  1 
C  TS  THERE  A  ^ITCH  MARK? 

IF     (DATBUF(  I)  ,LT.«+075)     GO     10    1220 
brio  KEa'PI     =    P(J)     *    2     !YtSi     GET     A    NEW    VALElF    OF    PITCH, 

u  =  J  +  1  I  BUMP  THE  PITCH  F'ERlOp  POINTER. 
C  ARE  WE  ENTERING  A  SILENT  AREA? 

IF  (NEWPI.LT.O)  GO  TO  1&00 
C  fOi  ARE  WE  ENTERING  Aim  Uj  VOICED  AREA? 

IF  (NEWpI.GT.4-00)  Gc  TO  1800 
C  NO »  TFEREFOKE  WE  *Rt  ENTfcPING  A  VniffD  ARfA. 
(  CHANGE  THE  FILTER. 

ASSIGN  350  TO  KXX 
C  TS  TEE  FjLTER  Rt'lMC  INITIALIZED  AFTER  AN  UNVOICED  AKFA? 

IF  (AT TEN)  G^  TO  800 

00  TO  9QC  ! THE  FILT[P  IS  NUT  BEING  Rp- J  NIT  1 AL IZfU, 
800     ATTFM  =  .FAL9E.   IClJAK  TH£  FLAG. 
C  TIME  TO  SLIDE  THE  B^FFER-p 

IF  ( I.GT. (NEWPl*2*K) )  GO  TU  820 

cHOLD    =    I  !YeSi     SA.Vt    THe     P01NTEK, 
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iHAr    =    NEWPI    *    2    *    y 

IRA    =     (  IHAH    -     l)/2bf 

IF     ((IRA*256).Nr.(IhAP-I))     IRA    =     IK A    *     1 

IVAR    =    IVAP    -    16    -    IP A 

LO    810     IDy    =    1»16 

REAn     (2»IVAR0     (PATllFdDX    +     (  (  ipY-1  )  *  2b6  )  )  «  T  D  X  =  1 ♦  2b6  ) 
I'ln  (  OMTlNJfc- 

C~  RESET    Trip    POlNTE  R 

]     =    JHOI D    +     (256    *     jPA) 
^^0  (ALL    FILIM    (K»     (NFWpI/2),     S,ISUM) 

CALL  BEG  ( ISTARTtSi l .K) 

GO  TO  1220 
VOP     KAA  =  2  *  K 
1000    IF(KAA.EO.I)  GO  10  H00 

S(KAA)  =  S(KAA  -  1)      ISHIFT  THE  SPACINGS 

KAA  -  KA/\  _x 

GO  TO  1 0  0  0 
3100    S<1)  =  NEWPI  jSEl  ltf  NEW  KITCH  PERIOD  SPACING 

isjm  =  o  i  compute-  The  length  or  the  filter 

CO    1P00    k 1    =    1 »  2 * K 

is Jm  =  isum  +  s(ki ) 

1^0  0  fOMTTNJE 

CALL    BEG     ( 1START*S» 1+ltK) 
1220  1=1+1 

C     IS     IT    TIF"F     To    SHIFT    THE    [  OFFER? 

IF      (  I.Ll  .11096)     GO    H     lfOO 


-236- 
C  YES.  SHIFT  I  HE  PUFFER. 

IR  r  I  SUM  /  ?56 

IF  (  (IK*256).N£.ISU|*)  IR  =  IP  +  1 

IV  AT  =  IV  AH  -  IR 

r 0  130  0  K  ?  =  1,16 

IF  (lVAR.Gt.?52)  GO  10  2900 

I, LAp  (2«IVAR)  (DATbuF(K3  +•  ( ( K2-1 ) *25&> ) . K3=l .256) 
150  0    CONTINUE 
C  PESEl  I 

13]0    I=(IR*256)+3 

C  FFSEI  THE  ISTARI  ARKAY. 

CALL  REBEG  (ISTART.I.K) 
1HO0    GO  TO  KXX,  (350.360.184  0) 
C  THIS  AREA  RESERVED  FOR  SiLENT  A«EA  LATER. 
1600    fxE/JRl  =  iAtJS(NEWHl)       JGET  R  Ip  OF  plNUS  SIGN 

pCVp  =  .TRUE.  !SET  POPE  FLAG. 

(-0  TO  If  0  0 
C  UNVOICED  ARE/\ 

1=1+1  IBbPP  THE  POINTER. 

IF  (NEWpl.LE.ISUM)  If  10  2Ufin 
16?0    K7  =  1  !S[T  UP  A  COUNTER. 

18?5    KA3  =  0 
C  APE  WE  THROUGH  SHIFTING? 
1830    IF  (KA8.EQ.  (  (2*K)-K7)  )  GO  JO  1835 

KNEU  =  2    *  K  -  KAP 

hCLO  =  KNEW  -  i 
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S(KNEW)  =  S(KCLP)         !  NU,  CONTINUE 

KAB  r  K/>B  f  1 

ro  to  i83n 

1P-35    1SJM  =  C 

CO  1836  1TE'  =  li2*K+3 
iSUM  =  ]SU^  +  S(  ITL  ) 

1836  (  OMTI.NUE 

(ALI  BEG  ( ISTAKTtSt l+l«K) 
K8  =  1 

1837  ASSIGN  ]8'40  TO  KXX 

GO  TO  3^0  !GO  COMPUTE  AN  OUTPUT, 
J  8  U  0    K  8  =  K  8  +  2 

IF  (K8.LT.S(i>)  00  TO  1837      !ARE  WE  THRU? 

K7  =  K7  +  l       !Y|_F«  BUMP  Yl 

IF  (K7  -  (2*K))  1825,1835»1850  ! AkE  wE  THRU? 
C  ENTER  THp  ATTENUATION  PRCCEDURE 
1850    ATTEN  =  .TPUT.    !Se.T  1  HE  ATTENUATION  FLAG. 

IF  (DONE)  GO  TO  190c      !A-<E  WE  ALMOST  THROUGH? 

ASSIGN  3f>0  TO  KXX  Iff,  this  is  only  an  unvoiced  apla. 

(-0  TO  3£,0 

i ^ro   imo  =  i 

c  ro  wf  have  0  suitable  inl? 

1920    IF  (IN0.LE.255)  GO  TO  1950 

1N0  r  I  r  Li  -  ?56  !NO,  TRY  AGAIN. 

GO  TO  1920 

1950    OUTPUT  (JNO)  =  2048        ISILEN1  AKEA,  IHEREFORE 
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!  St.  T  OUTPUT  r  2  018  (OR  0)  . 
CUTPUT  ( IMD  +1)  =  C 

IF  (  IND,Nf:.255)  GO  -|  (  1970  »TIP"E  W    pUMP  OUTPUT  BUh? 
URlTECMKvARtEND  =  ^000)  (  UUTPUTC  KB)  ,  KB=1 ,  256) 
JND  =  -l  IRE.SE1  POIMTFK 

1970    I NO  =  J  PC  +  ? 

TO  TO  1^20         ItCTT.  UNI  IP  END  OF  FILE. 
C  TITS  £RE/\  WILL  BE  DEALING  WITH  UNVOICED  SEGMENTS  WITH 

c  isuM>r .Ea/pt. 

20  0  0    URITE (7,2200) 

22P0    FORMAT  (•  »,»pKPGFAf'  FXITEU  FROM  2000') 

(-C  TO  50  0  0 
29Q0    TO  2910  ILAX  =  K2«l6 

LO  2905  JLAx  =  1*256 

CAlpUF  (  JLM  +  (  (ILAX-D*25b)  )  =  2048 
2905    CONTINUE 

29io   co.mtinue 

(-L    TO  13m 
30  0  0    ENDFILE  2 
FNDFILE  3 
ENDFILE  4 
STOP 
END 


BEG 
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(enter  ) 

' 

r 

Initiali  zation 

i 

r 

Set   PSTART 
array  to  zero 

i 

' 

Set    PSTART 
array  to  the 
initial  positions 

(re_ 

ruRN) 

FIGURE    A-16 
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TlTir:   BEG. FOR   75q323 


C 

c 

c: 
r: 

C    THIS    fROgRA1"!     IS    a    SUBROUTINE    pOK    THt     l\!ON\/02 .  FOR    PROGRAM. 

C    Tj     If  T 1  I A  |  I  ^  E  S     I  HF    IS'IAKI     AKRAY    WHLn    CALLED, 

C 

C        INITIALIZATION 

SUBROUTINE  BTG  (  PST/  I'T  »  S  ,  I  <  K  ) 

INTEGER  PSI'APT(15)»  S(lf) 
C  7fRO  THE  PSTaRT  A^Rmy 

CO  in  I U  =  l ,  2  *  K  +  -j 

I-  ST  APT  <  III)  =  0 
in      CONTINUE 
C  INITIALIZE  THE  ARPA*  WITh  THE  SPaCINGS. 

KkQ    -    I 

F  S  T  AR  T  ( 1 )  =  I 

LC  ?0  J  =  *% (2*K)  +  1 

KND  =  KND  -  S(d  -  1) 

PSTART  (J)  =  Ki'MO 
20      CONTINUE 

RETURN 

tND 
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REBEG 
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C  THIS  SUBROUTINE  RF-I^ITIaL IZES  I  HF.  ISTART  ARRAY  AFTER 

C    /••     SHUT     IN     1HE    BUFFi-R    HAS    BEE|\|    PERFORMED, 

( 

C       INITIALIZATION 

SUBROUTIME  RFPtF  (IsTART«ItK) 

INTEGER  ISTARK15) 
('  PERFORM  THE  EXCHANGE  OF  UE  ELEMENTS  IN  ISTART. 

TO  ?  0  J  =  1,  ?  *  K  +  l 

ISTART        (J)     =    I    -     11,096    -     1STAFT(J)) 
20  (.OMTINUF 

RETURN 

F  N  D 
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,r 

Initialization 
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TKIS  SUBROUTINE  IS  USLD  l-ITH  THE  P" /\  I M  PROGRAM  NONVO? 
AND  is  JSED  TO  CALCULATE  AN  OUTPUT  POIIM1.   IT  ALSO 
DETERflWrS  lp  THE  COEFFICIENTS  MEED  TO  BE  "TURNED  OFF". 

parameters: 

arrays: 

1.  A  -  CONTAINS  THt  COEFFICIENTS 

2.  DATbUF  -  CONTAINS  1  'HE  INPUT  DATA 

?.   S  -  CONTAINS  THfc  SPACINGS  BETWEEN  THE 

COEFFICIENTS. 
M.   ISTaRT  -  CONTAINS  THE  INITIAL  LOCATIONS 

of  the  roef.  f  cp  the  present  filter, 
variables: 

1.   i  -  pointer  of  the  current  processing  in  datbuf 

?.  K  -  value  that  letermines  the  no.  of  coef. 

3.  TOUT  -  A  VALUE  RETURNEU  TO  THE  MAIN  PROGRAM, 

THE  OUTPUT  OF  IFF  FILIER  FROM  ONE  POINT. 

SUBPOUTINE  CALC2  ( f  , r AT BUF * S ♦ J < K t TOUT » 1START ) 
1 1 «ENS ION  A  (15)  «  B(]M 

INTEGER  DATBUF  (HOSf  )  t  S{  l«f>  i  ISTARTC15) 
LOGICAL  SCALE 
•  SCAI  E  =  .FALCE. 
SI  E  IF  amy  BAD  SAMPLES  h/VE  BEE'M  REACHED. 


-245- 


^ND  =  I 

h(  1  )  =  /\(  l  ) 

LO  in  IF  =  2« (2*K)  +  1 

JMD  =  JND  -  S( IP-1) 
C  TC  Wf  HAyE  n  COt-FFlCItNT  THAT  HAS  PASSED  JTS 
c:  AREAS? 

II   ( JND.GE.ISTART(IF-I) )  GO    TC  p 

F  (  1 P )  =  A  (  I  P  ) 

r.c   to   lc 
e  F(IP)    =   0.0  !YESi    turiv>   off    the    cof.f. 

sen  f    =    .TRLIT. 
IP  comtinue 

IF  (.NOT.SCALE)  GO  IP  18  ! UO  WE  MEED  TO  RESCALE? 

SCALE  =  .FALSE,   !YLS,  RE  SCALE. 

SCAt  =0.0 
C  COMPUTE  THE  RESCAl.E  FACTf,R 

TO  15  IT  =  1, (2*K)  +  1 

SCAL  =  SCAl  +  b(  n ) 
lb      (OMTTiVJf 
C  RESCALE  THE  cOE»" FTEClENTs 

(0  17  IU  =  1,  (2*K)  +  1 

Fdu)  =  Bdin  /  scal 

17      COMTINUE 

C  rCMPUTATTOM  OF  The  OUTPUT  POINT 
IP      Y  =  0  •  0 

Y  =  P( 1 )  *  DATUUF ( I) 
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IK' DFX  =  I 

{  0  ?0  IY  =2,  (2*K)  +  i 

INOFX  =  INDEX  "  s<  3>  -1  ) 

Y  -  Y  *  (B( IY)*  DATf  IF( INDLX)  ) 

CONTINUE 

I  OUT  =  IFIX  (Y  +  .5) 

RETURN 

END 
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